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PREFACE 


I  gave  the  lecture  series  "Boundary— Layer  Theory"  in  the  winter 
semester  l^lA^  for  the  members  of  my  Institute  and  for  a  considerable 
number  of  collaborators  from  the  Hermann  Goring  Institute  for  Aviation 
Research.  The  series  embraced  a  total  of  sixteen  two-hour  lectures. 

The  aim  of  the  lecture  series  was  to  give  a  survey  of  the  more 
recent  results  of  the  theory  of  viscous  fluids  as  far  as  they  are  of 
importance  for  actual  applications.  Naturally  the  theory  of  the  boundary 
of  frictional  layer  takes  up  the  greatest  part.  In  view  of  the  great 
volume  of  material,  a  complete  treatment  was  out  of  the  question. 

However,  I  took  pains  to  make  concepts  everywhere  stand  out  clearly. 
Moreover,  several  important  typical  examples  were  treated  in  detail. 

Dr.  H.  Hahnemann  (LFA,  Institute  for  Motor  Research)  went  to 
considerable  trouble  in  order  to  perfect  an  elaboration  of  this  lecture 
series  which  I  examimed  and  supplemented  in  a  few  places.  Miss  Hildegard 
Munz  participated  in  the  illustration.  To  both  I  owe  my  most  sincere 
thanks  for  this  collaboration. 


Schlichting 


Aerodynamisch.es  Institut 

der  Technischen  Hochschule,  Braunschweig 

October  19^2. 
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INTRODUCTION 


Gentlemen:  In  the  lecture  series  starting  today  I  want  to  give 
you  a  survey  of  a  field  of  aerodynamics  which  has  for  a  number  of  years 
been  attracting  an  ever  growing  interest.  The  subject  is  the  theory  of 
flows  with  friction,  and,  within  that  field,  particularly  the  theory  of 
friction  layers,  or  boundary  layers. 

As  you  know,  a  great  many  considerations  of  aerodynamics  are  based 
on  the  so-called  Ideal  fluid,  that  is,  the  frictionless  incompressible 
fluid.  By  neglect  of  compressibility  and  friction  the  extensive  mathe¬ 
matical  theory  of  the  ideal  fluid  (potential  theory)  has  been  made 
possible. 

Actual  liquids  and  gases  satisfy  the  condition  of  incompressibility 
rather  well  if  the  velocities  are  not  extremely  high  or,  more  accurately, 
if  they  are  small  in  comparison  with  sonic  velocity.  For  air,  for 
instance,  the  change  in  volume  due  to  compressibility  amounts  to  about 
1  percent  for  a  velocity  of  60  meters  per  second. 

The  hypothesis  of  absence  of  friction  is  not  satisfied  by  any 
actual  fluid;  however,  it  is  true  that  most  technically  important  fluids, 
for  instance  air  and  water,  have  a  very  small  friction  coefficient  and 
therefore  behave  In  many  cases  almost  like  the  ideal  frictionless  fluid. 

Many  flow  phenomena.  In  particular  most  cases  of  lift,  can  be  treated 
satisfactorily,  -  that  Is,  the  calculations  are  in  good  agreement  with 
the  test  results,  -  under  the  assumption  of  frictionless  fluid.  However, 
the  calculations  with  frictionless  flow  show  a  very  serious  deficiency; 
namely,  the  fact,  known  as  d'Alembert's  paradox,  that  in  frictionless 
flow  each  body  has  zero  drag  whereas  in  actual  flow  each  body  experiences 
a  drag  of  greater  or  smaller  magnitude.  For  a  long  time  the  theory  has 
been  unable  to  bridge  this  gap  between  the  theory  of  frictionless  flow 
and  the  experimental  findings  about  actual  flow.  The  cause  of  this 
fundamental  discrepancy  is  the  viscosity  which  is  neglected  in  the  theory 

*"*Vortragsreihe  ’•Grenzschlchttheorie . '  Tell  A:  Laminar e  Str <5mungen . " 
Zentrale  fur  wissenschaftliches  Berichtswesen  der  Luftfahrtforshung  des 
Generalluftzeugmeisters  (ZWB)  Berlin-Adlershof,  pp.  1-153*  Given  in  the 
Winter  Semester  1941/42  at  the  Luftfahrtforschungsanstalt  Hermann  Goring, 
Braunschweig.  The  original  language  version  of  this  report  is  divided  into 
two  main  parts.  Tell  A  and  Tell  B,  which  have  been  translated  as  separate 
NACA  Technical  Memorandums,  Nos.  1217  and  1218,  designated  part  I  and  part  II, 
respectively. 
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of  the  ideal  fluid;  however,  in  spite  of  its  extraordinary  smallness  it 
is  decisive  for  the  course  of  the  flow  phenomenon.  As  a  matter  of  fact 
the  problem  of  drag  can  not  he  treated  at  all  without  taking  the  viscosity 
into  account. 

Although  this  fact  had  been  known  for  a  long  time,  no  proper  approach 
to  the  theoretical  treatment  of  the  drag  problem  could  be  found  until  the 
beginning  of  the  present  century.  The  main  reason  was  that  unsurmountable 
mathematical  difficulties  stood  in  the  way  of  theoretical  treatment  of  the 
flow  phenomena  of  the  viscous  fluid.  It  is  Professor  Prandtl's  great  merit 
to  have  shown  a  way  to  numerical  treatment  of  viscosity,  particularly  of 
the  technically  important  flows  under  consideration  and  thereby  to  have 
opened  up  new  vistas  on  many  important  perceptions  about  the  drag  problem 
and  related  questions.  Prandtl  was  able  to  show  that  in  the  case  of  most 
of  the  technically  important  flows  one  may  treat  the  flow,  as  a  whole,  as 
frictionless  and  utilize  the  simplifications  for  the  calculation  thus  made 
possible,  but  that  in  the  immediate  neighborhood  of  the  solid  walls  one 
always  had  to  take  the  friction  into  consideration.  Thus  Prandtl  subdivides, 
for  the  purpose  of  calculation,  the  flow  surrounding  a  body  into  two 
domains:  a  layer  subject  to  friction  in  the  neighborhood  of  the  body,  and  a 
frictionless  region  outside  of  this  layer.  The  theory  of  this  so-called 
"Prandtl* s  friction  or  boundary  layer"  has  proved  to  be  very  fruitful  in 
modem  flow  theory;  the  present  lecture  will  center  around  it. 

At  this  point  I  want  to  indicate  a  few  applications  of  the  boundary- 
layer  theory.  A  first  important  application  is  the  calculation  of  the 
frictional  surface  drag  of  bodies  immersed  in  a  flow,  for  instance,  the 
drag  of  a  flat  plate  in  longitudinal  flow,  the  frictional  drag  of  a  ship, 
a  wing  profile,  and  an  airplane  fuselage.  A  special  property  of  the 
boundary  layer  is  the  fact  that  under  certain  circumstances  reverse  flow 
occurs  in  the  immediate  proximity  of  the  surface.  Then,  in  connection 
with  this  reverse  flow,  a  separation  of  the  boundary  layer  takes  place, 
together  with  a  more  or  less  strong  formation  of  vortices  in  the  flow 
behind  the  body.  Thus  a  considerable  change  in  pressure  distribution, 
compared  with  frictionless  flow,  results,  which  gives  rise  to  the  form 
drag  of  the  body  immersed  in  the  flow;  The  boundary— layer  theory  therefore 
offers  an  approach  to  the  calculation  of  this  form  drag.  Separation  occurs 
not  only  in  the  flow  around  a  body  but  also  in  the  flow  through  a  divergent 
tunnel. 

Thus  flow  phenomena  in  a  diffuser,  as,  for  instance,  in  the  bucket 
grid  of  a  turbine,  may  be  included  in  boundary— layer  theory.  Furthermore, 
the  phenomena  connected  with  the  maximum  lift  of  a  wing,  where  flow 
separation  is  concerned,  can  be  understood  only  with  the  aid  of  boundary- 
layer  theory.  The  problems  of  heat  transfer  also  can  be  explained  only 
by  boundary-layer  theory. 

Aa  will  be  shown  in  detail  later,  one  must  distinguish  between  the 
two  states  of  boundary— layer  flow  —  laminar  and  turbulent;  their  flow  laws 
are  very  different.  Accordingly,  the  lecture  will  be  divided  into  three 
main  parts:  1.  Laminar  flows,  2.  Turbulent  flows,  3.  Laminar— turbulent 
transition.  Although  the  boundary  layer  will  be  our  main  consideration, 
it  will  still  be  necessary  as  preparation  to  discuss  to  some  extent  the 
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general  theory  of  the  viscous  fluid.  This  will  be  done  in  the  first 
chapter. 


CHAPTER  I.  VISCOSITY 


Every  fluid  offers  a  resistance  to  a  form  variation  taking  place  in 
finite  time  interval,  which  is  of  different  magnitude  according  to  the 
type  of  fluid.  It  is,  for  instance,  very  large  for  syrup  or  oil,  but 
only  small  for  the  technically  important  fluids  (water,  air). 

The  concept  of  viscosity  can  be  best  made  clear  by  means  of  a  test 
according  to  figure  1: 

Let  fluid  be  between  two  parallel  plates  lying  at  a  distance  h 
from  each  other.  Let  the  lower  plate  be  fixed,  while  the  upper  plate  is 
moved  with  the  velocity  Uq  uniformly  and  parallel  to  the  lower  one. 

For  moving  the  upper  plate  a  tangential  force  P  must  be  expended  which 
Is 

P  =  (1.1) 

h 


according  to  experiment,  where  F  is  the  area  of  the  upper  plate  and  u 
is  a  constant  of  proportionality.  (End  effects  are  not  included).  The 
quantity  u  is  called  the  viscosity  coefficient  or  the  dynamic  viscosity. 

Since  the  phenomenon  in  question  is  a  parallel  gliding,  the  transverse 
velocity  component  in  the  y-direction,  denoted  by  v,  equals  zero.  The 
fluid  adheres  to  the  upper  and  lower  surface,  respectively,  a  linear 
velocity  distribution  between  the  plates  is  set  up,  the  magnitude  of 
which  depends  solely  on  y. 


u(y) 


or: 


du  _  u0 
dy  ~  h 


Since  for  y  =  0:  u  =  0,  for  y  =  h:  u  =  Uq.  If  one  designates  P/F, 

the  tangential  force  per  unit  area,  as  the  frictional  shearing  stress 
there  follows: 
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The  dimensions  of  n  are  accordingly  kg  sec/m2.  A  flow  as  represented 
in  figure  1,  where  no  transverse  velocity  occurs  and  the  shearing  stress 
at  all  points  of  the  flow  is  therefore  given  by  equation  (1.2),  is  called 
simple  shear  flow.  In  the  special  cas9  described,  the  shearing  stress  is 
everywhere  of  equal  magnitude,  and  equal  to  that  at  the  surface.  Besides 
the  dynamic  viscosity  p.  the  concept, of  kinematic  viscosity  V  is 

required,  which  for  the  density  pH*  sec  2nr~^j  is  defined  as 


V  =  t  [m2/s] 


For  20°  C,  v  is,  for  instance,  for  water: 


V  =  1.01  X  ICT6  m2/s 


for  air: 


v  =  14.9  x  ICT6  m2/s  i  lcH*'  m2/s 


if  the  air  pressure  has  the  standard  value  p0  =  760  mm  hg. 

CHAPTER  II.  POISEUTLLE  FLOW  THROUGH  A  PIPE 


The  elementary  empirical  friction  law  of  the  simple  shear  flow 
derived  above  permits  the  immediate  determination  of  the  flow  and  the 
resistance  in  a  smooth  pipe  of  circular  cross  section  and  of  constant 
diameter,  d  =  2r.  At  a  very  large  distance  from  the  beginning  of  the 
pipe  one  cuts  off  a  piece  of  pipe  of  length  l  (fig.  2)  and  examines 
the  cylinder  of  diameter  2y,  the  axis  of  which  is  identical  with  the 
pipe  axis.  According  to  what  has  been  said  so  far,  the  velocity  probably 
will  be  again  a  function  of  y.  A  pressure  difference  p-j^  -  p2  is 

required  for  forcing  the  fluid  through  the  cylinder.  According  to 
practical  experience,  the  static  pressure  across  every  cross  section  may 
be  regarded  as  constant.  The  flow  is  assumed  to  be  steady  and  not 
dependent  on  the  distance  from  the  beginning  of  the  pipe.  Equilibrium 
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must  then  exist  between  the  pressure  and  the  frictional  shearing  stress 
which  attempts  to  retard  the  motion.  Thus  for  the  cylinder  of  radius  y 
the  following  equation  is  valid:  pressure  force  difference  acting  at 
the  cross  sections  =  frictional  force  acting  along  the  cylinder  wall,  or 


jry2  =  2itylT 


(2.1) 


or 


Pi  ~  p2  y 
l  2 


(2.1a) 


Since  flow  parallel  to  the  axis  is  to  be  expected,  one  takes  from  the 

previous  paragraph,  t  =  —  n—  (the  minus  sign  indicates  that  the 

dy 

velocity  diminishes  with  increasing  distance  from  the  axis;  thus  du/dy 
is  negative,  the  shearing  stresses  under  consideration,  however,  are 
positive),  and,  after  separation  of  the  variables,  du  becomes: 


du  =  — 


Pi  ~  *>2 

U2 


|  dy 


(2.2) 


and,  on  integration: 


u(y) 


Pi  ~  ?2 


From  the  fact  that  for  y  =  r  the  velocity  is  supposed  to  be  u(y)  =  0 
follows  that  the  constant  of  integration  C  has  to  be  C  =  r2/4.  Thus: 


^  (ra  - y2) 


(2.3) 
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This  equation  (2.3)  is  Poissuille’s  law  for  pipe  flow.  It  states  that 
the  velocity  u(y)  is  distributed  parabolically  over  the  pipe  cross 
section.  The  apex  of  the  parabola  lies  on  the  pipe  axis;  here  the 
velocity  is  greatest,  namely: 


umax 


P1  ~  P2  r2 


Therewith  one  may  write  (2.3): 


=  w  (x  -  4) 


(2.4) 


(2.3a) 


By  Poiseuille*s  law  (equation  (2.3))  the  drag  of  the  developed  laminar 
flow  (which  is  proportional  to  p-^  —  Pg)  is  directly  proportional  to 

the  first  power  of  the  velocity. 


This  statement  is  characteristic  of  all  kinds  of  laminar  flow 
whereas,  as  will  be  seen  later,  the  drag  in  turbulent  flow  is  almost 
proportional  to  the  second  power  of  the  velocity. 


The  flow  volume  for  the  present  case  remains  to  be  given.  With 


designating  an  area  element,  Q  is  Q  = 
velocity  paraboloid,  therefore 


u(y)  d  F  =  volume  of  the 
1 


dF 


u 

max 


(2.5) 


This  flow  law  is  often  used  for  determination  of  the  viscosity,  by 
measuring  the  quantity  flowing  subjected  to  a  pressure  gradient  (usually 
produced  by  gravity  in  a  vertical  capillary  tube).  Of  course,  the 
starting  losses- must  be  taken  into  consideration  which  due  to  the  mixing 
zone  (vortex  formation)  at  the  pipe  end  are  not  recovered  to  their  full 
extent. 

A  drag  coefficient  \  will  now  be  defined.  Since  turbulent  flows 
are  more  important  than  laminar  ones  and  since  the  drag  in  turbulent 
flow  increases  about  as  the  square  of  the  velocity,  \  will  also  be 
referred  zo  u^. 

For  flow  problems,  let  \  thus  be  defined  as:  ratio  of  the  pressure 
drop  along  a  test  section  of  a  specified  characteristic  length  to  the  dynamic 
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pressure  q  =  pu2/2,  with 
across  the  cross  section). 


u  =  =  the  mean  velocity 

itr2 
Then: 


(average  taken 


\ 


_  djD  A 


(2.6) 


with  A  =  characteristic  length,  thus,  for  the  present  case,  the  pipe 
Aiameter,  anA  with  \  =  Aimensionless  quantity.  For  the  present  AevelopeA 
laminar  pipe  flow,  accorAing  to  equation  (2.5) 


Thus: 


dp  =  £l 

Ax 


^2  _  8nQ 


l 


2A  8pQ 

_2  4 

pu  irr 


2d8p  Q 

...  2  2  _ 
pur  jtr  u 

=  1 


6W  P 
uA 


or 


\ 


64 

Be 


(2.7) 


with  the  dimensionless  quantity  Be  =  —  signifying  the  Beynolds  number 

of  the  circular  pipe.  Since  the  pressure  drop  which  is  only  linearily 
dependent  on  the  velocity  was  referred  to  u2,  then,  for  laminar  flow: 

\  ~  A  logarithmic  plot  of  A.  =  f(Ee)  or  \  =  f(u)  therefore  results 
u 

in  a  straight  line  inclined  45°  toward  the  Be— axis  (compare  fig.  82  Part  II.''. 

After  this  short  analysis  of  the  one-dimensional  case  of  viscous 
fluid  we  will  now  consider  the  three-dimensional  case. 
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CHAPTER  III.  EQUATIONS  OF  MOTION  OF  THE  VISCOUS  FLUID 

a.  State  of  Stress 


For  this  purpose  one  must  know  first  of  all  the  general  state  of 
stress  in  a  moving  viscous  fluid  and  must  then  connect  this  state  of 
stress  with  the  state  of  deformation.  For  the  deformation  of  solid 
todies  the  resistance  to  the  deformation  is  put  proportional  to  the 
magnitude  of  the  deformation  (assuming  the  validity  of  Hooke’s  law). 

For  flowing  fluids,  on  the  other  hand,  the  resistance  to  deformation 
will  depend  on  the  deformation  velocity,  that  is,  on  the  variation  of 
velocity  in  the  neighborhood  of  the  point  under  consideration.  (Solid 
bodies:  displacement  gradient  =  displacement  per  second.  Fluid: 
velocity  gradient). 

One  starts  from  the  basic  law  of  mechanics  according  to  which: 
mass  x  acceleration  =  sum  of  the  acting,  or  reBultant  force.  For  the 
mass— per— unit  volume,  that  is,  the  density  p,  one  may  write  the  law 


P 


*Dw 

—  =  K  +  R 
Dt  “  “ 


+  F 


(3.1) 


Dw 

—  =  substantial  acceleration 
Dt 

K  =  mass  forces 

R  =  surface  forces,  composed  of  pressure  forces  normal  to  the 
surface  and  frictional  forces  in  the  direction  of  the 
surface 

F  =  negligible  extraneous  forces 


In  order  to  formulate  the  surface  forces,  one  Imagines  a  small  rectangular 
element  of  volume  dV  =  dx  dy  dz  cut  out  of  the  flow  (fig.  3)  the  left 
front  corner  of  which  lies  at  the  point  (x,  y,  z).  The  element  is  to  be 
very  small  so  that  only  the  linear  variations  of  a  Taylor  development 
need  to  be  taken  into  consideration;  on  its  surfaces  dy  dz  act  the 
resultant  stresses  (vectors); 


a*x 


p  or  p  +  — > -  dx,  respectively 
ox 


(3-2) 


*  Throughout  the  text,  underscored  letters  are  used  in  place  of 
corresponding  German  script  letters  used  in  the  original  text. 
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(The  index  x  signifies  that  the  stress  tensor  acts  on  a  surface  element 
norma]  to  the  x— direction) . 


Analogous  terms  result  for  the  surfaces  dz  dx  normal  to  the 
y— axis  and  dx  dy  normal  to  the  z-axis,  if  x  in  equation  (3*2)  is 
replaced  everywhere  by  y  or  z,  respectively.  From  this  there  results 
as  components  of  the  resultant  force: 

^Px 

Force  on  the  surface  element  normal  to  the  x— direction:  -  dx  dy  dz 

dx 


Force  on  the  surface  element  normal  to  the  y-direction: 


% 

- —  dy  dz  dx 
oy 


dp 

Force  on  the  surface  element  normal  to  the  z— direction:  — -  dz  dx  dy 

dz 


The  total  resultant  surface  force  B  per  unit  volume  caused  by  the 
state  of  stress  is  therefore: 


B 


dx  dy  dz 


(3-3) 


p  ,  p  end  p  are  vectors  which  can  be  further  decomposed  into 
-x  -y  -z 

components.  In  this  decomposition  the  components  normal  to  every  surface 
element,  that  is,  the  normal  stresses,  are  designated  by  a  (indicating 
by  the  index  the  direction  of  this  normal  stress);  the  other  components 
(tangential  stresses)  ace  denoted  by  t  (with  double  index:  the  first 
indicates  to  which  axis  the  surface  element  is  perpendicular,  the  second, 
the  axial  direction  of  the  stress  t).  With  these  symbols  there  is: 


ia  +  Jt  +  kr 
x  xy  xz 

iT  +  ja  +  kT 
yx  y  yz 

It  +  Jt  +  ka 


i 


zx 


z 

) 


(3.4) 


This  state  of  stress  represents  a  tensor  with  nine  vector  components, 
which  can  be  characterized  by  the  stress  matrix  (stress  tensor): 
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ffx 

V 

Txz  \ 

Tyx 

T'yz  j 

(3.5) 

''yZX 

Tzy 

that 

those 

of  the  six  tangential 

stresses 

which  have  the  same  Indices,  although  in  interchanged  sequence,  must  be 
equal.  This  follows  for  a  homogeneous  state  of  stress  from  the  equilibrium 
of  the  small  cube  di  djr  dz  with  respect  to  rotation: 


Since  t  djr  di  is  the  force  attempting  to  rotate  the  cube  counter— 
xy 

clockwise  about  the  z-axis,  (seen  from  above  in  fig.  3).>  with  the  lever 
arm  dx,  and  since,  correspondingly,  the  force  -t  dx  dz  attempts  to 

Jf  Jw 

rotate  the  cube  clockwise  about  the  z~axis,  with  the  lever  arm  dy  the 
balance  of  moments  requires: 


T  dy  dz  dx  -  Txy  it  dz  dy  *  0,  thus  Txy  *  t^. 


Correspondingly,  because  of  freedom  from  rotation  about  the  x-axis 

t  =  t  ,  and  because  of  freedom  from  rotation  about  the  y-axis 
yz  zy* 

t  =  t  .  the  nine  components  of  the  stress  tensor  are  reduced  to  six 
zx  xz* 

and  the  stress  matrix  (equation  (3»5))  is  converted  into  the  stress 
matrix  symmetrical  with  respect  to  the  principal  diagonal: 


(stress  matrix)  (3»6) 


For  the  frictional  force  one  obtains  according  to  equation  (3*3)  by 
insertion  of  the  components  from  equation  (3.4)  and  by  reduction  to  the 
six  remaining  terms  according  to  equation  (3.6): 
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do  St  St 

R  =  H  - +  — -  +  - 1  +  J 

V  Sx  by  Sz , 


St  Sa  St 
*y  +  _y_  +  yz 

Sx  Sy  Sz 


=  x— component 


+  k 


+  y-component 


xz 


St, 


yz 


Sx 


by 


ba2 

Sz, 


(3-7) 


+  z-component. 

For  the  case  of  the  frictionless  (ideal)  fluid,  all  shearing 
stresses  disappear 


t  =  t  =  t  =  0 
xy  yz  zx 


(3.8) 


and  only  the  normal  stresses  remain,  which  in  this  case  are  all  equal. 
Since  the  normal  stresses  from  within  toward  the  outside  are  denoted  as 
positive,  the  normal  stresses  equal  the  negative  fluid  pressure: 

ax  =  <iy  =  az  =  -p  (3-9) 

The  static  pressure  equals  the  negative  arithmetic  mean  of  the  normal 
stresses: 

-p  =  —  ( a  +  a  +  o  \ 

3  V  x  y  z) 


b.  State  of  Deformation 

The  state  of  stress  treated  so  far  is,  alone,  not  very  useful. 
Therefore  we  will  now  consider  the  state  of  deformation  (that  is  the 
field  of  velocity  variations)  and  then  set  up  the  relations  between 
state  of  stress  and  state  of  deformation. 
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Let  the  velocity  w^  with  the  components  u^,  v^,  in  the 

directions  of  the  axes  exist  at  the  point  A  the  coordinates  of  which 
are  XA*  yA»  za* 


If  one  limits  oneself  to  the  points  x,  y,  z  in  th9  immediate 
neighborhood  of  A  with  the  velocity  w  =  iu  +  jv  +  kw,  and  if  one 
limits  oneself  —  as  also  in  setting  up  the  state  of  deformation  —  to 


linear  terms  only,  one  obtains  for  the  deformation  the  relative  change 


in  position  between  the  points  x,  y,  z  and  x  ,  y  , 

A  A 


z 

A 


per  unit 


time,  that  is,  the  difference  of  the  velocities  at  the  points  x,  y,  z 

V  yA,  V 


dw  as  i  du  +  J  dv  +  k  dw  =  distortion  of  the  fluid  region  in  the 
neighborhood  of  the  point  A. 

Omitting  the  index  A  one  obtains  therefore: 


dw  =  —  dx  + 
dx 


dtf  ,  dw,  ./du,  du,  du, 

-=■  dy  +  —  dz  =  i  (  —  dx  +  —  dy  +  —  dz 

dy  dz  \dx  dy  dz 


(3.11a) 


dx  +  ^  dy  +  ^  dz 
dy  dz 


Thus  the  velocity  variation  (and  hence,  on  integration,  the  velocity 
itself)  in  the  neighborhood  of  the  point  A  is  known  if  th9  nine 
partial  derivatives  of  the  velocity  components  with  respect  to  the  space 
coordinates  are  known.  Corresponding  to  the  stress  matrix,  one  may  form 
a  deformation  matrix: 
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du 

du 

du 

dx 

dy 

dz 

dv 

d v 

dv 

dx 

dy 

dz 

dw 

dw 

dx 

dy 

dz 

deformation  matrix 


(3.12) 


The  friction  forces  of  the  viscous  fluid  are  given  by  a  relation  (which 
will  have  to  be  determined)  between  these  two  matrices.  First,  the 
deformation  matrix  is  to  be  somewhat  clarified. 

1.  Case  of  pure  elongation. 


One  assumes  u  —  u_  =  a(x  —  xA),  with  a  =  ^1  =  constant.  Let  all 

a  a  dx: 

other  terms  of  the  matrix  disappear;  the  matrix  will  then  appear  as 
follows : 


Then  the  velocity  variation  is  simply  du  *  a  dx,  and  u  =  ax. 

All  points  of  the  y— axis  remain  at  rest,  the  points  to  the  right  and 
left  of  it  are  elongated  or  compressed,  according  to  whether  a  >  0 
or  a  <0  (fig.  5).  The  equation  u  =  ax  therefore  represents  an 
elongation  or  expansion  parallel  to  the  x— axis.  Corresponding  relations 
apply  for  the  other  terms  of  the  principal  diagonal  of  the  matrix. 

2.  Case  of  pure  translation 


All  terms  disappear;  the  matrix  then  reads: 


0  0  0 
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In  this  case,  which,  as  a  matter  of  fact,  should  have  teen  mentioned 
first,  u  —  u  =  0,  du  =  0;  u  =  constant.  The  velocity  component 

a, 

parallel  to  the  x-axis  is  uniform  (correspondingly  for  the  other  axes). 
3.  Case  of  angular  deformation. 


One  assumes  u  —  u  =  e(y  —  y  ) ,  with  e  =  —  =  constant.  All 

a  w  a" 

other  terms  equal  zero,  and  the  matrix  reads: 


0  e  0 
0  0  0 
0  0  0 


du  =  e  dy  and  u  =  ey;  that  is,  all  points  of  the  x-axis  retain  their 
position;  all  points  of  the  y-axis  shift  to  the  right  (left),  when  e  >  0 
(e  <  0);  for  e  >  0  the  y-axis  is  rotated  clockwise  by  the  angle  e 
(because  of  the  linearity).  The  y-axis  is  simultaneously  elongated.  The 
phenomenon  in  question  is  therefore  a  shearing  (fig.  6),  with  tan  €  =  e. 

Correspondingly  there  results  for  v  —  v  =  f  (x  —  x.)  and 

fit  A 


0 

0 

0 


dv  =  f  dx;  v  =  ix 


All  points  of  the  y— axis  retain  their  position;  the  points  of  the  x-axis 
are  rotated  by  the  angle  5;  tan  5  =  f  (fig. 7)*  Terms  outside  of  the 
principal  diagonal  of  the  matrix  result  therefore  in  a  deformation  of 
the  right  angle  with  axis-elongation  (shearing).  The  right  angle  between 
the  x—  and  y— axes  is,  therefore,  for  e  >  0  and  f  >  0  deformed  by 


_  ov  +  ou  _  ^  _  Defo:rmation  about  the  z— axis 

dx  dy 


€  +  6 
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Correspondingly : 


du  dw 
221  dz  dx 


deformation  about  the  y-axis 


7  =  Hi  +  HL  =  deformation  about  the  x— axis 

yz  dy  dz 


(The  deformation  angles  are  herein  regarded  as  small  so  that  the  tangent 
may  be  replaced  by  the  argument). 


c.  Navier— Stokes  Formulation  for  the  Stress  Tensor 


One  now  proceeds  to  relate  the  stress  matrix  (equation  (3.6))  with 
the  deformation  matrix  (equation  (3.12)).  The  former  is  symmetrical  with 
respect  to  the  principal  diagonal,  bub  not  the  latter.  However,  one 
obtains  a  symmetrical  deformation  matrix  by  adding  to  equation  (3.12) 
its  reflection  in  the  principal  diagonal.  Furthermore,  one  first  splits 
off  the  pressure  p  (contribution  of  the  ideal  fluid)  from  the  stress 
matrix  and  sets  the  remaining  stress  matrix,  according  to  Stokes, 
proportional  to  the  deformation  matrix  made  symmetrical: 


x  xy  xz 


t  „  <t  r , 


!  T  T  a 

i  xz  yz  uz 


i  p 

0 

0 

du 

du 

du 

£u 

dz 

dw  , 

I  * 

I 

! 

dx 

dy 

dz 

dx 

dx 

dx 

i 

i 

i  =  -  0 

p 

0 

■ 

+  u 

dl 

dl 

dz 

+  Li 

& 

dz 

dw 

Sr 

dx 

dy 

dz 

i 

dy 

dy 

dy  ' 

i 

j  0 

0 

p 

dw 

dw 

dw 

du 

dz 

dw 

Sr 

dx 

dy 

dz 

dz 

dz 

dz 

(3.13) 


From  equation  (3. 13)  each  stress  component  may  be  given  immediately  by 
coordinating  the  homologous  parts  of  the  matrices  to  each  other.  For 
instance : 


—  p  +  2  |i  —  =  static  pressure  +  pressure  due  to 
dx  velocity  variation,  or: 
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T 


Thus  for  one-dimensional  flow  t__  =  t__  = 

■*v 


Furthermore,  there  follows  tram,  equation  (3.I3): 


or 


because 


du 

dx 


+ 


=  div  v  =  0 
oz 


(3.14) 


for  the  incompressible  flows  free  of  sources  and  sinks  under  considera¬ 
tion.*  Thus  for  the  viscous  incompressible  flow,  as  for  the  ideal  fluid 
the  pressure  equals  the  arithmetic  mean  of  the  normal  stresses. 

With  these  results  the  components  of  the  friction  force  may  be 
expressed  according  to  equation  (3*7)  as  follows: 


*  The  compressibility  manifests  itself  as  normal  stress,  since  it  can 
be  interpreted  as  a  pressure  disturbance,  for  instance  due  to  variation 
in  density,  which  attempts  to  spread  in  all  directions  -  considered 
infinitesimally . 
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*x  = 


da  dr 
x  xy  .  xz 


dx 


+  _^  =  _3i+  J&  +  &+  sfa 

*  3z  Sx  W  3y2  Sz2 


,  ( a2u  d‘:v  d£w_\ 

y  0^2  dxdj  / 


or,  since  div  w  =  0, 


x  =  "  s + 


Ey  =  -|+^T 


E_  =  — 


4-  iiAv 

Sz 


(3.15) 


in  which  An  =  &  +  &  +  & 
dx2  dy2  dz2 

If  one  finally  designates  the  mass  forces  by  K  =  p(iX  +  JY  +  kZ),  and 
assumes  the  decomposition  of  the  substantial  derivative  into  a  local 
and  convective  part  as  known  from  Euler* s  equation,  one  obtains  for  the 
components  of  the  equation  of  motion  of  the  non— stationary,  incompressible, 
and  viscous  fluid  from  equation  (3.1): 


T  &  + 

w&' 

)=  PX  - 

dp 

mAu 

dx 

ay 

az  j 

ax 

av 

V  —  + 

dv  ^ 
w  - 

)=  PY  - 

-  &  + 

pAv 

ax 

ay 

az  / 

dy 

aw 

ax 

V  &  + 

oy 

=  pZ  - 

_  + 

az 

mAw 

> 


(3.16) 


In  addition,  the  continuity  equation 
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is  used.  Written  in  vector  form,  the  Navier-Stokes  differential  equation 
and  the  equation  of  continuity  read 


—  =  K  —  i  grad  p  +  VAw 
Dt  “  P 


(3.18) 


div  w  =  0 


(3.19) 


Due  to  the  friction  terms,  therefore,  terms  of  the  second  order  enter 
the  differential  equation. 

Boundary  conditions  are  attached  to  these  equations.  If  all  friction 
terms  on  the  right  side  are  cancelled,  that  is  V  =  0,  the  differential 
equations  Become  equations  of  the  first  order  and  one  Boundary  condition 
is  sufficient,  namely  the  boundary  condition  of  the  potential  flow: 
vn  =  0  on  the  Bounding  walls. 

This  means  that  the  normal  component  vn  of  the  velocity  at  the 

Bounding  surface  must  disappear  on  the  surface  itself  whereas  the  fluid 
still  can  glide  parallel  to  the  Boundary  (tangential  velocity 

parallel  to  the  surface  /  0). 

For  viscous  flow  where  the  differential  equation  is  of  the  second 
order,  two  Boundary  conditions  are  required,  namely: 


vn  =  0  and  vt  =  0  (condition  of  no  slip)  (3-20) 

that  is,  the  fluid  must  in  addition  adhere  to  the  surface. 

Second  lecture  (Dec.  8,  1941) 


CHAPTER  17.  GENERAL  PROPERTIES  OF  NAVIER-STOKES  EQUATIONS 


These  Navier-Stokes  differential  equations  represent  together  with 
the  equation  of  continuity  a  system  of  four  equations  for  the  four 
unknown  quantities  u,  v,  w,  p.  On  the  left  side  of  the  Navier-Stokes 
differential  equations  are  the  inertia  terms,  on  the  right  side  the  mass 
forces,  the  pressure  forces,  and  the  friction  forces. 

Since  Stokes’  formulation  is,  of  course,  at  first  purely  arbitrary. 
It  is  not  a  priori  certain  whether  the  Navier-Stokes  differential 
equations  describe  the  motion  of  a  fluid  correctly.  They  therefore 
require  verification,  which  is  possible  only  By  way  of  experimentation. 
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Unfortunately,  due  to  un surmountable  mathematical  difficulties,  a 
general  solution  of  the  differential  equation  is  not  yet  known,  that  is, 
a  solution  where  inertia  and  fric'ion  terms  in  the  entire  flow  region 
are  of  the  same  order  of  magnitude.  However,  known  special  solutions 
(for  instance,  the  pipe  flow  with  predominant  viscosity  or  cases  with 
large  inertia  effect)  agree  so  well  with  the  experimental  findings, 
that  the  general  validity  of  Navier— Stokes  differential  equation  hardly 
seems  questionable. 

The  plane  problem: 

By  far  the  greatest  part  of  the  application  of  Navier-Stokes 
differential  equations  concern  "plane"  cases,  that  is,  the  cases  where 
no  fluid  flows  in  one  direction.  The  velocity  vector  w  is  then  given 


w  =  iu  (x,  y,  t)  +  Jv  (x,  y,  t) 


(4.1) 


since  w  =  0.  The  equation  system  (equations  (3.16)  and  (3*17))  then  is 
transformed  into  the  3  equations 


/d  u  du  du 

Pl^+U^+T^ 

V  dt  3x  by 


=  pX  - 


1p+  +  asA 

dx1 2  by2/ 


bx 


Vdt  Sr  S 7/  Sy  ^2 


(4.2) 


$2  +  =  0 

bx  by 


with  the  three  unknown  factors  u,  v,  p  (X  and  Y  are  the  components 
of  the  mass  force  K  per  unit  volume).- 

After  various  minor  transformations  the  equation  system  may  be 
written  as  a  single  equation.  To  this  end  one  introduces  the  rotational 
vector  rot  j£  which  for  the  plane  case  has  only  one  component  not 
equalling  zero: 


1  rot.v  .  <■>.  =  i  -  &) 

2  z  L  2  \dx  by/ 


(4.3) 
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Furthermore,  the  mass  force  in  equation  (4.2)  is  put  equal  to  zero. 
This  is  permissible  in  all  cases  where  the  fluid  is  homogeneous  and  no 
free  surfaces  are  present.  In  order  to  introduce  <x>z  into  equation  (4.2), 
the  first  equation  of  (4.2)  is  differentiated  with  respect  to  y,  and 
the  second  with  respect  to  x;  then  the  first  is  subtracted  from  the 
second  and  one  obtains: 


or 


(t.t) 


\  dt  dx  dy  / 


(*.5) 


With  this  transformation  the  pressure  terms  have  been  eliminated. 
Equation  (4.5)  may  now,  with  (i/p  =  V,  be  written: 


=  v&u 
Dt 


(vorticity  transport  equation) 


(4.6) 


with  a)  =  o>z  being  denoted  as  the  vorticity. 

This  equation  signifies:  The  convective  (substantial)  variation 
of  the  vortex  strength  equals  the  dissipation  of  vorticity  by 
friction. 

Equation  (4.6)  forms  with  the  equation  of  continuity  a  Bystem  of 
two  equations  with  two  unknowns,  namely  u  and  v,  the  derivatives  of 
which  define  co. 

By  introducing  a  flow  function  \|r(x,y)  one  may  finally  introduce 
a  single  equation  with  the  unknown  t|c.  The  flow  function  represents 
the  Integral  of  the  equation  of  continuity.  One  sets: 


That  is,  therefore,  the 
equation  of  continuity 
is  identically 
satisfied  by 


u  = 


v  =  — 


dx 
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Moreover, 


to  =  i  _  2u\  =  -  +  a£i\  =  _  it# 

2U  *)  2^2  tf)  2 


That  is: 
as  large 
equation 


The  Laplacian  of  the  flow  function  is  exactly  minus  two  times 
as  the  vorticity  (angular  velocity).  With  this  result 
(4.5)  "becomes,  after  division  by  P: 


+  u 


+  v  jAijf  =  VAA\|r 
ox  dy  / 


(k.9) 


or  expressed  only  in  i|r  with  equation  (4.7): 


at  ay  Sx  ax  ay  v 


(4.10) 


This  one  equation  with  the  unknown  t  is  the  vorticity  transport 
equation,  but  written  in  terms  of  t. 

The  inertia  terms  are  again  on  the  left,  the  viscosity  terms  on  the 
right  side.  Equation  (4.10)  is  a  differential  equation  of  the  fourth 
order  for  the  flow  function.  Again,  its  general  solution  is  extremely 
difficult  because  of  the  non-linearity.  For  very  slow  (creeping)  motions 
the  friction  terms  very  strongly  predominate.  Then  one  may  set: 


AAt  =  0 


(4.11) 


This  simplification  is  permissible  only  because  the  differential  equation 
remains  of  the  fourth  order,  so  that  no  boundary  condition  is  lost. 
However,  being  linear,  thiB  equation  is  at  least  solvable.  It  appears 
also  in  the  theory  of  elasticity  where  it  is  designated  as  the  bipotential 
equation.  There  exists  a  solution  of  equation  (4.11)  by  Stokes  for 
moving  droplets.,  which  was  extended  by  Cunningham  to  very  small  drop 
diameters  (comparable  to  the  mean  free  path  of  the  molecules). 

Herewith  we  shall  conclude  the  more  general  considerations  and  turn 
to  the  boundary  layer  problem  proper,  limiting  ourselves  to  fluids  of 
very  small  viscosity  V. 

A  few  preparatory  considerations  will  lead  up  to  the  boundary  layer 
problem.  One  might  conceive  the  notion  of  simply  eliminating  all  the 
friction  terms  of  the  Navi er— Stokes’  differential  equation  in  the  case  of 
small  viscosity.  However,  this  would  be  fundamentally  wrong  as  will  be 
proved  below. 
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An  equation  which  is  completely  analogous  to  equation  (4.5)  occurs 
in  the  theory  of  heat  transfer: 


\dt 


\ 

—  +  v  — 

bx  by  J 


=  \ 


(4.12) 


where  the  velocity  components  are  retained  whereas  the  rotation  o 
replaces  the  temperature  $ ,  the  density  p  the  specific  heat  cp 

per  unit  volume,  and  the  viscosity  p.  the  thermal  conductivity  X..  On 
the  left  of  equation  (4.12)  stands  the  temperature  change  due  to 
convection,  on  the  right  the  change  due  to  heat  transfer. 

The  temperature  distribution  around  a  heated  body  immersed  in  a 
flow  with  the  free  stream  velocity  uQ  (for  instance  fig.  8)  is 

determined  by  the  differential  equation  (4.12).  One  perceives  intuitively 
that  for  small  uQ  the  temperature  increase  starting  from  the  body 
extends  toward  the  front  and  all  sides  far  into  the  flow  (solid  contour) 
whereas  for  large  uQ  this  influence  is  mainly  limited  to  a  thin  layer 
and  a  narrow  wake  (dashed  contour). 

The  analogy  of  equations  (4.12)  and  (4.5)  indicates  that  the 
friction— rotation  distribution  in  question  muBt  be  similar:  For  small 
free  stream  velocity  the  rotation  is  noticeable  at  large  distance  from 
the  body,  whereas  for  large  uQ  the  rotation  is  limited  to  the  immediate 
neighborhood  of  the  body. 

Thus  for  rapid  motions,  that  is,  large  Reynolds  numbers  (compare 
next  section),  one  expects  the  following  solution  of  Navier-Stokes' 
differential  equations: 

1.  In  the  region  outside  of  a  thin  boundary  layer  <*>  =  0,  that  is, 
potential  flow 

2.  Inside  this  thin  boundary  layer  cu  ^  0,  thus  no  potential  flow. 

Therefore,  one  must  not  set  o>  =  0  in  this  boundary  layer,  even 
for  small  viscosity. 

It  is  true  that  the  potential  flow  is  also  a  solution  of  Navier- 
Stokes'  differential  equations,  but  it  does  not  satisfy  the  boundary 
layer  condition  v^  =  0. 

Proof:  The  potential  flow  may  be  derived  from  potential  <J>(x,y,z) 
as: 

w  =  grad  $,  with  A<1>  =  2 —  +  2— _  +  — -  =  0 

dx2  dy2  bz2 


(4.13) 
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However ,  if  AD  =  0,  then  also  grad  AD  =  A  gradD  =  0,  that  is.  Aw  =  0 
for  potential  flow.  According  to  equation  (3.18)  this  fact  signifies 
that  in  the  Navier— Stokes  differential  equations  the  friction  terms 
vanish  identically,  and  hence  that  the  potential  flow  actually  satisfies 
the  Navier— Stokes  differential  equations.  However,  it  satisfies  only 
the  one  boundary  condition  vr  =  0. 

Thus,  for  the  limiting  case  of  small  viscosity,  one  obtains  useful 

solutions  for  the  limiting  process  v  - ^  0  not  by  cancelling  the 

friction  terms  in  the  differential  equation,  since  this  reduces  its 
order  (the  differential  equation  of  the  fourth  order  for  the  flow 
function  would  turn  into  an  equation  of  the  second  order;  the  Navier— 
Stokes  differential  equations  would  change  from  the  second  to  the  first 
order),  so  that  one  can  satisfy  only  correspondingly  fewer  boundary 
conditions. 

Thus  the  limiting  process  V - 0  must  not  be  performed  in  the 

differential  equation  itself,  but  only  in  its  solution. 

This  can  be  clearly  demonstrated  on  an  example  (referred  to  for 
comparison  by  Prandtl)  of  the  solution  of  an  ordinary  differential 
equation.  Consider  the  damped  oscillation  of  a  mass  point.  The 
differential  equation 


2 

m  +  k  +  cx  =  0  (4.1k) 

dt2 


applies  in  which  m  represents  the  oscillating  mass,  k  the  damping 
constant  and  c  the  spring  constant,  (x  =  elongation,  t  =  time). 

Let  for  instance  the  two  initial  conditions  be: 


t  =  0;  x  =  0;  dx/dt  =  1 


In  analogy  to  the  case  in  question  one  considers  here  the  limiting 
case  of  a  very  small  mass  m,  since  then  the  term  of  the  highest  order 
tends  toward  zero.  If  one  would  simply  put  m  =  0,  one  would  treat 
nothing  but  the  differential  equation 

k  JT  +  cx  =  0  (4.15) 

du 

which  by  aesuming  the  solution  to  be  of  the  form  x  =  A  ekk  is 
transformed  into  k  X  +  c  =  0,  whence  X,  =  — c/k.  That  is,  the 
solution  reads: 


24 


NACA  TM  No.  1217 


-  £t 

x  =  A  e  k  (4.16) 

However,  the  two  initial  conditions  x  =  0  and  dx/dt  =1  at  the 
time  t  =  0  cannot  be  satisfied  with  this  solution.  But  if  one  treats 
the  complete  differential  equation  (4.14)  in  the  same  manner  there  results 


+  k  X,  +  c  =  0 


and  hence: 


-  k  +  k 


2m 


or  the  square  root  might  be  developed  into  a  series  and  (since  now  the 

limiting  process  m - >0  is  to  be  performed)  broken  off  after  the 

second  term: 


Thus  X.-^  corresponds  to  the  previous  solution  of  the  firBt  order 
differential  equation,  where,  however,  X^  heen  lost.  For  very 

small  m,  X.^  ~  -k/m;  therewith  the  general  solution  becomes,  by 
combination  of  the  particular  solutions. 


X,i  t 

=  A^e  +  A^  e  =  A^ 


X-^t 


-  £t  -  -t 

e  k  +  A^  s  ® 


(4.17) 


Since  for  t  =  0,  x  is  also  supposed  to  equal  zero,  there  follows: 
A^  =  —  A-^,  thus: 


x  = 


(4.18) 


This  equation  is  plotted  schematically  in  figure  9*  The  first  term  of 
equation  (4.18),  which  alone  cannot  satisfy  the  boundary  conditions. 
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starts  from  the  value  at  the  time  t  =  0  and  decreases  expo¬ 

nentially.  The  second  term  is  important  only  for  small  t— values  and 
plays  no  role  for  large  t.  It  is  very  rapidly  variable  and  assures 
that  the  total  solution  (solid  line)  satisfies  the  boundary  conditions, 
"he  slowly  variable  solution  (in  corresponds  to  the  potential  flow, 

the  second,  rapidly  variable  particular  solution  (in  indicates,  as 

it  were,  the  narrow  region  of  the  boundary  layer;  the  smaller  m,  the 
narrower  this  region. 

Herewith  we  shall  conclude  the  general  remarks  and  turn  to  the  law 
of  similarity. 


CHAPTER  V.  REYNOLDS*  LAW  of  SIMILARITY 


So  far  no  general  methods  for  the  solution  of  the  Navier-Stokes 
differential  equations  are  known.  Solutions  that  are  valid  for  all 
values  of  the  viscosity  are  so  far  known  only  for  a  very  few  special 
cases  (for  instance,  Poiseuille*s  pipe  flow).  Meanwhile  the  problem 
of  flow  in  a  viscous  fluid  has  been  tackled  by  starting  from  the  limits, 
that  is,  one  has  treated  on  the  one  hand  flowB  of  very  great  viscosity, 

•  on  the  other  hand  flows  of  very  small  viscosity,  since  one  obtains  in 
this  manner  certain  mathematical  simplifications.  However,  starting 
from  these  limiting  cases  one  cannot  possibly  interpolate  for  flows  of 
average  viscosity. 

The  theoretical  treatment  of  the  limiting  cases  of  very  great  and 
very  small  viscosity  is  mathematically  still  very  difficult.  Thus 
research  on  viscous  fluids  was  undertaken  largely  from  the  experimental 
side.  The  Navier-Stokes  differential  equations  offer  very  useful  indi¬ 
cations,  which  permit  a  considerable  reduction  of  the  volume  of  experi¬ 
mental  investigation.  The  rules  in  question  are  the  so-called  laws  of 
similarity. 

The  problem  is:  Under  what  conditions  are  the  forms  of  flows  of 
any  liquids  or  gases  around  geometrically  similarly  shaped  bodies  them¬ 
selves  geometrically  similar^  Such  flows  are  called  mechanically  similar. 

Consider  for  instance  the  flows  of  two  different  fluids  of  different 
velocities  around  two  spheres  of  different  size  (fig.  10).  Under  what 
conditions  are  the  flows  geometrically  similar  to  each  other?  Obviously 
this  is  the  case  when  at  points  of  similar  position  in  the  two  flow 
patterns  the  forces  acting  on  volume  elements  at  these  points  have  the 
-  same  ratio.  Depending  on  what  kinds  of  forces  are  in  effect,  various 

laws  of  similarity  will  result  from  this  requirement. 

Most  important  for  this  investigation  is  the  case  where  all  forces 

*  except  the  inertia  and  friction  forces  are  negligible.  Furthermore,  no 
free  surfaces  are  to  be  present,  so  that  the  effect  of  gravity  is 
compensated  by  the  hydrostatic  pressure.  In  this  case  the  flow  around 
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the  two  spheres  is  geometrically  similar  when  the  inertia  and  friction 
forces  have  the  same  ratio  at  every  point. 


The  expressions  for  the  inertia  and  viscosity  forces  acting  on  the 
volume  element  will  now  he  derived:  there  is  as  friction  force  per 


unit  volume 


dT 


P  u 


6u 

dx 


The  ratio 


whereas  the  inertia  force  per  unit  volume  is 


inertia  force 


pu^ 

dx 


friction  force 


(5.1) 


must,  therefore,  he  the  same  at  all  points  of  the  flow.  One  now  inquires 
as  to  the  variation  of  these  forces  with  variation  in  the  quantities 
characteristic  of  the  phenomenon:  free  stream  velocity  V,  diameter  d, 
density  p,  and  viscosity  p.  For  variation  of  V  and  d  the  indi¬ 
vidual  quantities  in  equation  (5*1)  at  similarly  located  points  vary 
as  follows: 


u  *vV; 


du  Y.  a2u  V 
31  d  ay2  d2 


Therewith  equation  (5.1)  "becomes: 


Inertia  force 
friction  force 


(5-2) 


The  law  of  mechanical  similarity  is  therefore:  The  flows  around 
geometrically  similar  bodies  similarly  located  and  alined  with  respect 
to  the  flow  have,  for  equal  p  Y  d/p,  geometrically  similar  stream 
lines  as  well.  If  the  flows  in  question  are,  for  instance,  two  flows 
of  the  same  fluid  of  equal  temperature  and  density  (p  and  p  equal)  around 
two  spheres,  one  of  which  has  a  diameter  twice  that  of  the  other,  the  flows 
are  geometrically  similar  provided  that  the  free  stream  velocity  for  the 
larger  sphere  has  half  the  magnitude  of  that  for  the  smaller  sphere. 
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The  quantity  p  V  d/p  is, as  a  quotient  of  two  forces,  a  dimension¬ 
less  number.  This  fact  is  immediately  recognized  by  substituting  for 
the  quantities  their  dimensions: 


2  — U 
kg  sec  m 

m 

TTI 

kg  sec 

>  v 

sec 

>  a 

ill 

>  4 

2 

TTI 

—  - 

-  . 

pYd  _  kg  sec^  m  m^ 

4  nf^  sec  kg  sec 


This  law  of  similarity  was  discovered  by  Osborne  Reynolds  in  his 
studies  of  fluid  flows  in  a  pipe.  The  dimensionless  quantity  is  called 
after  him: 


pV  d/p  =  V  d/v  =  Re  =  Reynolds'  number 


The  introduction  of  this  dimensionless  quantity  helped  greatly  in 
advancing  the  development  of  modern  hydrodynamics. 


Connection  between  Similarity  and 
Dimensional  Considerations 

As  is  known,  all  physical  laws  can  be  expressed  in  a  form  free  of 
the  units  of  measure.  Thus  the  similarity  consideration  may  be  replaced 
by  a  dimensional  analysis.  The  following  quantities  appearing 
in  the  Navier— Stokes  differential  equations  are  essential  for  the  stream 
line  pattern:  V,  d,  p,  p.  The  question  is  whether  there  is  a  combination 


dp  p?  p& 

which  is  a  Reynolds  number  and  therefore  has  the  dimension  1.  This 
amounts  to  determining  a,  3,  7,  5  in  such  a  manner  that 

p?  =  K°  L°  T°  =  1  (5.3) 

with  K,  L,  T  representing  the  symbols  for  force,  length, and  time, 
respectively.  Without  limiting  the  generality  a  may  be  set  equal  to 
unity  (a  =  1)  since  any  power  of  a  dimensionless  quantity  is  still  a 
pure  number.  With  a  =  1  there  results  from  equation  (5*3) 
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[v4%ru6]  =  |L3  (n^B J 


=  K°  L°  T° 


(5.4) 


By  equating  the  exponents  of  L,  T,  K  on  the  left  and  right  sides  one 
obtains  the  three  equations: 


K: 

7  + 

8=0 

L:  1+0 

-ky  — 

25  =  0 

> 

(5-5) 

T: 

2  7  + 

6  =  1 

The  solution  gives: 


0  =  1;  7  =  1;  8  =  -1 


(5-6) 


Accordingly  the  only  possible  dimensionless  combination  of  V,  d,  p,  p 
is  the  quotient 


2S£  =  Ee 


(5.7) 


This  dimensional  analysis  lacks  the  pictorial  quality  of  the  similarity 
consideration;  however,  it  offers  the  advantage  of  applicability  even 
when  knowledge  of  the  exact  equation  of  motion  is  still  missing,  if 
there  is  only  known  what  physical  quantities  determine  the  phenomenon. 


CHAPTER  VI.  EXACT  SOLUTIONS  OF  THE  NAVIER— STOKES  EQUATIONS 


In  general,  the  problem  of  finding  exact  solutions  of  the  Navier- 
Stokes  differential  equations  encounters  insurmountable  difficulties, 
particularly  because  of  the  non-linearity  of  these  equations  which 
prohibits  application  of  the  principle  of  superposition.  Nevertheless 
one  can  give  exact  solutions  for  a  few  special  cases,  mostly,  when  the 
second  power  terms  vanish  automatically.  A  few  of  these  exact  solutions 
will  be  treated  here. 

One  investigates  first  layer  flows  in  general, _ that  is,  flows 
where  only  one  velocity  component  exists  which,  moreover,  is  not 
dependent  on  the'  analagous  position  coordinate,  whereas  the  two  other 
velocity  components  vanish  identically;  thus  for  instance: 
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u  =  u(y,  z,  t) 
v  =  w  s  0 


(6.1) 


The  Navier— Stokes  differential  equations  (3.16)  are  thereby  transformed 
into: 

\ 

(p  p  N 

Q  U  O  U 

-  +  - 

ay2  az2 


°  =  p  Y  - 

ay 


0  =  p  z  -  — 

az 


> 


(6.2) 


while  the  continuity  equation  is  identically  satisfied. 


a.  Pipe  Flow,  Steady  and  Starting 


1.  Steady  pipe  flow. 

For  the  case  of  the  pipe  lying  horizontally  the  mass  forces  are 
everywhere  constant,  and  the  equation  system  (equation  (6.2))  yields: 


(6.3) 


with  the  solution 


p=ax+by+cz 


Thus  the  pressure  is  a  linear  function  of  the  position. 

In  the  case  of  the  pipe  standing  vertically  (fig.  11)  the  mass 
force  is  constant  in  the  y—  and  z— direction  and  increases  in  the 
x— direction  corresponding  to  the  hydrostatic  pressure  of  the  field  of 
gravity.  Moreover,  if  one  puts  Y  =  0  and  Z  =  0: 


dp/dy  -  dp/dz  =0:  v  =  w  =  0 
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Thus  there  remains  from  equation  (6.2) 


Sfp 

ax2 


=  0  or 


—  =  constant 
cfcx 


The  system  (equation  (6.2))  uhen  becomes,  under  the  further  assumption 
of  steady  flow: 


&  +  *?“  =  i  &  =  constant  (6.1* ) 

.2  2  H  &c 

dy  dz 


One  tries  a  solution  of  the  form 


u  =  u  (a2*  by2  +  cz2) 
max 


(6.5) 


with  Ujnax  representing  the  velocity  at  the  pipe  center.  With  the 
condition  u  =  0  at  the  pipe  wall  y2  +  z2  =  r2  there  results: 

a  =  1;  b  =  c  =  — l/r2 


and  therefore 


with 


u(y,  z) 


^^max 


"N 


> 

r2  6p 
4p.  dx 

) 


(6.6) 


This  solution  is  identical  with  equations  (2.3)  and  (2.4)  if  one  bears 
in  mind  that  in  those  equations  y  represented  the  radial  distance  from 
the  pipe  center.  Thus  Poiseuille'B  pipe  flow  was  found  as  an  exact 
solution  of  the  Navier— Stokes  differential  equations. 


2.  Starting  flow  through  a  pipe. 

By  the  expression,  "Starting  flow  through  a  pipe,"  the  following 
problem  is  meant:  Let  the  fluid  in  a  circular  pipe  of  infinite  length 
be  at  rest  until  the  time  t  =  0.  At  the  time  t  =  0  let  a  pressure 
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difference  invariable  with  time,  —  p^,  suddenly  be  established 

across  a  piece  of  pipe  of  length  2,  so  that  for  the  entire  pipe  the 
following  equation  applies: 


t  >  0: 


Sp 

dx 


=  constant 


(6.7) 


Here  —  /  0, 
St 


and  the  inertia  and  friction  terms  balance  each  other. 


The  solution  of  this  problem,  which  will  not  be  further  discussed  here, 
was  given  by  Szymanski  (reference  12).  The  phenomenon  is  not  dependent 
on  the  longitudinal  coordinate  x.  The  velocity  profiles  at  various  times 
can  be  seen  from  fig.  12.  It  is  characteristic  that,  first,  at  the 
pipe  center,  the  velocity  remains  locally  constant  and  the  friction  is 
noticeable  only  in  a  thin  layer  near  the  surface.  Only  later  does  the 
friction  effect  reach  the  pipe  center.  Poiseuille*8  parabolic  profile 
of  the  steady  pipe  flow  is  attained  asymptotically  for  t->  »  . 


One  must  clearly  distinguish  between  the  non— steady  starting  flow 
through  a  pipe  discussed  here  and  the  steady  pipe  inlet  flow.  This 
latter  is  the  flow  at  the  inlet  of  a  pipe.  The  rectangular  velocity 
profile  present  in  the  entrance  cross  section  is,  with  increasing 
distance  x  from  the  inlet,  gradually  transformed  under  the  effect  of 

friction  into  Poiseuille*s  parabolic  profile.  Since  here  /  0, 

dx 

this  is  not  a  layer  flow.  This  inlet  flow  was,  for  the  plane  problem, 
exactly  calculated  from  the  differential  equations  by  H.  Schlichting 
(reference  14)  and  for  the  rotationally-symmetrical  problem,  according 
xo  an  approximate  method,  by  L.  Schiller  (reference  32). 


Third  lecture  (Dec.  15,  19^1) 
b.  Plane  Surface;  a  Surface  Suddenly  Set  in  Motion 
and  an  Oscillating  Surface 
1.  Oscillating  Surface. 

Let  a  plane  surface  of  infinite  extent  perform  in  its  plane  recti¬ 
linear  oscillations  in  the  x— direction  (fig.  13)«  The  y-direction  is 
assumed  normal  to  the  surface.  Let  the  oscillation  take  place  with  the 
velocity  uQ  =  A  cos  n  t,  with  A  denoting  the  amplitude,  n  the 

frequency  of  the  oscillation.  The  fluid  near  the  surface  is  carried  along 
by  the  friction.  Since  the  surface  is  of  infinite  extent,  the  state  of 
flow  is  Independent  of  x  and  z.  The  flow  in  question  is  therefore  a 
non— steady  plane  layer  flow  for  which 
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u  =  u(y,t); 


v  =  v  =  0; 


dp 

8x 


=  0 


If,  finally,  one  puts  the  mass  forces  equal  to  zero,  the  Navier-Stokes 
differential  equations  (equation  (3.16))  are  reduced  to 


£&  =  y 

^t  .  2 
dy 


(6.8) 


with  the  boundary  condition  y  =  0,  u  =  uQ  cos  n  t.  This  equation  has 

the  same  structure  as  the  differential  equation  for  thermal  expansion  in 
a  rod  (linear  equation  of  heat  conduction).  It  becomes  identical  with 
it  if  v  is  replaced  by  the  thermal  conductivity  a  and  u  by  the 
temperature  (compare  equation  (4.12)). 

For  the  solution  of  equation  (6.8)  one  uses 


u  =  Ae^y 


cos  (nt  -  ky) 


(6.9) 


where  k  is  a  constant  to  be  determined. 
Then 


=  —  n  Ae~ty  sin  (nt  —  ky) 
ot 


=  —  k  Ae^y  cos  (nt  —  ky)  +  k  Ae-^  sin  (nt  —  ky) 
dy 


Tv 


— -  =  —  2k2Ae  sin  (nt  -  ky) 

8y2 


Insertion  in  equation  (6.8)  gives  nA  =  2k^AV,  and  thence 

k  =  \! -2-  (6.10) 

\  2  v 


Thr  velocity  distribution  u(y,  t)  is  therefore  an  oscillation  with 

_  ILpr 

amplituie  decreasing  toward  the  outside  Ae  4  ;  the  layer  at  the 
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distance  y  from  the  surface  has  a  phase  lag  of  y  y  with  respect  to 

the  motion  of  the  surface.  The  wave  length  X.  in  the  y-direction 
follows  from  k\  =  2n  as 


\ 


(6.11) 


Thus  the  co-oscillating  layer  is  thinner ,  the  greater  the  frequency  n 
and  the  smaller  the  viscosity  V.  The  result  \  ~  \[v  is  to  be  noted. 
The  velocity  profiles  for  various  times  are  given  in  figure  13. 

2.  The  Surface  Suddenly  Set  in  Motion. 

Equation  (6.8)  yields  another  exact  solution  of  the  Navier— Stokes 
differential  equations,  namely  the  flow  in  the  neighborhood  of  a  plane 
surface  which  suddenly  starts  moving  in  its  own  plane  with  the  constant 
velocity  uQ.  In  this  case  the  boundary  conditions  are 


t  —  0,  for  all 
t  >  0,  y  =  0: 


(6.12) 


An  appropriate  variable  for  the  solution  is  the  dimensionless  quantity 


T)  = 


.7 

s\Jvi 


(6.13) 


As  the  solution  of  equation  (6.8)  with  the  boundary  conditions 
equation  (6.12)  one  obtains 


u(y,  t) 


(6.14) 


The  correctness  of  the  solution  is  readily  confirmed  by  substitution. 
The  variation  of  the  flow  with  time  is  indicated  in  figure  14. 

The  probability  integral  appearing  in  equation  (6.14)  — - 

Jn 

has,  for  q  =  q  =1.9,  the  value  0.01.  Therefore,  for  q  -  qK, 

u  =  0.01  uQ.  The  thickness  of  the  layer  carried  along  by  the  friction 
is,  therefore. 


3^ 
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&  =  2  tjb  \'vt  (6.15) 


Thus  here  again  the  friction  layer  thickness  is 

5  ~ 

At  this  point  we  conclude  the  discussion  of  layer  flows  and  turn  to  a 
few  other  exact  solutions  of  Navier-Stokes  differential  equations. 


c.  Plane  Stagnation— Point  Flow 

The  plane  flow  in  the  neighborhood  of  a  stagnation  point  on  a  smooth 
wall  is  considered.  With  the  cosrdinate  system  according  to  figure  15 
the  corresponding  potential  flow  has  the  potential 

®  =  |  (x2  -  y2)  (6.i6) 


and  the  stream  function 


t  =  a  x  y 


The  velocity  components  are: 

U  =  ax;  V  =  -ay  (a  =  constant)  (6.17) 


This  is  a  potential  flow  which,  coming  from  the  direction  of  the  y-axis, 
encounters  the  solid  wall  y  =  0,  divides,  and  flows  off  parallel  to 
the  x— axis.  Whereas  the  potential  flow  glides  along  the  wall,  the  viscous 
flow  must  adhere  to  it.  If  one  deBigiates  the  velocity  components  of 
the  viscous  flow  by  u(x,  y)  and  v(x,  y),  the  boundary  conditions  for 
them  are: 


y=0;  u=0;  v=0 

y  =  °°;  u  =  U;  v  =  V 


(6.18) 


For  the  stream  function  of  the  viscouB  fluid  one  uses  the  equation: 
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3^ 


1 1  (x,  y)  -  x  f(y) 


(6.19) 


Thence  one  obtains: 


u  =  x  —  =  x  f»(y) 

dy 


=  -  r(v) 


(6.20) 


The  boundary  conditions  (Equation  (6.18))  require 


y  =  0:  r  =  f*  =  0 
y  =  “ :  f  =  ay;  f  *  =  a 


(6.21) 


At  large  distances  from  the  surface,  that  1b,  in  the  potential  flow,  the 
pressure  is  calculated  from  Bernoulli* s  equation 

p  =  PQ  ~  |  W2  =  p0  -  |  a?  (x2  +  y 2)  (6.22) 

with  W  =  V 2  signifying  the  magnitude  of  the  velocity  of  the 

potential  flow,  given  by  equation  (6.17).  Let  pQ  be  the  total  pressure 
of  the  potential  flow. 

For  the  viscous  flow  one  formulates  the  analogous  equation 


p  =  pQ  -  |  a2  [x2  +  F(y)] 


(6.23) 


With  these  equations  one  turns  to  the  Navier-Stokes  equations  which  read, 
for  vanishing  mass  forces: 


u 


u 


8u 

—  _L 

8u 

V  —  — 

1 

*E  +  J 

(  v2 

a  u 

A. 

a2u 

a y 

p 

ax  \ 

i- 

CVi 

ay2 

+ 

v  = 

_  1 

^  +  V  1 

, 

a2v' 

ax 

ay 

p 

ay  ' 

I  2  + 

2 

\ax 

ay  - 

(6.24) 


36 


NACA  TM  No.  1217 


The  continuity  equation  has  already  been  integrated  by  introduction  of 
the  stream  function  By  substitution  of  equations  (6.20)  and  (6.23) 
into  equation  (6.24)  one  obtains 


x(f«2  _  f  f")  =  a2x  +  vrf»" 

f  ft  =  F»(y)  -  v  f" 

2 


(6.25) 


These  are  two  differential  equations  for  the  two  unknown  functions 
f(y)  and  F(y)  which  determine  the  velocity  distribution  and  the 
pressure  distribution,  respectively.  The  component  of  the  velocity  u 
parallel  to  the  surface,  that  is,  the  function  f’(y),  is  of  particular 
interest.  Since  F  (y)  does  not  appear  in  the  first  equation,  one 
solves  first  the  first  equation  and  then,  after  substitution  of  this 
solution,  the  second.  Thus  the  differential  equation  to  be  solved  reads 
at  first 


f2- 


f  f"  =  ar  +  v  f,M 


(6.26) 


with  the  boundary  conditions  according  to  equation  (6.21). 

This  non— linear  differential  equation  cannot  be  solved  in  closed 
form.  If  one  introduces  instead  of  y  the  variable 


\  =  a  y 


( 6 . 26a ) 


and  in  addition  the  similarity  transformation 


f(y)  =  Aip(S) 


( 6 . 26b ) 


Then  the  inhomogeneous  term  in  equation  (6.26)  becomes  equal  to  1,  and 
the  solutions  therewith  become  independent  of  the  specific  data  for  the 
flow.  Thereby  equation  (6.2 6)  becomes 


a2A2 


dqA2  ___ 

\d  %) 


<P 


d2<p 

I? 


.  a2  ♦  vta3  12 
di3 


If  one  now  equates  a2 A2  =  a2  and  vAa^  =  a2,  that  is. 


(6.27) 
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(6.27a) 


the  differential  equation  for  <p(£)  with  i  =  y  reads: 

<p»"  +  <p"  =  cp*2  +1  =  0  (6.28) 


with  the  boundary  conditions 


|  =  0:  cp  =  q>*  =  0 

|  =  »:  cp*  =  1 


The  solution  found  by  series  development  can  be  found  in  the  thesis  of 
Hiemenz  (reference  10),  compare  table  1*.  The  velocity  component 
parallel  to  the  surface  is 

5  =  1  f*(y)  =  <p’U)  (6.29) 

U  a 


It  is  indicated  in  figure  l6.  The  curve  cp’  (  I)  increases  linearly  at 
1  =  0  and  approaches  one  asymptotically.  For  about  |  =  2.6,  <p*  ~  0.99; 

thus  within  about  one  percent  of  the  final  value.  If  one  again  designates 
the  corresponding  distance  from  the  surface  y  =  6  as  the  boundary  layer 
thickness  (friction  layer  thickness),  then 

6  =  I  h  =  S-6H  (6-30) 


Thus  in  this  flow,  as  in  the  former  ones, 

5  ~  \Jv 


It  is  also  remarkable  that  the  dimensionless  velocity  distribution 
according  to  equation  (6.29)  and  the  boundary  layer  thickness  according 
to  equation  (6. 30)  are  independent  of  x,  thus  do  not  vary  along  the  wall. 


*The  tables  appear  in  appendix,  chapter  XII. 
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For  later  applications  the  characteristics  important  for  the  friction 
layer,  displacement  thickness  5*  and  momentum  thickness  d,  are 

introduced  here;  they  are  defined  by 


US*  =  f  (U  -  u)  dy  (6.31) 

y=o 


2  P00 

U  d  =  J  u(TJ  -  u)  dy  (6.32) 

7=0 


The  displacement  thickness  gives  the  deflection  of  the  stream  lines  of 
the  potential  flow  from  the  surface  by  the  friction  layer;  the  momentum 
thickness  is  a  measure  of  the  momentum  loss  in  the  friction  layer.  By 
insertion  of  equation  (6.29)  in  (6.31)  and-  (6.32)  and  calculation  of  the 
definite  integral  one  finds 


(6.33) 


(6.34) 


and  hence 


-*=  2.218 

d 


(6.35) 


The  quantity  S*  is  indicated  in  figure  l6.  For  comparison  with  a 
later  approximate  solution  one  also  notes  the  numerical  value  of  the 


dimensionless  quantity 


5»2  dU 
V  dx* 


One  finds  from  equations  (6.17)  and  (6.33) 


—  22  =  0.4202  (6.36) 

v  dx 


The  exact  solution  of  the  Navier— Stokes  differential  equations  found 
here  gives,  therefore, • for  large  Reynolds  numbers  a  friction  layer 

thickness  decreasing  with  \/~  and  a  transverse  pressure  gradient 
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decreasing  with  pa\/va.  Both  confirm  the  boundary  layer  assumptions* 
to  be  discussed  later. 


d.  Convergent  and  Divergent  Channel 

A  further  class  of  exact  solutions  of  the  Navier— Stokes  differential 
equations  exists  for  the  convergent  and  divergent  channel  with  plane  walls 
(fig.  17),  as  given  by  G.  Hamel  (reference  11). 

Without  entering  into  the  details  of  the  rather  complicated  calcu¬ 
lations  the  character  of  the  solutions  will  be  briefly  sketched: 

The  velocity  distributions  for  convergent  channels,  plotted  against 
distance  along  the  surface  for  various  included  angles  a  and  for  various 
Re-numbers  appear  as  indicated  in  figure  18.  At  the  tunnel  center  the 
velocity  is  almost  constant,  and  at  the  surfaces  it  suddenly  declines  to  zero. 

In  the  case  of  divergent  tunnels  one  obtains  greatly  differing 
forms  for  the  velocity  profiles,  depending  on  the  included  angle  and 
the  Re-number.  Here  all  velocity  profiles  have  two  inflection  points. 

For  small  Re-numbers  and  small  included  angles  the  velocity  is  positive 
over  the  entire  cross  section  (solid  curve  in  fig.  19);  for  larger  angles 
and  larger  Re-numbers,  on  the  other  hand,  the  velocity  profiles  have 
reverse  flow  at  the  surface  (dashed  curve  in  fig.  19).  The  reverse  flow 
is  the  initial  phase  of  a  vortex  formation  and  therefore  of  the  separation 
of  the  flow  from  the  surface.  Generally,  the  separation  does  not  occur 
symmetrically  on  both  surfaces;  the  flow  separates  from  one  side  and 
adheres  to  the  other  surface  (fig.  20). 

These  examples  also  confirm  the  theory  that  exact  solutions  have  the  same 
character  as  approximate  solutions  of  boundary  layer  theory;  in  particular, 
they  confirm  that  for  the  convergent  channel  a  very  thin  layer  with  con¬ 
siderable  friction  effect  is  present  near  the  surface  (here  also  the 
calculation  shows  that  the  layer  thickness  *v  ^y)  and  that  for  the 
divergent  channel  reverse  flow  and  separation  occur. 

We  here  conclude  the  chapter  on  the  exact  solutions  of  the  Navier— 

Stokes  differential  equations  and  turn  to  the  approximate  solutions. 

By  exact  solutions  have  been  understood  those  where  in  the  Navier-Stokes 
differential  equations  all  terms  are  taken  into  consideration  that,  in 
the  various  cases,  are  not  identically  zero.  By  approximate  solutions 
of  the  Navier— Stokes  differential  equations  will  be  understood,  in 
contrast,  solutions  where  terms  of  small  magnitude  are  neglected  in  the 
differential  equations  themselves.  However,  by  no  means  are  all  the 
friction  terms  to  be  neglected  simultaneously,  since  this  would  represent 
the  case  of  potential  flow. 


*The  rotationally— symmetrical  stagnation— point  flow  has  been 
calculated  by  Homann  (reference  17).  Instead  of  equation  (6.28)  one 

obtains  the  differential  equation  cp1"  +  2cpcp"  —  cp'^  +  1  =  0. 
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CHAPTER  VH.  VERY  SLOW  MOTION  (STOKES,  OSEM) 


The  exact  solutions  of  the  Navier-Stokes  differential  equations 
discussed  in  the  previous  chapter  are  of  a  very  special  kind.  Most  of 
them  dealt  with  flows  along  a  plane  surface,  where  the  stream  lines 
are  rectilinear.  Most  flows  existing  in  practice,  as  for  instance  flows 
around  arbitrary  bodies,  cannot  be  calculated  exactly  from  the  Navier— 
Stokes  differential  equations,  but  must  be  treated  by  approximate  methods. 
Two  kinds  of  such  approximations  are  possible: 

1.  For  predominant  viscosity,  completely  neglecting  the  inertia  terms 
suggests  itself  (very  small  Re-number;  He <  1). 


2.  For  very  small  viscosity  and  therefore  predominant  inertia  one 
takes  the  viscosity  into  consideration  only  in  a  very  thin 
layer  in  the  neighborhood  of  the  solid  wall;  for  the  rest, 
the  flow  is  regarded  as  frictionless.  Here  the  Re-number  is 
very  large  (Prandtl*s  boundary  layer  theory). 

The  first  limiting  case  with  very  small  Re-number  will  be  discussed 
in  this  chapter.  A  small  Re-number  indicates  small  velocities,  small 
body  dimensions,  and  large  viscosity.  Since  the  inertia  terms  depend 
on  the  square  of  the  velocity  whereas  the  friction  terms  are  linear, 
all  inertia  terms  in  the  Navier-Stokes  differential  equations  are,  for 
very  small  Re-numbers,  negligible.  It  is  to  be  expected  that  an 
approximation  will  thereby  be  obtained  for  very  slow  (creeping)  motion, 
as  for  instance  the  falling  of  a  minute  fog  particle*)  or  the  slow 
motion  of  a  body  in  a  very  viscid  oil. 


Neglecting  all  inertia  terms  one  obtains  from  the  Navier-Stokes 
differential  equations  (3.16)  the  following: 


dx 


nAul 


=  mav  > 

=  HAW 

dz 

du  |  fry  |  dw 
dx  dy  dz 


(7-1) 


(7.2) 


*For  a  sphere  falling  in  air  (v  *  14  x  10“^  m^/sec)  for  instance: 
Re  =  V  d/v  =  1,  for  d  =  1  mm;  V  =  l.tO  cm/sec. 
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The  same  boundary  conditions  apply  to  this  system  of  equations  as  apply 
to  the  complete  Navier— Stokes  differential  equations,  namely  vanishing 
of  the  normal  component  vn  -  0  and  the  tangential  component  v^_  =  0 

at  the  hounding  surfaces. 

The  neglect  of  all  inertia  terms  in  Navier-Stokee  differential 
equations  does  not  represent  as  serious  an  inaccuracy  as  the  neglect  of 
all  friction  terms  when  transforming  the  Navier— Stokes  differential 
equations  into  Euler’s  differential  equations  of  the  frictionless  flow. 
That  is,  by  neglecting  the  inertia  terms,  the  order  of  the  differential 
equations  is  not  lowered  so  that  in  the  simplified  differential  equations 
the  same,  boundary  conditions  as  in  the  Navier— Stokes  complete  differential 
equations  can  still  be  satisfied. 

Furthermore  one  obtains  from  the  equations  (7*1),  taking  into 
account  the  continuity,  by  differentiating  the  first  with  respect  to  x, 
the  second  with  respect  to  y,  the  third  with  respect  to  z,  the  following 
equation  for  the  pressure  p 

=  =  °  (7.3) 

chc  dy  dz 


that  is,  for  creeping  motions  the  pressure  function  p(x,  y,  z)  is  a 
potential  function. 

The  details  of  the  calculation  will  not  be  discussed  more  thoroughly, 
particularly  since  the  creeping  motion  is  technically  not  very  important. 
However,  at  least  Stokes*  famous  solution  for  the  sphere  will  be  discussed 
briefly  (fig.  21).  The  drag  of  a  sphere  for  creeping  motion  consists  of 
the  contributions  of  the  pressure  drag  (form  drag)  and  the  surface  fric¬ 
tion  drag.  The  latter  is  obtained  by  integration  of  the  wall  shearing 
stress  over  the  entire  sphere  surface.  Stokes  performed  the  integration 
of  the  equation  systems  (7.1)  and  (7.2)  for  a  sphere  in  a  uniform 
flow  of  velocity  UQ.  There  results,  according  to  Stokes,  for  the  entire 

drag  of  the  sphere  of  radius  R: 

W  =  W_  +  W_  =  6  itpR  U  (7.4) 

Dr  R  o 


The  drag  is,  therefore,  proportional  to  the  first  power  of  the  velocity. 

If  one  introduces  for  the  sphere  a  drag  coefficient  cw  which,  in  the 
customary  manner,  is  referred  to  the  frontal  area  and  the  dynamic  pressure 
of  the  free  stream  velocity 
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W  =  cw 


n  R'  %  Ur 


(7.5) 


there  results  for  the  drag  coefficient  according  to  Stokes'  formula: 


24  „  UoD 

‘*■5?  Ee=  — 


(7.6) 


One  can  state  immediately  that  the  stream  line  pattern  of  this 
creeping  motion  must  be  the  same  ahead  of  and  behind  the  sphere  since 
for  reversal  of  the  initial  flow  (sign  reversal  of  the  velocity 
components)  the  equation  system  (equation  (7.1))  goes  over  into  itself. 
The  stream— line  pattern  for  the  viscous  sphere  flow,  as  it  presents 
itself  to  an  observer  who  is  at  rest  relative  to  the  flow  at  infinity, 
is  shown  in  figure  22.  The  fluid  particles  are  pushed  aside  by  the 
sphere  in  front  and  come  together  again  behind  it. 

As  shown  by  a  comparison  of  Stokes'  drag  formula  equation  (7.6) 
with  test  results  (reference  33),  this  formula  is  valid  only  for  the 
region  Re  <  1. 


Correction  by  Oseen 

In  Oseen' s  later  improvement  of  Stokes'  solution  for  the  sphere 
the  inertia  terms  in  the  differential  equations  are  partly  taken  into 
consideration.  Oseen  formulates  the  velocity  components  u,  v,  w: 


u  =  UQ  +  u';  v  =  v*;  w  =  w' 


(7.7) 


where  u' ,  v' ,  w'  may  be  considered  as  disturbance  velocities  which  in 
general  are  small  compared  with  the  free  stream  velocity  UD.  This 

assumption  is  not  actually  correct  for  the  immediate  proximity  of  the 
sphere  surface.  With  the  formulation  (equation  (7*7))  the  inertia  terms 
in  equation  (3.16)  are  divided  into  two  groups,  for  instance: 


du*  TT  dvf 
o  t  >  uo  t  > 
ox  ox 


and 


The  second  group  of  second  order,  as  compared  with  the  first  group,  is 
neglected.  Therewith  one  then  obtains  "rom  the  Navies— Stokes  equations 
of  motion  (7,.l6)  the  following  equations  of  motion,  which  are  taken  as 
a  basis  by  Oseon. 
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^3 


p  UQ  +  &  =  uAu* 
ox  ox 


TT  ^V*  ,  At  l 

p  u°  &T  +  *  -  > 


TT  <3w*  dp  .  . 

p  U  -  +  —  =  uAw* 

0  dx  dz 


(7.8) 


In  addition,  one  usee  the  continuity  equation: 


du*  dv1  +  dv1  _  q 

dx  dy  d  z 


(7.9) 


and  the  same  boundary  conditions  as  in  the  Navier-Stokes  differential 
equations.  One  calls  the  contributions  of  the  convective  terms  in  these 
equations  that  were  taken  into  consideration,  for  Instance  UQ  ^1,  th? 

dx 


semi -quadratic  terms.  These  differential  equations  of  Oseen  and  Stokes1 
differential  equations  are  both  linear.  The  stream  line  pattern,  as  it 
results  for  this  sphere  flow  according  to  Oseen,  is  given  In  figure  23. 
Here  again  the  observer  is  at  rest  relative  to  the  fluid  at  large  distance 
from  the  sphere.  Thus  the  sphere  is  dragged  past  the  observer  with  the 
velocity  Up.  The  stream  line  pattern  ahead  of  and  behind  the  sphere  are 


now  not  the  same,  as  was  the  case  in  Stokes*  solution.  Ahead  of  the 
sphere  exists  almost  the  same  displacement  flow  as  in  Stokes’  pattern; 
behind  the  sphere,  however,  the  stream  lines  are  closer  together,  that  is 
the  velocity  is  greater  here  than  in  Stokes'  case.  A  wake  is  present 
behind  the  sphere  similar  to  that  from  test  results  for  large  Reynolds 
numbers. 


For  the  sphere  drag  calculated  by  Stokes  there  results  with  the 
drag  coefficient  cw  introduced  in  equation  (7*5)  the  formula: 


c 


w  ~ 


2k 

Re 


(1  + 


-2-  Re 

l6 


) 


Re 


(7.10) 


The  test  results  (reference  33)  show  that  Oseen* s  formula  is  fairly 
accurate  up  to  about  Re  =  5* 

With  these  brief  remarks  we  conclude  the  limiting  case  of  small 
Reynolds  numbers  and  turn  to  the  case  which  is  of  foremost  interest  in 
practice:  the  case  of  very  large  Reynolds  number. 


CHAPTER  VIII.  FRANDTL’S  BOUNDARY  LAYER  EQUATIONS 


The  other  extreme  case  of  very  small  viscosity  or  of  very  large 
Reynolds  number  will  now  be  treated.  In  this  case  the  inertia  effects 
are  predominant  within  the  main  body  of  the  fluid  whereas  the  viscosity 
effects  there  are  almost  negligible. 

A  significant  advance  in  the  treatment  of  motion  of  fluids  for 
large  Reynolds  numbers,  that  is,  in  general,  of  fluids  of  very  small 
viscosity,  was  attained  by  L.  Prandtl  in  1904  (reference  7).  Prandtl 
demonstrated  in  what  way  viscosity  is  essential  for  large  Reynolds 
numbers  and  how  one  can  simplify  the  Navier— Stokes  differential  equations 
in  order  to  obtain  at  least  approximate  solutions. 

Let  us  consider  the  motion  of  fluid  of  very  small  viscosity,  for 
instance  of  air  or  water  surrounding  a  cylindrical  streamline  body 
(fig.  24).  Up  to  very  near  the  surface  the  velocities  are  of  the  order 
of  magnitude  of  the  free  stream  velocity  UQ.  The  stream  line  pattern 

as  well  as  the  velocity  distribution  agree  to  a  large  extent  with  those 
of  the  frictionless  fluid  (potential  flow).  More  thorough  investigations 
show,  however,  that  the  fluid  by  no  means  glides  along  the  surface  (a'5 
in  potential  flow)  but  adheres  to  it.  The  transition  from  zero  velocity 
at  the  surface  to  the  fully  developed  velocity  as  it  exists  at  some 
distance  from  the  body,  is  effected  in  a  very  thin  layer.  Thus  one  must 
distinguish  between  two  regions  which,  it  is  true,  cannot  be  rigorously 
separated : 

1.  A  thin  layer  in  the  immediate  proximity  of  the  body  where  the 

velocity  gradient  normal  to  the  surface  ^  is  very  large 

dn 

(boundary  layer).  Here  the  viscosity  4,  though  very  small, 
plays  an  essential  role  inasmuch  as  the  frictional  shearing 
&w 

stress  t  =  u  —  can  assume  considerable  values, 
on 

2.  In  the  remaining  region  outside  of  this  layer  velocity  gradients 

of  such  magnitude  do  not  occur,  so  that  there  the  effect  of 
viscosity  becomes  insignificant.  Here  frictionless  potential 
flow  prevails. 

In  general  one  may  say  that  the  boundary  layer  is  thinner,  the 
smaller  the  viscosity  or,  more  generally,  the  larger  the  Re-number.  It 
was  shown  before  on  the  basis  of  exact  solutions  of  the  Navier— Stokes 
differential  equations  that  the  boundary  layer  thickness  Is 

6  ~ 

The  approximations  to  the  Navier— Stokes  differential  equations  to  be  made 
below  are  more  valid  the  thinner  the  boundary  layer.  Thus  the  solutions 
of  the  boundary  layer  equations  have  an  asymptotic  character  for  infinitely 
increasing  Reynolds  numbers. 
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Let  us  now  make  the*  simplifications  of  the  Navler-Stokes  differential 
equations  for  the  boundary  layer.  To  this  end  the  order  of  magnitude  of 
the  separate  terms  of  the  Navler-Stokes  differential  equations  must  he 
estimated.  One  considers  the  flow  around  a  cylindrical  body  according 
to  figure  24.  One  imagines  the  Navler-Stokes  differential  equations 
written  non-dimensionally,  by  referring  all  velocities  to  the  free  stream 
velocity  U0  and  the  lengths  to  a  body  length  2.  The  pressure  will  be 
made  dimensionless  with  p  U02,  the  time  with  2/UQ.  Furthermore 

Re  =  -2-  represents  the  Reynolds  number.  Accordingly,  the  Navler-Stokes 

differential  equations  become  —  omitting  the  mass  forces  according  to 
equation  (k.2),  by  writing  the  same  letters  for  the  dimensionless  quantities 
as  for  the  dimensional  ones  — 


du  du  du 

—  +  u  —  +  v  —  = 

dt  bx  by 

1  11  5  i 

8 


bv 

dv 

dv 

—  + 

u  —  + 

V  —  s 

bt 

dx 

dy 

8 

1  8 

8  1 

dp  +  l 

a2u 

a2u 

dx  R 

dy2 

82 

i 

l/8‘ 

dp  1 

a2v 

by  R 

a/ 

dy2 

82 

8 

l 

8 

(8.1) 


bu.  [  dv 
dx  by 
1  1 


0 


The  estimation  gives:  Longitidinal  velocity  u  is  of  the  order  of 
magnitude  1.  Dimensionless  boundary  layer  thickness  B/z«l.  Therefrom 
follows : 


du  _1  _  d^u  ^  _1_ 

*  5’  3y2  62 


whereas  the  derivatives  with  respect  to  x  are  of  normal  order  of 
magnitude,  thus 

o 

du  , ,  d  u 

- —  i;  — 

ax2 


1 


be 
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From  the  continuity  equation  follows  therefore:* 


dv 

by 


1: 


v 


dv 

dy 


6 


The  transverse  velocity  v  in  the  boundary  layer  is  therefore,  to  the 
first  order,  small  in  comparison  with  the  longitudinal  velocity.  Further 
there  follows: 


dx 


S; 


/w 


1 

5 


Thus  there  result  for  the  single  terms  of  Navier— Stokes  differential 
equations  the  orders  of  magnitude  noted  in  equation  (8.1):  On  the  right 
side  in  the  first  equation 


so  that  it  can  be  neglected 

One  now  obtains  within  the  boundary  layer  friction  terms  which  are 
of  the  same  order  of  magnitude  as  the  inertia  terms,  if  l/R  is  of  the 
order  of  magnitude  5^,  or,  if  the  dimensional  quantities  are  again 
written. 


6 

Z 


(8.2) 


In  the  second  equation  of  motion  all  terms  then  are  of  the  order  of 

rta 

magnitude  5,  including  the  transverse  pressure  gradient  — •  In  the 

dy 

boundary  layer,  as  long  as  it  is  thin,  the  dependence  of  the  pressure 
on  y  may  therefore  be  neglected.  Thus  approximately  the  same  pressure 
prevails  within  the  boundary  layer  as  at  its  edge,  that  is,  the  pressure 
of  the  potential  flow.  The  pressure  within  the  boundary  layer  is  therefore, 
as  it  were* impressed  by  the  potential  flow. 


The  second  equation  of  motion  is  therewith  exhausted  and  does  not 
have  to  be  considered  further.  Using  again  the  dimensional  quantities 
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one  now  makes  the  Navier— Stokes  differential  equations  assume  the  following 
simplified  form: 


— 

St 


Su 

+  u  —  +  v 
Sx 


Su 

dy 


i  jp  +  y  a2u 

p  a*  3y2 


Sx  Sy 

Boundary  conditions:  y  =  0:  u=v=0 

y  =  oo  :  u  =  U 


> 


(8.3) 


The  pressure  in  the  Boundary  layer,  which  is  dependent  only  on  x,  is 
determined  from  the  potential  flow  U  (x,  t),  assumed  to  be  known 
according  to  Bernoulli’s  equation. 


or 


=  p  +  £  U  2 
o  2  ° 


l  S£  =  n  3u  +  Su 

p  Sx  Sx  St 


for  stationary  flow  j 


>  (8-4) 


for  nonstationary  flow 


/ 


With  the  potential  flow  known,  the  equation  system  (8.3)  represents 
a  system  of  two  equations  with  the  two  unknowns  u  and  v. 


Numerical  example:  In  order  to  help  clarify  the  concepts  a  numerical 
example  is  given  for  the  thickness  of  the  boundary  layer.  The  problem  is: 
What  is  the  boundary  layer  thickness,  for  instance,  in  the  case  of  the 
plate  in  longitudinal  flow  at  the  distance  l  =  100  centimeters  from  the 
leading  edge?  Let  the  velocity  be  uQ  =  30  meters  per  second  and  the 

kinematic  'viscosity  for  air  V  -  0.14  X  10~^  meters  square  per  second; 

then  the  Reynolds  number  is  E  =  UQ  |  =  2.1  X  10^  and  -Tb  =  1.45  X  10^, 

A  numerical  factor  is  still  missing  in  the  formula  (8.2)  for  the  boundary 
layer  thickness.  For  the  plate  in  longitudinal  flow  it  is,  as  later 
calculations  will  Bhow,  five,  provided  one  understands  by  the  boundary 
thickness  S  that  distance  from  the  surface  where  the  velocity  has  the 
value  0.99  uQ.  Thus  a  calculation  by  the  formula 
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results  for  the  present  case  in  a  boundary  layer  thickness  5  =  3.45  mm. 

It  should  he  added  that  the  Reynolds  number  is  already  so  large  that  the 
boundary  layer  at  the  end  of  the  plate  would  be  turbulent.  The  transition 
from  laminar  to  turbulent  lies  further  upstream,  and  at  this  point  the 
boundary  layer  thickness  would  then  be  somewhat  smaller  than  the  value 
determined  above. 


Fourth  lecture  on  December  22,  1942 


Physical  Summary  and  Conclusions 


The  physical  content  of  the  considerations  so  far  can  be  summed  up 
in  the  following  sentences: 

1.  In  a  very  thin  layer  on  the  body,  the  boundary  layer,  the 

velocity  passes  from  the  value  zero  at  the  surface  to  the 
value  which  the  potential  flow  would  have  in  the  neighborhood 
of  the  surface. 

2.  The  pressure  in  the  boundary  layer  is  practically  independent 

of  the  coordinate  normal  to  the  surface  and  equals  the 
pressure  of  the  potential  flow  along  the  surface. 

3.  In  the  boundary  layer  the  only  friction  force  to  be  taken  into 

consideration  is  the  shearing  stress  t  =  p.  — . 

8y 

4.  (without  proof)  The  curvature  of  the  surface  may  be  neglected 

in  the  boundary  layer  as  long  as  the  radius  of  curvature  is 
large  compared  with  the  boundary  layer  thickness  (Boltz, 
Thesis,  (reference  9)). 

5.  All  these  considerations  are  valid  only  as  long  as  no  separation 

of  the  flow  from  the  surface  occurs. 

Without  integration  of  the  boundary  layer  equations  one  can  draw  from 
these  sentences  important  physical  conclusions  as  to  the  flow  pattern: 

In  particular  separation  occutb  if  a  transport  of  boundary  layer  material 
into  the  interior  of  the  fluid  takes  place.  If  a  region  with  pressure 
increase  exists  along  the  body  contour,  the  retarded  fluid  in  the  boundary 
layer  is  in  general,  because  of  its  small  kinetic  energy,  not  able  to 
penetrate  too  far  into  the  region  of  higher  pressure.  It  then  withdraws 
laterally  from  the  region  of  higher  pressure,  separating  from  the  body, 
and  is  deflected  into  the  interior  (fig.  25).  As  the  point  of  separation 
one  defines  the  boundary  between  forward  flow  and  reverse  flow  of  the 
layer  nearest  the  surface,  thus 
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point  of  separation: 


(8.5) 


Determining  when  and  where  separation  occurs  requires  for  each  case  the 
integration  of  the  boundary  layer  equations. 


One  can  readily  understand  that  for  the  velocity  profile  u(y)  at 
the  separation  point  and  for  all  velocity  profiles  in  the  decelerating 

flow  >  0^  an  inflection  point**  must  be  present.  From  equation  (8.3) 

namely,  for  the  surface  y  =  0  there  follows  immediately  the  relation 


(8.6) 


The  curvature  of  the  velocity  profile  at  the  surface  therefore  exchanges 
signs  with  the  pressure  gradient.  Thus  for  flow  with  pressure  decrease 


^accelerated  flow 


M 

dx 


°> 


'afu\ 

■Sy24 


<0  1b  valid  and  therefore 


~  2  '  •  surface 

also  <0  in  the  entire  boundary  layer  (fig.  26).  For  the  region 

of  the  pressure  increase  (decelerating  flow,  ^  >  0^  f  )  ->  0. 

\  dx  /  VdyV 

surface 

However,  since  in  any  case  at  larger  distances  from  the  surface 

!i  <o,  there  must  exist,  for  decelerating  flow,  within  the  boundary 

8y2  .2 

layer  a  point  where  — ^  =  0  (inflection  point)  (fig.  27).  For 

8y 

decelerating  potential  flow  the  boundary  layer  profile  has,  therefore, 
an  inflection  point.  Since  the  separation  profile  with  vanishing  surface 
tangent  must  necessarily  have  an  inflection  point,  it  follows  that 
separation  can  occur  only  when  the  potential  flow  is  decelerating 

(i  >°y 


*  The  velocity  profile  at  the  point  of  separation  therefore  startB 
with  a  vanishing  tangent  (fig.  25).  Velocity  profiles  behind  the  point 
of  separation  have  reverse  flow  in  the  neighborhood  of  the  surface  (fig.  25). 

**  The  presence  of  an  inflection  point  is  significant  for  the 
stability  of  the  velocity  profile  (transition  from  laminar  to  turbulent, 
compare  chapter  XXI  Part  II.) 
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If  separation  Is  present,  the  potential  flow  can  no  longer  envelop 
the  body  closely  everywhere.  Thus  the  pressure  distribution  sometimes 
deviates  considerably  from  that  given  by  the  potential  theory.  In  such 
cases  the  pressure  variation  impressed  on  the  boundary  layer  can  in 
most  cases  be  determined  only  empirically,  because  the  frictionless  outer 
flow  itself  depends  on  the  complicated  phenomena  connected  with  the 
separation. 

Thus  the  boundary  layer  theory  explains  also  the  fact  that  in 
addition  to  the  frictional  drag  a  pure  pressure  drag,  called  "form  drag,” 
appears. 

In  regard  to  later  calculations  the  following  explanation  shall  be 
given:  If  equation  (8.3)  is  differentiated  with  respect  to  y,  there 
results  for  stationary  flow 


2  2 
d  u  8u  8u  8  u  8v  c3u 

dxdy  8y  dx  V  ^ 2  8y  8y 


1  8gp  +  v  8^u 

P  8xdy  +  ^3 


(8.7) 


Due  to  T-n  =  0  and  to  the  boundary 

oy 

from  equation  (8.7)  for  the  surface 


conditions  u  =  v  =  0 
y  =  0  the  relation 


one  obtains 


(8.8) 


which  is  valid  for  all  stationary  boundary  layer  profiles  (pressure 
increase  and  pressure  decrease). 


Frictional  Drag 

As  a  result  of  the  integration  of  the  boundary  layer  equation  one 
obtains  the  velocity  distribution  and  the  separation  point  and  can  there¬ 
from  calculate  the  particularly  interesting  surface  friction  drag  in  the 
following  manner.  The  friction  drag  Wp  results  from  the  integration  of 
the  surface  shearing  stress  over  the  surface  of  the  body.  For  the  plane 
case  one  obtains,  with  the  symbols  according  to  figure  28,  for  the  friction 
drag 


Wp  =  2b 


>2* 


T  COS  Cp,  ds 
o 


'8=0 


(8.9) 
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The  Integration  has  to  be  extended  ov<-r  both  sides  of  the  surface  from 
the  stagnation  point  to  the  separation  point  (s  =  2’),  b  signifies  the 


height  of  the  cylindrical  body.  Because  cos  <p  ds  =  dx  and  tq 
one  obtains  for  the  friction  drag 


W, 


R 


(8.10) 


This  integration  also  has  to  be  extended  along  both  sides  of  the  body. 
For  calculation  of  the  friction  drag  one  needs,  therefore,  the  velocity 
gradient  at  the  surface.  The  latter  can  only  be  obtained  by  integration 
of  the  boundary— layer  differential  equations.  If  the  separation  point 
appears  ahead  of  the  trailing  edge,  the  formula  has  to  be  applied  only 
up  to  the  separation  point  or,  sometimes,  up  to  the  point  of  laminar— to— 
turbulent  transition  which  is  located  further  upstream.  Behind  this 
transition  point  exists  turbulent  surface  friction  drag.  In  order  to 
obtain  the  total  drag,  the  form  drag  has  to  be  added  to  this  friction 
drag;  however,  the  form  drag  cannot  be  obtained  from  the  boundary  layer 
calculation  in  a  simple  manner. 


CHAPTER  IX.  EXACT  SOLUTIONS  OF  THE  BOUNDARY  LAYER 
EQUATIONS  FOR  THE  PLANE  PROBLEM 
a.  The  Flat  Plate  in  Longitudinal  Flow 


One  of  the  simplest  examples  for  the  application  of  the  boundary 
layer  equations  (8.3)  is  the  flow  along  a  flat  plate.  Let  the  plate 
begin  at  x  =  0,  extend  parallel  to  the  x— axis,  and  be  infinitely 
long  (fig.  29).  Let  the  stationary  flow  of  the  free  stream  velocity  UG 
be  treated.  The  calculations  for  it  were  made  by  H.  Blasius  (reference  8) 
in  his  Gottingen  thesis.  In  this  case  the  velocity  of  the  potential  flow 
is  constant,  thus  dp/dx  =  0.  The  boundary  layer  equations  (8.3)  become 
therefore 


u 


du 

dx 


v 


du 

dx 


+ 


du 

d.y 

dv 

dy 


dy‘ 


> 


0 


(o.n 
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with  the  boundary  conditions 


y  =  0:  u  =  v  =  0 


y  =  oo  : 


a 


U 


o 


(9-la) 


Since  the  entire  system  has  no  characteristic  length,  the  assumption 
suggests  itself  that  the  velocity  profiles  at  various  distances  x  are 
affine  to  each  other,  so  that  one  may  write  u/UQ  =  <p  (y/&).  Let  therein 
5  =  6(x)  represent  a  measure  of  the  boundary  layer  thickness,  increasing 
with  the  length  of  run.  One  arrives  at  an  estimate  of  the  boundary  layer 
thickness  in  the  following  manner: 

According  to  the  former  exact  solutions  of  Navier-Stokes  equations 
(chapter  VI),  for  instance  for  the  non— stationary  problem  of  the  surface 
suddenly  set  in  motion. 


with  t  denoting  the  time  since  the  heginning  of  the  motion.  Applied 

to  the  present  stationary  problem  one  may  substitute  for  the  time  t  the 

time  required  by  a  fluid  particle  to  travel  from  the  leading  edge  of  the 

][ 

plate  to  the  point  x.  This  time  t  =  —  and  one  haB,  therefore,  for 

uo 


the  present  case 


5~ 


(9-2) 


Thus  it  is  useful  to  introduce  for  the  distance  from  the  surface  y  the 
new  dimensionless  coordinate  t)  =  y/5  or  according  to  equation  (9.2) 


n  =  y 


(9.3) 


For  further  calculations  one  observes  that 


^  =  _  i  1. 

chc  2  x’ 


The  continuity  equation  is  again  integrated  by  introduction  of  the  stream 
Function  \|f ;  for  th ;  latter  one  assumes 
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\jr  yvx  Ur  f  ( q ) 


(9.J4) 


Accord  ir.gly  f ( tj )  is  a  dim-nsiunl-ss  stream  function;  for  the  velocity 
components  or."  obtains  from  eq  lat  icn  (Q.t) 


6t  1*  oq  ,  x 

u  -  -  —  =  U  i  1  ( q  ) 

oy  dq  dy  u 


(9 • Ir  ) 


v  =  — 


V  =  I’ 

p 


V  Ur 


>  (9.6) 


(T]f  -  f) 


Furthermore  one  obtains 


8u 

8x 


1  Uo  „«.  8u  TI 

57,f'  a7  =  y0' 


f".  =  !^2_  n" 

vx  *  ^ 2  Vx 


(9.7) 


By  insertion  of  these  values  in  equation  (9.1)  there  results 


un  U  *  U 

- —  Tif«f"  +  -2_  (nf'-f)  f"  =v-2.  f»" 

2x  2x  xy 


or,  after  simplifying,  the  following  differential  equation  for  the  stream 
function  f(q) 


f  f"  +  2  f»"  =  0 


19.8) 


Because  of  equation  (9. la)  and  of  equations  (9-5)  and  (9.6)  the  boundary 
conditions  are 


q  =  0:  f  =  f1  =  0;  q  =  f*  =  1 


(9.9) 


For  fh  ■  present  case  therefore,  there  results  from  the  two  partial 
differential  equations  (9.1),  by  the  similarity  transformation  (equation 
(9.)),  an  ordinary  non-linear  differential  equation  of  the  third  order. 
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The  three  boundary  conditions  according  to  equation  (9*9)  are  sufficient 
for  a  unique  determination  of  the  function  f(rj)  from  this  equation 

A  particular  solution  of  the  differential  equation  (9*8)  is  the 
solution 

f  =  t]  +  constant 

This  solution  corresponds  to  the  potential  flow;  we  shall  revert  to  it 
presently. 

The  general  solution  of  the  differential  equation  (9.8)  cannot  be 
given  in  closed  form.  Therefore  one  must  calculate  it  either  by  numerical 
methods  or  by  series  developments.  Blasius  obtained  the  solution  by  a 
power  series  development  near  T)  =  0  and  an  asymptotic  development 
near  tj  =  00 ,  which  are  combined  at  an  appropriate  point.  Since  the 
method  of  calculation  is  characteristic  of  the  solution  of  the  boundary 
layer  differential  equations  it  will  be  described  in  more  detail.  The 
power  series  around  tj  =  0  is  formulated  in  the  form 

\ 

Ao  p  At 

f(n)  =  A  +  A^  +  — ■  T}^  +  -f  V  +  .  .  . 


f"(T)) 


Ah  9  A, 

=  A_  +  A_ti  +  —  tij  +  .  . 

2  3'  2!  '  31 


(9.10) 


f«"(T))  =  A3  +  A4T) 


a  ^5  ?  ^6  7 

+  Ai.TI  +  — T|*-  +  —  T1J  +  •  » 


21  3'- 


Because  of  the  boundary  conditions  for  tj  =  0  one  has  immediately 


A  =  A  =  0 
o  1 


By  insertion  of  the  equations  (9*10)  dnto  the  differential  equation  (9*8) 
one  obtains 


2A  +  tj  2AU  +  (A22  +  2A^  +  ^7  (^A2  A3  +  2A^  +  .  .  .  =  0 


In  order  to  make  the  equation  (9.10)  represent  a  solution  of  the  differ¬ 
ential  equation  it  is  required  that  in  the  last  equation  all  coefficients 
of  the  single  powers  of  rj  vanish.  First,  one  has  immediately  A^  =  A^ 
and  further  _ 


*Thls  also  follows  at  once  from  equations  (8.6)  and  (8.8). 
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*5  =  -  §  A/;  Ag  =  0;  A?  =  0 

*8  -  -  f  *2  H  =  T  V 


Thus  only  the  coefficients  A^y  A^,  Ag.  .  .  .  are  different  from  zero. 

The  coefficient  A2  remains  at  first  undetermined  since  the  third  boundary 

condition  for  T)  =  »  was  not  yet  satisfied.  The  remaining  coefficients 

A  ,  A_,  A  .  .  .  can  all  be  expressed  in  terms  of  A0.  One  therefore 
5  o  11  2 

sets  up,  with  Aq  =  a,  a  series  for  f(q)  which  progresses  by  powers 
of  T)^,  in  the  following  form: 


00 


n=o 


(9.11) 


The  results  for  the  first  coefficients  are: 

CQ  =  lj  C-l  =  1;  C2  =  11; 

C3  =  375;  C4  =  27897;  =  3  817137 

The  asymptotic  development  near  tj  =  »  is  formulated  in  the  form 

f  =  fl  +  f2  +  .  .  •  (9.12) 


where  the  higher  approximations  are  to  be  small  in  comparison  with  the 
lower  approximations,  for  instance  f^  «  f^.  The  first  asymptotic 

approximation  to  correspond  to  the  potential  flow  is,  as  was  mentioned 
above, 

f1  =  ti  —  3  (9.13) 

For  this  approximation  f^"  =  0,  and  one  obtains  therefore  by  replacing 
the  quantity  ff”  in  the  differential  equation  (9*8)  by  f^f2"  the 

following  equation  for  the  second  asymptotic  approximation: 
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(q  -  3)  f”  +  2fJ*  =  0 


ft " 

-f-  -  i  O  -n) 

f£  2 


log  ^2  -  \  01  -  J  l2  +  0 


If  one  substitutes  for  the  integration  constant  C  =  — 3  /k  +  log  7, 
(7  is  thus  a  new  integration  constant)  one  obtains 


f"  =  7  e 
2 


-H-i-sr 


and  after  integration 


(9.14) 


Because  f^  (°°)  =  1  and  f*>  (*>)  =  0,  the  solution  f*  =  f^  +  f^ 

satisfies  the  third  boundary  condition  f*  (*•)  =  1.  Another  integration 
of  equation  (9.1*0  gives  as  second  asymptotic  solution  for  f  =  f^  +  fg 


f  =  T)  -  3  +  7 


-  y  <1  - 

e  dq 


(9.15) 


This  solution  still  contains  two  integration  constants  3  and-  7 
corresponding  to  the  fact  that  only  one  of  the  three  boundary  conditions 
was  satisfied.  The  asymptotic  solution  can,  in  the  same  manner,  be 
improved  still  farther  by  equating  f  =  f-j_  +  fg  +  ty  The  differential 

equation  for  f^  was  solved  by  Blasius;  a  more  detailed  discussion  is 
unnecessary . 


These  two  solutions,  the  power  series  near  q  =  0  according  to 
equation  (9,11)  and  the  asymptotic  solution  near  q  =  <»  according  to 
equation  (9.15);  now  have  to  be  Joined  together  and  the  three  integration 
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constants  a,  0  and  y  have  thereby  to  be  determined.  This  is  effected 
in  the  following  way:  At  a  point  q  =  q  ,  where  both  solutions  are 

serviceable,  f,  f1  and  f"  from  the  power  series  arid  the  asymptotic 
solution  are  brought  to  agreement.  Because  of  the  differential 
equation  (9*8)  the  higher  derivatives  will  then  automatically  agree. 

In  this  manner  one  obtains  three  equations  for  the  three  unknown 
integration  constants.  The  numerical  calculation  gives 


a  =  0.312;  0  =  1.73?  7  =  0.231 

A  table  for  f,  f*,  f"  taken  from  a  treatise  by  L.  Howarth  (reference  18) 
is  given  in  table  2*.  The  velocity  distribution  in  the  boundary  layer 
u/UQ  =  f 1  ( tj )  according  to  equation  (9*5)  is  represented  in  figure  30. 

In  comparison  with  the  stagnation  point  profile  (fig.  16)  it  is  striking 
that  the  velocity  profile  of  the  plate  flow  near  the  surface  has  a  very 
slight  curvature .  At  the  surface  itself  it  has  an  inflection  point,  that 

is,  for  y  =  0:  2-H  =  o. 

3y2 

The  transverse  component  of  the  velocity  in  the  boundary  layer 
given  by  equation  (9.6)  is  plotted  in  figure  31.  It  is  noteworthy  that 
outside  of  the  friction  layer,  for  q-^» 


v  =  v„  =  O.865  Uc 


The  fact  that  on  the  outer  edge  of  the  boundary  layer  the  transverse 
component  v  /  0  is  caused  by  the  deflection  of  the  potential  flow  from 
the  body  due  to  the  boundary  layer  thickness  Increasing  downstream.  For 
very  large  distances  from  the  wall  (far  In  the  potential  flow)  the 
boundary  layer  solution  does  not  go  over  exactly  into  the  undistrubed 
potential  flow  u  =  U0;  v  =  0.  This  has  to  be  tolerated  as  a  (very 
slight)  deficiency  of  the  boundary  layer  solution. 

For  the  present  case  a  separation  of  the  boundary  layer  does  not 
exist  since  the  pressure  gradient  equals  zero. 


Friction  Brag 

From  the  solution  given  above  the  surface  friction  drag  of  the 
plate  in  longitudinal  flow  is  to  be  calculated.  According  to  equation  (8.9) 
the  friction  drag  for  one  side  of  the  plate  is 


*  See  appendix* chapter  XII. 


58 


NACA  TM  No.  1217 


NACA  TM  No.  1217 


59 


Displacement  Thickness  of  the  Boundary  Layer 

By  the  development  of  the  boundary  layer  on  the  plate,  which 
increases  downstream  with  \Jx ,  the  potential  flow  is  deflected  outward 
from  the  surface  by  an  amount  5*,  which  is  called  the  displacement 
thickness  of  the  boundary  layer.  It  can  be  easily  calculated  from  the 
velocity  distribution  in  the  boundary  layer,  as  follows:  Let  y^  denote 
a  point  outside  the  boundary  layer;  then  according  to  the  definition 
for  8* 


lyi 

udy  =  U0('y1  -8*^ 

J  y=o 


or 


(9.20) 


According  to  equation  (9.5) 

•*y=o  Jt)=o 


Since  the  point  t|  =  lies  outside  of  the  boundary  layer,  one  can  put 

for  f(q)  the  first  approximation  of  the  asymptotic  solution  according 
to  equation  (9.13)*  thus 


f(T)l)  =  “  P  =  TJi  -  1.73 


Thus  one  finds  for  the  displacement  thickness  of  the  boundary  layer 


6o 
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5*  =  1.7 i 

^ o 


(9.21) 


The  distance  from  the  surface  y  =  5*-  is  also  shown  in  figure  30.  Thus 
the  stream  lines  of  the  potential  flow  are,  because  of  the  friction 
effect,  deflected  outward  by  this  amount. 

The  actual  boundary  layer  thickness  5  cannot  be  given  accurately 
since  the  friction  effect  in  the  boundary  layer  thickness  decreases 
asymptotically  toward  the  outside.  The  component  of  velocity  parallel 
to  the  surface  u  is  asymptotically  converted  into  the  velocity  UQ 
of  the  potential  flow  (the  function  f* •  ( tj )  asymptotically  approaches  the 
value  1).  If  one  wants  to  define  the  boundary  layer  thickness  as  the 
point  where  the  velocity  u  =  O.99  UQ  (full  value),  one  obtains  for  it 
according  to  table  2,  T)  =  5*0*  Therewith  one  has  for  the  thus  defined 
boundary  layer  thickness 


(9.21a) 


The  thus  defined  boundary  layer  thickness  equals  about  three  times  the 
displacement  thickness  of  the  boundary  layer. 


Let  here  also  be  introduced  the  value  for  the  momentum  thickness  a, 
needed  later.  This  latter  is  a  measure  for  the  momentum  loss  due  to 
friction  in  the  boundary  layer  and  is,  as  indicated  before  in  equation  (6.32), 
defined  by  the  equation 


a  = 


1** 

Jy=0 


The  calculation  results,  because  of  equation  (9*5)  >  in 


a 


1  00 

fd-fUd,^-  0.664  ^ 

Ut)=o 


(9.21b) 


Finally  the  form  parameter  becomes  therewith 


NACA  TM  No.  1217 


61 


—  =  2.605 
d 


(9.21c) 


Experimental  investigations  of  the  laminar  "boundary  layer  on  the  flat 
plate  were  performed  by  B.  G.  van  der  Hegge  Zynen  (reference  19)  and 
M.  Hansen  (reference  20).  In  all  essential  points  the  theoretical  results 
were  well  confirmed.  The  measurements  showed  that  the  laminar  "boundary 
layer  exists  to  about  the  Reynolds  number  (U  x/v)  =  3.5  to  5  x  10^, 

0  crit 

if  x  denotes  the  length  of  run  of  the  "boundary  layer.  For  larger 
Reynolds  numbers  transition  to  the  turbulent  state  of  flow  takes  place. 


Fifth  lecture  on  January  5*  19^2. 
b.  The  Boundary  Layer  on  the  Cylinder  (symmetrical  case) 

The  integration  method  of  Blasius  given  in  the  previous  section 
was  used  by  Hiemenz  (thesis  Gottingen  1911)  for  calculating  the  boundary 
layer  on  the  circular  cylinder.  The  same  method  was  later  further  extended 
by  Howarth  (reference  15)  to  the  general  case  of  a  cylindrical  body  of 
arbitrary  cross  section.  This  method  will  be  briefly  presented  for  the 
symmetrical  case.  One  considers  (fig.  32)  a  cylindrical  body  with 
symmetrical  cross  section  in  a  flow  approaching  in  the  direction  of  the 
symmetry  axis  with  the  velocity  U0.  Let  x  be  the  arc  length  along 

the  contour,  measured  from  the  front  stagnation  point,  y  the  vertical 
distance  from  the  surface.  Let  the  potential  flow  U(x)  be  given  by 
its  power  series  development  in  x.  At  the  stagnation  point  (x  =  0), 

U(x)  =  0,  and  for  the  symmetrical  case  only  the  odd  terms  of  the  power 
series  are  different  from  zero.  Therefore: 


U(x) 


3  5 

ux  +  ux  +ux  + 

1  3  5 


(9.22) 


The  coefficients  u-^,  u^,  .  .  .  depend  solely  on  the  shape  of  the 

body  and  are  therefore  quantities  known  from  the  potential  flow. 


The  stationary  boundary  layer  equations  according  to  equation  (8.3) 
are  also  valid  for  this  case  with  a  curved  surface  and  therefore  read: 

\ 


8u  6u  TT  dU  ,,  c^u 

u  —  +  v  —  =  U  —  +  V  - 

dx  oy  dx  ^y2 

^  ^  =  0 

dx  dy 


s 


(9.23) 
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From  equation  (9.22)  one  obtains  for  the  pressure  term: 


U  ffl  = 
dx 


u. 


u^x 


Uu^x^  + 


(9.24) 


The  continuity  equation  is  again  integrated  by  the  stream  function: 


u 


dx 


(9-25) 


It  is  now  necessary  to  find  a  suitable  formulation  for  the  velocity 
distribution  u(x,  y),  v(x,  y)  and  therewith  for  the  stream 
function  +(xj  y).  In  analogy  to  equation  (9.22)  a  power  series  in  x 
suggests  itself  for  u(x,  y)  as  well,  with  coefficients,  however,  which 
are  dependent  on  y.  It  is  important  to  find  a  form  where  the  coeffi¬ 
cients  (or  functions)  dependent  on  y  have  a  universal  character,  that 
is,  need  not  be  calculated  anew  for  each  shape  of  body,  but  may  be  calcu¬ 
lated  once  for  all.  Howarth  (reference  15)  succeeded  in  finding  such  a 
formulation. 


For  the  distance  from  the  surface  one  introduces  the  dimensionless 
variable: 


*(9.26) 


The  expression  (9.2U)  for  the  pressure  term  suggests  that  the  following 
equation  be  selected  for  i < 


+  =l/~<|ulxfl(T1)  +  +  6x5 


“585(1)  +  h^(T))  j 


+  •  •  • 


1 


(9- 


This  yields: 


(*  =  differentiation  with  respect  to  rj) 

2 


u  = 


u^xfj  +  4u^x^f’  +  6x; 


u 


u5g5  +  ^  h5 


(9.28) 


♦One  obtains  this  equation  from  that  of  Blasius  according  to 
equation  (9.})  by  substituting  for  TJ0  the  first  term  of  the  series 

equation  (9.22). 
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du 

8x 


u  f*  +  12u  x2f»  +  30x^ 
11  3  3 


u 


3 


u  g*  +  -2-  h* 


5  5  u-|  5 


(9-29) 


1 - 

1 _ 

S' 

J 

2 

-N 

^l> 

II 

3  |  >» 

so 

u  xf"  +  4u  x^f"  +  Gyp 
11  3  3 

U3 

u  gM  +  —  h" 

5  5  u±  5 

+  «  •  •  > 

1 

V 

J 

S.P.  U, 

O  U  _  _1 

^y2  V 


r 


<u-jXf|"  +  4u^x3f«w  +  6x5 


u. 


cug*"  +  - h»" 

5  5  ux  5 


+  .  .  .  >  (9.31) 


J 


' 

u  2 

T 

<  Ui^i  +  12u^x2f^  +  30x^ 

u5®5  +  —  h5 

+  .  .  A 

> 

(9.32) 


After  insertion  of  the  expressions  (9.24)  and  (9.28)  to  (9.32)  into  the 
first  equation  (9.23)  one  obtains  by  comparison  of  the  coefficients  a 
system  of  ordinary  differential  equations  for  the  unknown  functions  f-^ 
f-j,  g^,  hs,  .  ,  which  appears  as  follows: 


terms  with 

u  x 
1 


4u^u^x3 


6u^uc_x' 


6u  2x^ 
3 


gives  the  differential 


f2  -  f 

'  f”  =  1  +  f»” 

1 

1  1 

1 

4f»f! 

-  3f"f.  -  f.f"  = 

1  3 

1  3 

1  3 

6fi«5 

-  -  fig;  = 

6f*h» 

-  5f"h 

-  f  h"  = 

1  5 

1  5 

1  5 

equation 

1  +  f»" 

1  +  8^" 

—  +  h,w  —  8(f»2  -  f  f") 

2  5  3  3  3 


>  (9-33) 


|  J 

Formulation  of  the  flow  function  according  to  equation  (9.27)  nas 
thus  accomplished  the  elimination  of  the  coefficients  depending  on  body 


64 


NACA  TM  No.  1217 


shape  (u-^,  u-^,  .  .  .  )  from  the  differential  equations'  fd*r  the  functions 
f-|_,  fy  .  •  .  which  thus  now  have  a  universal  character. 


The  boundary  conditions  for  the  functions 

from 


follow 


y  =  0:  u  =  v  =  0 
y  =  ••:  u  =  U 


by  comparison  of  equation  (9.28)  with  equation  (9.22);  they  are 


q  =  0:  f  = 


T)  = 


f*  =0;  f  =  f»  =  0; 

13  3 


g5  g5  ° 


h  =  h*  =  0; 

6  5 


f* 

1 


(9-3*0 


1  •  f*  I  —  i  *  nr  I  —  V)  f  —  A* 

’  L  i.  >  si  n  -  u,  .  .  . 

3  4  5  9  5 


The  differential  equations  (9.33)  are  all  of  the  third  order,  and 
equation  (9.34)  gives  three  boundary  conditions  for  each.  The  differen¬ 
tial  equation  for  f,(q)  is  non-linear  and  is  identical  with  the  differ¬ 
ential  equation  (6.28)  obtained  in  chapter  VT  for  the  stagnation  point 
flow:  ^  =  cp;  q  =  |,  as  follows  by  comparison  of  equation  (6.26a) 

with  equation  (6.26).  All  the  remaining  differential  equations  in 
equation  (9 . 33)  are  linear,  with  the  coefficients  determined  by  the 
functions  of  the  preceding  approximations. 

The  solutions  of  the  differential  equations  (9*33)  are  best  obtained 
by  numerical  integration.  The  functions  f1  and  f^  were  already 

calculated  by  Hiemenz  (reference  10).  Howarth  (reference  15)  improved 
the  tables  for  f^  and  recently  Nils  Frossling  (reference  16)  calculated 
g  and  h  as  well.  The  f!  which  is  identical  with  cp*  according  to 

equation  (6.28)  was  already  given  in  figure  16.  The  function  fl,  can  bo 
seen  from  figure  33,  the  functions  g^  and  h*.  from  figure  34.  The 
numerical  values  are  compiled  in  table  3. 

Concerning  the  applicability  of  this  calculation  method  it  mist  b' 
mentioned  that  for  slender  body  shapes  the  series  for  U(x)  and  u(x,y) 
converge  poorly.  The  reason  is,  that  for  these  body  shapes  U(x)  has 
a  very  steep  ascent  In  the  neighborhood  of  the  stagnation  point  (fig.  i r- )  > 
while  showing  a  rather  flat  curve  further  on.  Gucn  a  function  cannot.  I)', 
developed  readily  into  a  Taylor  series.  For  such  body  shapes  many  more  of 
the  functions  of  the  differential-equation  system  (9*33)  would  be  required 
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than  have  been  calculated  bo  far.  For  blunt  bodies,  as  for  Instance  the 
circular  cylinder,  the  convergence  is  considerably  better  so  that  one 
proceeds  rather  far  with  this  calculation  although  not  always  up  to  the 
point  of  separation. 

Howarth  (reference  lb)  also  performed  the  corresponding  calculation 
for  the  unsyrametri cal  case  where  the  even  coefficients  also  appear  in 
the  power  series  for  U(x).  This  is  the  case  for  a  symmetrical  body 
at  an  angle  of  attack  and  quite  generally  for  any  un symmetrical  body. 

Frossling  (reference  l6)  made  the  application  to  the  rotationally 
symmetrical  case. 


Circular  Cylinder 

The  boundary  layer  on  the  circular  cylinder  will  be  treated  as  an 
example  for  the  application  of  this  method.  Whereas  Hiemenz  (reference  10) 
took  a  pressure  distribution  measured  by  him  as  the  basis  for  this  case, 
we  shall  here  calculate  with  the  potential— theoretical  pressure  distribution. 
The  velocity  distribution  of  the  potential  flow  reads,  with  the  symbols 
according  to  figure  3&> 


U(x)  =  2U0  sin  <p  =  2U0  sin  £ 

E 


(9.35) 


The  power  series  development  gives: 
U(x)  =  2  U 


r 

X  _ 

X.  { 

+  JL 

R 

3'.  V 

R  /  5’. 

[s)  i 

(9.36) 


In  comparison  with  equation  (9-22): 


U  U  U 

_  o  2°  2  0 

u,  =  2  — ;  u0  - - - - :  uR  =  — - : 

1  *  3  5  5!r5' 


(9.36a) 


Therewith  follows  from  equation  (9.26): 


2U£  =  y  .U£ 
V  V 


(9-37) 


It  follows  that  for  the  velocity  distribution  from  equation  (9.28) 
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1  u(x,y)  =  x  f, 

2  U0  R  1 


4-  •  •  • 


(9. 38) 


One  further  calculates  the  position  of  the  separation  point  xA  which 

is,  according  to  equation  (8-5)  determined  by  (—1  =0.  Therewith 

\oy/ 

y=o 

results  from  equation  (Q.36): 


f ”  ( 0 )  ^  -‘Jl  f" 

I  B  3!  3 


+  •  • 


.  =  0 


(9.38a) 


With  the  numerical  values 


f|(0) 


1.23 264; 


f”(0j  =  0.7246 


one  finds: 


^  =  1.60;  cp 
R  A 


(9.39) 


Hiemenz  (reference  10)  based  his  calculation  on  his  experimental 
pressure  distribution;  he  calculated  the  separation  point  to  be  at 

cp  =  82°,  whereas  his  measurements  gave  cp  =  81°.  This  result  is 
A  A 

considerably  different  from  that  obtained  with  the  potential— theoretical 
pressure  distribution.  The  reason  is  that  for  a  body  as  blunt  as  the 
circular  cylinder  the  experimental  and  the  potential— theoretical  pressure 
distribution  in  the  neighborhood  of  the  separation  point  differ  greatly. 

The  method  described  here  of  calculating  the  boundary  layer  by  a 
power— 3eri.es  development  starting  from  the  stagnation  point  has  found 
but  little  acceptance  because  of  the  extensive  calculation  required. 

It  is,  however,  indispensable  for  fundamental  considerations,  since  there 
exist  no  other  exact  solutions  of  the  differential  equations  of  the 
boundary  layer  for  the  flow  about  a  body. 

Thus  approximation  methods  came  into  use  for  the  practical  per¬ 
formance  of  boundary  layer  calculations;  they  will  be  discussed  in  the 
following  chapter.  It  is  true  that  their  accuracy  is  sometimes  consider¬ 
ably  lower  than  that  of  the  previously  discussed  exact  solutions. 


c.  Wake  behind  the  Flat  Plate  in  Longitudinal  Flow 


The  application  of  the  boundary  layer  equations  is  not  absolutely 
limited  to  the  presence  of  solid  walls.  They  may  also  be  applied  if 
tier"  is  present  within  the  flow  a  layer  in  which  the  friction  effect  is 


•This  result  varies  somewhat  if  in  the  series  development 
equation  (0.38a),  one  tak<-*a  further  terms  into  consideration.  However 
for  this  purpose  one  would  have  to  calculate  at  least  up  to  the  term 


,rj 
x  1 
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predominant.  This  Is  the  case  for  instance  when  within  the  flow  two 
layers  of  different  velocities  adjoin,  as  for  instance  in  the  wake  region 
behind  a  body  or  at  the  outflow  from  an  opening.  In  this  and  the  following 
chapter  we  shall  treat  two  examples  of  such  flows  which  we  shall  later 
encounter  again  in  the  discussion  of  turbulent  flows. 

The  wake  flow  behind  the  flat  plate  in  longitudinal  flow  is  chosen 
as  the  first  example  (fig.  37).  At  the  trailing  edge  of  the  plate  the 
two  boundary  layer  profiles  grow  together  and  form  a  "wake  profile"  the 
width  of  which  increases  with  distance  while  the  velocity  decrement  at 
its  center  decreases.  The  size  of  the  "wake"  is  directly  connected  with 
the  drag  of  the  body.  Otherwise,  however,  the  form  of  the  velocity 
distribution  in  the  wake  at  a  large  distance  from  the  body  is  not  dependent 
on  the  shape  of  the  body,  whereas  the  velocity  distribution  very  close 
behind  the  body  naturally  depends  on  the  boundary  layer  of  the  body  and 
on  any  existing  separation. 

From  the  velocity  distribution  measured  in  the  wake  one  may  cal¬ 
culate  the  drag  by  means  of  the  momentum  theorem  in  the  following  manner. 
The  momentum  theorem  states:  The  variation  of  the  momentum  with  time 
(  =  momentum  flow  through  a  control  area  fixed  in  space)  equals  the 
sum  of  the  resultant  forces.  By  resultant  forces  one  has  to  understand: 

1.  Pressure  forces  on  the  control  area,  2.  Extraneous  forces,  which  are 
transferred  from  solid  bodies  to  the  flowing  fluid,  for  instance  the 
shearing  stress  at  the  surface  which  gives  the  friction  drag.  For  the 
present  case  the  control  area  AAjBB^  is  placed  as  Indicated  in  figure  38. 

Let  the  boundary  A^B1  which  is  parallel  to  the  plate  be  so  far  distant 

from  it  that  it  lies  everywhere  in  the  undisturbed  velocity  U0.  Further¬ 
more,  the  rear  cross  section  BB^  1b  to  lie  so  far  behind  the  plate  that 

the  static  pressure  there  has  the  same  undisturbed  value  as  in  front  of 
the  plate.  Then  the  pressure  is  constant  on  the  entire  control  area,  so 
that  the  pressure  forces  make  no  contribution.  In  calculating  the 
momentum  flow  through  this  control  area  one  has  to  consider  that,  due  to 
continuity,  fluid  must  flow  out  through  the  boundary  A-jB^,  namely  the 

difference  between  the  larger  quantity  flowing  through  the  cross 
section  AAj  and  the  smaller  quantity  flowing  through  the  cross 

section  BB^.  The  cross  section  AB  does  not  make  a  contribution  to  the 
x— momentum,  since  for  reasons  of  symmetry  the  transverse  velocity  on  it 
equals  zero.  The  momentum  balance  is  given  in  the  table  below,  with 
entering  momentums  counted  a3  positive,  outgoing  ones  as  negative. 
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Cross  section 

Mass 

m^/secj 

Momentum  flow  in  x— direction 
=  volume  X  density x velocity 

A  B 

0 

0 

1*i 

A  A± 

j 

uG  dy 

P  b 

UD2  dy 

«ic 

5 

U( 

3 

1h 

1*1 

B  B1 

-  b 

u  dy 

-  p  b 

u2  dy 

1 

) 

Jc 

3 

1*1 

A1B1 

-  \ 

j  (Uo  ~ 

Jo 

-pbU, 

5  (UD  -  u)  dy 

Jo 

z- 

Mass 

Momentum  Flow 

Control  area 

=  0 

=  W 

The  total  momentum  loss  of  the  flow  for  the  present  case  equals  the 
drag  W  of  one  side  of  the  plate.  Thus  one  obtains 


W  =  b 


dy 


(9-40) 


The  integration  therein  may  be  extended  from  y  =  0  to  instead  of 

from  y  =  0  to  h,  since  for  y>  h  the  integrand  in  equation  (9.40) 
vanishes.  For  the  drag  of  the  plate  wetted  on  both  sides  one  obtains 
therefore : 
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+  eo 


2W  =  bp 


u(U0-u)  dy 


(9-41) 


In  equations  (9.40)  and  (9.4l)  the  integrals  are  to  be  extended  over  the 
wake  as  indicated  at  a  distance  behind  the  plate  where  the  static  pressure 
has  its  undisturbed  value.  However,  one  may  naturally  apply  equations 
(9.40)  and  (9.41)  also  to  the  boundary  layer  on  the  plate  at  a  certain 
point  x;  then  they  give  the  drag  of  the  part  of  the  plate  from  the 
leading  edge  to  this  point.  The  definite  integrals  in  equations  (9.40) 
and  (9.41)  represent  physically  the  momentum  loss  due  to  the  friction 
effect.  As  mentioned  before,  it  is  customary  to  introduce  for  this 
integral  also  the  momentum  loss  thickness  d  by  the  following  equation 
(compare  equation  (6.32)). 


r 

!  u  (Uo  -  u)  dy 


dy=D 


(9.42) 


Therewith  the  formula  for  the  drag  may  also  be  written,  by  comparison 
with  equation  (9.40): 


2 

W  =  b  p  UQ  d 


(9-^3) 


The  velocity  distribution  in  the  wake,  particularly  at  large 
distance  x  behind  the  plate  in  longitudinal  flow  (fig.  37)  is  to  be 
calculated  next.  This  calculation  must  be  performed  in  two  steps: 

1.  By  a  development  "from  the  front",  that  is,  by  a  calculation  which 
follows  the  further  development  of  the  Blasius— boundary  layer  profile 
present  at  the  trailing  edge  of  the  plate.  2.  By  a  development  "from 
the  rear"  that  is,  by  an  asymptotic  calculation  for  the  wake,  under  the 
assumption  that  the  difference  velocity 


ux(x,  y)  =  U0  -  u(x,  y) 


(9.44) 


is  s:nall. 

The  first  calculation  was  performed  by  S.  Goldstein  (reference  21) 
and  is  very  troublesome;  the  second  was  indicated  by  Tollmien  (reference  3) 
and  yields  rather  simple  laws,  in  particular  also  an  exact  solution  of  the 
«  differential  boundary  layer  equation.  Since  such  exact  solutions  are 
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comparatively  rare  and  since,  moreover,  the  asymptotic  law  for  the  wake 
applies  not  only  for  the  flat  plate  in  longitudinal  flow  but  for  any 
arbitrary  body  shape,  this  asymptotic  solution  will  be  treated  here 
somewhat  more  thoroughly.  The  wake  velocity  u^(x,  y)  introduced  in 

equation  (9.1*4)  is  assumed  to  be  so  small  in  comparison  with  the  free 
stream  velocity  UQ  that  the  second— power  terms  {u^/Uo)2  are  negligible 

relative  to  1.  Moreover,  the  pressure  term  dp/dx  in  the  boundary  layer 
equation  is  to  be  set  equal  to  zero  for  the  first  asymptotic  approximation. 
Therewith  the  differential  boundary  layer  equation  (8.3)  becomes,  for  the 
present  case: 


With  the  boundary  conditions: 


cku 

y  =  0:  —  =  0 

ay 

y  =  »:  u±  =  0  (u  =  UQ) 


(9.45a) 


For  the  solution  of  the  differential  equation  (9.45)  one  introduces  as 
before  in  the  plate  flow  according  to  Blasius*  equation  (9.3)  the  new 
variable 


(9.46) 


Further,  one  uses  for  u  the  equation: 


u-l  =  UD  C 


(9.47) 


The  distance  x  from  the  trailing  edge  of  the  plate  is  thus  made 
dimensionless  by  dividing  by  the  plate  length.  The  power  -  -  for  x 

1b  given  bv  the  fact  that  the  momentum  integral  which,  according  to 
equation  (9. hi)  gives  the  plate  drag  must  be  independent  of  x.  With 
the  second-power  terms  neglected  the  drag  of  the  plate  wetted  on  both 
sides  is,  according  to  equation  (9.41): 
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I  +  OO 


2W  =  b  p  U0  j  u.^  dy 

>y=-, 


(9-48) 


One  finds  further: 


)+0» 


\  +  00 


=  Uo°(f)  2\|f|  8(T,)dn- 


n  +« 


uo°K 


U- , 


»T)=—  M 


g(i))dTl  (9.^9) 


and  therewith: 


2W  =  b  p  U0  Cy -jk  |  «(tj)  dT, 


(9.50) 


The  calculation  of  the  single  terms  in  equation  (9.45)  gives: 

\ 

^1  _  TT  C  /x  i~  2  1/  ,x 

S - uo  2  (jJ-  j  (e  +  ne') 


du 

dy 


1=uoc(!rJ&i 


*2 
O  Ui 


8y' 


-  U  C  (y 
o  V  Z 


_  1 

2  1 


7  8 


>  (9.5D 


/ 


By  insertion  into  equation  (9.45)  one  obtains  after  division  by 

2  -  \ 

C  U0  (x/l)  x  the  following  differential  equation  for  the  velocity 
distribution  g(t)): 


g"  +  |  18'  +  |  g  =  0 


(9.52) 
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with  the  boundary  conditions: 


T)  =  0:  g’  =  0;  t]  =  oo  :  g  =  0 


A  single  integration  gives: 


S*  +  ^rjg  =  K-  0 


where,  because  of  the  boundary  condition  at  q  =  0,  the  integration 
constant  K  must  be  zero.  Repeated  Integration  gives  the  solution: 


g 


(9-53) 


where  a  multiplicative  integration  constant  may  be  put  equal  to  one 
without  limiting  the  generality  since  the  velocity  distribution  u^ 

still  contains,  according  to  equation  (9. 47),  the  multiplicative  free 
constant  C.  This  constant  C  is  determined  from  the  consideration 
that  the  plate  drag  obtained  from  the  momentum  loss  (equation  (9. 50)) 
must  be  the  same  as  the  frictional  drag  of  the  plate.  First, 


n  +» 


n+‘ 


8(tj)  dr)  = 


U  — 


-  n£ 


e  dq  =  2  \J~n 


and  therewith  from,  equation  (9-50): 


2W  =  p  b  UQ 


2 


C 


On  the  other  hand,  the  friction  drag  of  the  plate  wetted  on  both  sides 
was  according  to  the  solution  of  Blasius  (equation  (9.18)) 


2W  =  I.328  p  UQ 


2  JvT 

H 


Therefrom  follows  2C  \Jn  =  I.328  and 
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0.664 

\f  * 


(9.54) 


Thus  the  final  solution  for  the  wake  velocity  for  the  flat  plate  in 
longitudinal  flow  becomes: 


_  1  1  7  Uo 

_1  =  0 » 664  /  x  \  ?  4  xv 

uc  \pr  [i) 


The  velocity  distribution  of  this  asymptotic  law  is  represented  in 
figure  39-  It  is  noteworthy  that  the  function  for  the  velocity  distribution 
is  identical  with  Gauss’  error  function.  In  keeping  with  the  hypothesis 
the  law  according  to  equation  (9-55)  is  valid  only  for  large  distances 
behind  the  plate,  according  to  the  calculations  by  Tollmien  (reference  3) 
for  x  -  32. 

The  development  of  the  wake  from  the  front,  performed  by  Goldstein 
is  valid  only  for  comparatively  small  x/Z.  However,  for  intermediate 
x/z  both  solutions  can  be  Joined  to  some  extent,  so  that  one  then  obtains 
the  velocity  distribution  in  the  entire  wake.  Such  a  figure  is  given  by 
Tollmien  (reference  3). 


Sixth  Lecture  on  January  12,  1942. 
d.  The  Plane  Jet 

A  further  example  of  a  flow  without  bounding  wall  to  which  the 
boundary  layer  theory  is  applicable  is  the  outflow  of  a  Jet  from  a  hole. 
The  problem  to  be  treated  is  the  plane  stationary  one  where  the  Jet  goes 
out  from  a  long  narrow  slot  and  mixes  with  the  surrounding  fluid  at  rest. 
This  is  one  of  the  rare  caseB  where  the  differential  boundary  layer 
equations  may  be  integrated  exactly.  The  calculations  were  performed  by 
H.  Schlichting  (reference  22)  and  W.  Bickley  (reference  30)  and  will  be 
discussed  a  little  more  thoroughly. 

Due  to  the  friction  effect  the  Jet  entrains  a  part  of  the  fluid  at 
rest  and  sweeps  it  along.  There  results  a  stream-line  pattern  like  the 
one  drawn  in  figure  40.  Let  the  x— direction  coincide  with  the  Jet  axis 
and  the  origin  lie  in  the  slot.  It  then  immediately  becomes  clear  that 
the  width  of  the  Jet  increases  with  the  distance  x  and  the  mid— velocity 
decreases  with  the  distance  x.  For  the  calculation  the  slot  is  assumed 
to  be  infinitely  narrow.  In  order  to  make  the  volume  of  flow,  together 
with  its  momentum,  finite,  the  velocity  in  the  slot  is  then  infinitely 
large.  Again,  as  in  the  previous  example  concerning  the  wake  flow,  the 
pressure  term  dp/dx  may  be  neglected  since  the  pressure  varies  only 
very  little  in  the  x-direction.  The  3maliness  of  the  pressure  term  can 


7U 


NACA  TM  No.  1217 


also  be  shown  subsequently  from  the  finished  solution.  Thus  the  differ¬ 
ential  boundary  layer  equations  for  the  present  case  read,  according  to 
equation  (8.3): 


bu  8u  S^u 

u  —  +  v  —  =  V  - 

ck  by  .2 

dy 


8u  ,  bv 

SF  +  ¥ 


o 


r  (9.56) 


with  the  boundary  conditions: 


J 


y  =  0: 


y  =  oo;  u  =  0 


> 


(9.56a) 


Since  the  pressure  is  constant,  the  entire  momentum  flow  in  the  x— 
direction  for  control  area  fixed  in  space  (compare  figure  VO)  must  be 
independent  of  the  distance  from  the  hole  x.  If  one  chooses  the  lateral 
boundaries  of  the  control  area  at  so  large  a  distance  from  the  axis  that 
there  u  =  0,  then 


J  =  p 


n  +00 


u2dy 


=  Constant 


(j— 00 


(9.57) 


It  is  to  be  noted  for  the  integration  of  the  equation  of  motion 
equation  (9.96)  that  for  this  problem,  as  before  for  the  plate  in  longi¬ 
tudinal  flow,  no  special  length— dimension  exists.  Thus  affinity  of  the 
velocity  profiles  u(x,  y)  is  suggested,  that  is:  with  b  signifying 
a  suitable  width  of  the  Jet,  the  velocity  profiles  are  only  functions 
of  y/b*  Accordingly  one  uses  the  following  expression  for  the  stream 
function  \|r: 
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The  two  —  at  first  unknown  —  exponents  p  and  q  are  determined  from 
the  conditions  that 

1.  the  momentum  flow  for  the  x-direction  ' s  independent  of  x 

according  to  equation  (9*57)*  and  that 

2.  the  acceleration  terms,  for  instance  u  — ,  and  the  inertia 

dx 

term  in  equation  (9.56)  are  of  the  same  order  of  magnitude 
and  hence  must  "be  of  the  same  degree  in  x. 

This  yields  the  two  equations 


2p  —  q  -  0  I 

2p-2q-l=p-3q 

) 

It  follows  that 


P  =  1/3;  q  =  2/3 


(9-59) 


(9.60) 


Therewith  the  final  equations,  after  addition  of  suitable  constant  factors, 
read  as  follows: 


1/2  2/3 
3v  x 


(9-61) 


^  =  v1/2  x1//3  f(n) 


(9.62) 


Therefrom  one  obtains,  with 


h  =  _  ^  a 

ax  3  X' 


ih  = 
dy 


-2/3 

3v 


the  following  expressions  for  the  velocity  components  and  their  derivatives 
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u 


v 


1 


f’(T)) 


-  I  Vl/2  X-2/3 

3 


(f-2rjf»)  I 


(9.63) 


y 


|U  =  x  1  1  v  1/2  f" 

oy  9 


-\2  -5/3 

v  £Ji  =  x  J_  ft" 
V  .  2  27 
ay 


By  substitution  into  the  differential  equation  (9.56)  there  results, 

1  -5/3 

after  cancelling  the  factor  x  ,  the  following  differential  equation 
for  the  flow  function  f(q): 


f ' 


2 


+ 


ff"  +  f,n  =  0 


(9.64) 


with  the  boundary  conditions: 


tj  =  0:  f  =  0;  f"  =  0 
t)  =  «;  f*  =  0 


(9.64a) 


As  for  Blasius*  plate  flow  here  also  an  ordinary  differential  equation  (9.64) 
was  obtained  from  the  two  partial  differential  equations  (9*56)  by  means  of 
the  similarity  transformation  equation  (9.61).  Here  also,  as  in  most 
boundary— layer  problems,  the  differential  equation  is  non— linear  and  of 
the  third  order. 

The  integration  of  this  differential  equation  (9.64)  is  attained  in 
a  surprisingly  simple  manner.  First,  one  obtains  by  a  single  integration 


ff’  +  f"  =  Constant  =  0 


(0.6M 
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The  constant  of  integration  therein  is  zero  because  of  the  boundary 
condition  f"(o)  =  0.  The  second  order  differential  equation  now  obtained 
could  be  integrated  once  more  if  a  factor  2  were  present  in  the  first 
term.  This  factor  can  be  obtained  by  performing  the  following  further 
similarity  transformation: 

One  puts:  4  =  a  t]  (9.66) 

f  =  2a  F(4)  (9.67) 

a  is  a  free  constant  which  is  determined  later.  With  the  equations  (9.66) 
and  (9.67)  one  obtains  from  equation  (9.65) ,  the  prime  now  signifying 
differentiation  with  respect  to  4, 


F"  +  2FF*  =  0  (9.68) 


with  the  boundary  conditions: 


4=0:  F  =  0;  4  =  F»  =  0 


(9»68a) 


This  differential  equation  can  now  be  integrated  again  and  yields 

F’  +  F2  =  K 


The  constant  of  integration  K  is  obtained  from  the  boundary  conditions, 
equation  (9.68a),  as  K  =  1,  if  one  stipulates  F*(o)  =  1,  which  is 
possible  without  limiting  the  generality  because  a  is  still  present  as 
a  free  constant  in  F.  One  now  has  for  F  the  first  order  non— linear 
differential  equation 


F«  +  F 


(9.69) 


which  is  a  Riccati  differential  equation.  The  integration  yields 


t 


arc  tanh  F 


and  therewith  for  the  inverse  function 
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F  =  tanh  I 


-2? 

1  —  e 

-2t 

1  +  e 


(9.70) 


Furthermore,  77-  =  1  —  tan2  If  one  inserts  the  solution  found  nto 
d! 

equation  (9.63),  one  obtains  for  the  velocity  distribution 


u 


2  2 

=  —  or 


_  1 

x  3 (1  — tanh2! ) 


(9.71) 


The  velocity  distribution  over  the  width  of  the  Jet  calculated  from  this 
equation  is  represented  in  figure  41.  The  free  constant  a  remains  to 
be  determined.  This  can  be  done  from  the  condition  (equation  (9. 57)) 
according  to  which  the  momentum  flow  in  the  x-direction  is  independent 
of  x.  From  equation  (9.71)  and  (9*57)  one  obtains 


J  =  |  p  a,3 v1/2  (l  -taniPl)  d|  =  ^  p  c^v1/2 


(9-72) 


Let  the  momentum  J  for  the  Jet  be  a  prescribed  constant  which  is,  for 
instance,  directly  related  to  the  excess  pressure  under  which  the  Jet 
flows  from  the  slot.  Then  a  becomes,  according  to  equation  (9.72), 


a  =  0.8255 


(9.73) 


and  therewith  the  velocity  distribution 

\l/3 


U  =  0.4543 


2 

vx 


(l — tanh2 I ) 


1/3 


v  =  0.5503  i  2  j  j  2^1- tanh2|  tanh 


I  =  0.2752  I  -±Y% 

P  V  /  X 


(9.74) 
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The  value  of  the  transverse  component  of  the  velocity  at  the  edge  of  the 
jet  (y  =  °°)  also  is  noteworthy.  From  equation  (9.7!)  one  finds  for 
this  lateral  inflow  velocity 


=  -  0.550 


(9.75) 


One  can  further  calculate  the  flow  volume  for  a  layer  of  unit  height 


One  finds 


Q 

P 


1/3 

3.3019  (££) 


^9.76) 


The  flow  volume  increases  downstream,  since  fluid  at  rest  is  carried 
along  from  the  side. 

The  solution  indicated  here  naturally  has  a  singular  point  at 
x  =  0,  since  an  infinitely  narrow  slot  with  infinitely  large  exit 
velocity  was  assumed.  Actually  for  a  narrow  slot  of  finite  width  one 
has  immediately  behind  the  slot  opening  a  velocity  distribution  that  is 
rectangular  across  the  Jet  cross  section  but  which  at  some  distance  is 
trarisformed  into  the  bell— shaped  distribution  found  here  with  width 

2/3  -  l/3 

b  -w  x  and  mid— velocity  ~  x 

Finally  it  should  be  mentioned  that  the  corresponding  rotationally— 
symmetrical  problem  where  the  jet  goe3  out  from  a  very  small  circular 
hole  also  can  be  solved  in  closed  form.  (compare  H.  Schlichting 
(reference  22)).  In  this  case  the  width  of  the  Jet  is  proportional  to  x 
arid  the  midvelocity  proportional  to  x~^. 


o.  The-  Boundary  Layer  for  the  Potential  Flow  U  =  u^xm. 

Another  class  of  exact  solutions  of  the  boundary  layer  equations 
will  b--  discussed  briefly  which  includes  the  plate  in  longitudinal  flow 
•t:  d  ih- •  stagnation  point  flow  as  special  cases.  Falkner  and  Shan 
'  ... u\ n  37)  have  shown  that,  just  as  for  Blasius*  plate  flow,  the 
:  ;n-y  lay-  r  di  ff.-r-ntial  •  quat  Ions  for  the  potential  flow 
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can  be  reduced  by  a  similarity  transformation  to  an  ordinary  differential 
equation  (u-^  =  constant,  m  >  0  accelerated,  m<  0  retarded  flow).  For 
m  >  0,  x  =  0  is  the  stagnation  point  of  the  flow.  For  m  =  0  one 
obtains  U  =  u-^  =  U0,  therefore  the  plate  flow;  m  =  1  gives  U  =  u-^x, 

therefore  the  stagnation  point  flow  according  to  equation  (6. If). 

The  differential  equations  of  the  boundary  layer  read 


u 


8u 

v  — 

8y 


U  ♦  V  — 

31  .Sy2 


8u  8v 
8x  8y 


0 


The  pressure  term  becomes 


u  an 

dx 


m  u 


2m-l 

x 


As  a  new  independent  variable  one  introduces 


i  = 


m  +  1 


m— 1 


(9-78) 


and  the  continuity  equation  is  integrated  by  Introduction  of  a  stream 
function  for  which  one  uses  the  equation: 


♦ 


\pl 


mfl 


) 


(9.79) 


One  has 


m— 1 

81  _  m  -  1  84  _  ,  I m  +  1  lu  1  x  2 

8x  2  x^  8y  y  2  \  v 


and  one  obtains 
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m 


u  =  u^x  qp*  ( I. )  =  U<p'(l) 


(9.80) 


'u. 


3at-l 


*  ■  rf^Vt1  v  2 


krb-2 


8  u  m  +  1  2  2  „»« 

-  =  -  u,  x  c  cp* 


Sy 


2  2 


CP*"(0 


^  =  u,  x111-1  J  m  cp'  +  - - -  |  cp" 

8x1  ]  2  5 


v  =  — 


m  + 


7 


a-1 
x  2  ) 


m  +  1  _  .  m  —  1  .. 
-  cp  +  -  i  cp* 


(9-81) 


After  insertion  into  the  first  equation  of  motion  and  division  by 
m  u-j^x2®-1  one  obtains,  when 


2m 

m  +  1 


the  following  differential  equation  for  cp(|): 

<p«"  =  -  qxp"  +  0  (cp*2  -  1) 


(9.82) 


(9.83) 


Boundary  conditions: 


|=0:  cp  =  cp  =  0;  |  =  00:  cp*  =  1 


The  differential  equation  (9.83)  was  solved  for  different  values  by 
Hartree  (reference  38).  The  result  is  given  in  figure  4la.  The 
corresponding  values  of  3  and  m  are  given  in  the  following  table. 
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m 

3 

-0.065h 

-0.1k 

0 

0 

0.111 

0.2 

0.333 

0.5 

1 

1 

k 

1.6 

For  accelerated  flow  (m  >  0,  B  >  0)  one  obtains  velocity  profiles 
without  inflection  points,  for  retarded  flow  (m  <  0,  B  <  0)  velocity 
profiles  with  inflection  points.  Separation  occurs  for 


3  =  -  0.199,  that  is,  m  =  -  0.091 


.  ,  -0.091 

Separation  takes  place  for  the  potential  flow  U(x)  =  u  x  ,  thuB 

for  very  weak  retardation.  Compare  chapter  XI  a,  where^almost  the  same 
result  is  obtained  from  an  approximation  calculation. 


Special  cases: 


1.  Stagnation  point  flow:  It  is  obtained  for 


m  =  1 ;  p  =  1 


Then  %  =  \j  —  y;  •vjr  =  \j\j  x  cp(  |) .  These  are  the  same  expressions  as 

for  the  stagnation  point  flow,  equation  (6.26a)  and  (6.27a),  also  (6.19) 
and  (6.26b).  The  differential  equation  (9*83)  also  j s  transformed  into 
the  equation  of  the  stagnation  point  flow  (equation  (6.28)). 

2.  Plate  in  longitudinal  flow:  This  case  is  obtained  for 


m  =  0; 


3  =  0 
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Then 


variable  according  to  equation  (9-3)' 


=  -2-  with  Ti 
Furthermore, 


signifying  Blasius' 
t  becomes 


i(r  =  \[2  \Jv\J  x  Cp(l):  thus  ®  =  —  compared  with  the  expression 

0  l2 

.  dtp  df 

for  \|r  for  the  plate  flow  equation  (9.4).  Because  of  —  =  — , 

1 —  dT) 

q>"(S)  =  \J  2  f  *'  (  tj  ),  and  cp*"(|)  =  2f*"(T)),  the  differential  equation  (9.83) 

is  for  this  case  transformed  into  2f'"(T))  +  ff"(T])  =  0.  This  is  identical 
with  equation  (9.8). 


CHAPTER  X.  APPROXIMATE  SOLUTION  OF  THE  BOUNDARY  LAYER  EQUATION 
BY  MEANS  OF  THE  MOMENTUM  THEOREM 
(KARMAN-POHLHAUSEN  METHOD,  PLANE  PROBLEM) 
a.  The  Flat  Plate  in  Longitudinal  Flow 


The  examples  of  exact  solutions  of  the  boundary  layer  equation 
discussed  in  the  previous  chapter  give  sufficient  proof  of  the  rather 
considerable  mathematical  difficulties  in  solving  the  differential 
equation.  Yet  the  examples  treated  were  selected  as  simple  as  possible. 

In  some  other  cases  the  mathematical  calculations  are  still  more  difficult 
than  in  those,  examples.  Particularly  the  problem  of  flow  about  a  body 
of  arbitrarily  prescribed  shape  cannot,  in  general,  be  solved  by  exact 
solution  of  the  differential  equations  of  the  boundary  layer.  Just  this 
problem,  however,  is  of  special  practical  importance,  for  instance  for 
the  calculation  of  the  boundary  layers  on  airplane  wings. 

There  exists  therefore  a  tendency  to  apply  at  least  approximate 
methods,  even  if  their  accuracy  is  sometimes  not  quite  satisfactory, 
for  cases  where  the  exact  solution  cannot  be  obtained  with  a  reasonable 
expenditure  of  calculation  time.  Such  simpler  approximate  solutions  can 
be  attained  if  one  does  not  attempt  to  satisfy  the  differential  boundary 
layer  equation  for  every  particle  of  fluid.  Instead  one  selects  a 
plausible  expression  for  the  velocity  distribution  in  the  boundary  layer 
and  satisfies  merely  the  momentum  equation  which  is  obtained  by  inte¬ 
gration  from  the  equation  of  motion.  A  method  on  this  basis  for  the 
plane  problem  of  flow  about  an  arbitrary  body  was  worked  out  by  von  Karman 
and  Pohihausen  (references  23  and  24).  We  shall  try  out  this  method  in 
this  chapter  at  first  on  the  simpler  case  of  the  flat  plate  in  longi¬ 
tudinal  flow,  where  no  pressure  variations  exist.  For  this  special  case 
the  momentum  theorem  yields  the  statement  that  the  momentum  loss  of  the 
flow  through  a  control  area  fixed  in  space  according  to  figure  42  equals 
thp  fried. ion  drag  W(x)  of  the  plate  from  the  leading  edge  (x  =  0)  to 
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the  point  x.  Application  of  the  momentum  theorem  for  this  case  was 
discussed  in  detail  in  chapter  IX;  for  the  drag  of  the  plate  wetted  on 
one  side  according  to  equation  (9.40)  it  had  resulted  in  the  formula: 


W(x)  =  h  p 


u  (U0  -  u)  dy 


Uy=0 


(10.1) 


On  the  other  hand  the  friction  drag  can  also  he  expressed  as  the  integral 
of  the  surface  shearing  stress,  namely: 


W(x)  =  b 


nx 


tq  (x)  dx 


Ux=o 


(10.2) 


In  forming  the  integral  (equation  (10.1))  it  is  to  be  noted  that  the 
integrand  outside  of  the  boundary  layer,  where  u  =  U0,  does  not  make 

a  contribution.  By  comparison  of  equations  (10. 1)  and  (10.2)  it  follows 
that : 


T 

O 


d 

P  — 
dx 


u  (U 


Uy=o 


u)  dy 


(10.3) 


If  one  introduces  moreover  the  momentum  loss  thickness,  as  defined  in 
equation  (9.42),  equation  (10. 3)  can  also  be  written  in  the  form: 


tt  2  di3  _ 

0  dx  "  p 


(10.4) 


This  is  the  momentum  theorem  of  the  boundary  layer  for  the  special  case 
of  the  flat  plate  in  longitudinal  flow.  Physically  it  states  that  the 
surface  shearing  stress  equals  the  momentum  loss  in  the  friction  layer, 
since  in  the  present  case  the  pressure  gradient  makes  no  contribution. 
The  next  chapter  will  acquaint  us  with  the  extension  of  equation  (10.4) 
to  include  the  more  general  case  o^  a  boundary  layer  tfith  pressure 
dif  f erence . 

Equations  (10.4)  and  (10. 3),  respectively,  will  now  be  used  for 
approximate  calculation  of  the  friction  layer  on  the  flat  plate  in 
longitudinal  flow.  A  comparison  of  the  results  of  this  approximate 
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calculation  with  the  exact  solution  according  to  chapter  IX a  will  give 
information  about  the  usefulness  of  the  approximation  method.  For  the 
approximate  calculation  one  selects  a  suitable  expression  for  the 
velocity  distribution  in  the  boundary  layer  in  the  form: 


u  =  UQ  f 


=  U0  f(ri) 


(10.5) 


with 


T )  =  -i  5 
1  5 


=  &(x) 


(10.6) 


5  represents  the  boundary  layer  thickness,  undetermined  at  first.  For 
the  flat  plate  it  may,  moreover,  be  assumed  again  that  the  velocity 
profiles  at  various  distances  from  the  leading  edge  of  the  plate  are 
affine  to  each  other.  This  assumption  is  contained  in  equation  (10.5) 
if  f(q)  there  stands  for  a  function  which  no  longer  contains  any  free 
parameters.  Furthermore,  f ( tj )  must,  for  large  values  rj ,  assume  the 
constant  value  1. 

The  velocity  distribution  being  given  by  equation  (10. 5),  the 
momentum  integral  in  equation  (10.3)  may  be  evaluated.  It  yields: 


u  (IT  -  u)  dy  = 


UD2  5 


f  (1  -  f)  dl) 


(10.7) 


The  definite  integral  in  equation  (10.7)  can  be  calculated  immediately 
if  a  definite  formulation  is  given  for  f  ( tj ) .  Thus  one  puts 


-  f)  cLtj 


(10.8) 


Hence 


2  2 
U  ■a  =  u(U  -  u)  dy  =  a  5  U 


(10.9) 


■3  =  a  5 


(10.10) 
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Furthermore,  the  displacement  thickness  of  the  boundary  layer  6*  becomes, 
according  to  equation  (9-20): 


P1 


5* 


( 1  —  f )  dn  = 


(10.10a) 


On  the  other  hand,  the  shearing  stress 


at  the  surface  is: 


To.  _  y  ( du  \ 

p  Vdy/ 

y=o 


f • (o)  =  3 


(10.11) 


if  one  introduces  the  further  simplification 


3  =  f*(0) 


(10.12) 


By  introduction  of  these  expressions  into  the  momentum  equation  (10.4) 
there  results: 


VU°  _  tj  2  dB 

“HI 


or 


d5  _  £  y_ 
dx  a  U 

o 


(10.13) 


(10.14) 


The  integration  with  the  initial  value  6=0  for  x  =  0  yields,  as 
first  result  of  the  calculation: 


6 


(10.15) 


Hence  the  shearing  stress  becomes,  according  to  equation  (10.11): 
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and  furthermore 


2 


and  hence,  finally,  the  total  drag  of  the  plate  wetted  on  both  sides 
according  to  equation  (10.2): 


2W 


2b  \[ZW  \jt 


MpUn32 


(10.17) 


A  comparison  of  the  results  for  boundary  layer  thickness,  surface  shearing 
stress,  ana  total  drag,  which  were  thus  found,  with  the  corresponding 
formulas  for  the  exact  solution  according  to  equation  (9*  •  * )  shows  that 
the  approximate  calculation  according  to  the  momentum  theorem  reproduces 
the  characteristics  of  the  formulas  with  perfect  correctness  in  all  cases, 
that  is,  the  dependence  of  the  boundary  layer  quantities  on  the  length 
of  run  x,  the  free  stream  velocity  UQ,  and  the  viscosity  coef¬ 
ficient  V.  The  numbers  a,  3,  unknown  at  first,  can  be  determined  only 
after  making  special  assumptions  for  the  velocity  distribution,  that  is, 
explicitly  prescribing  the  function  f(Tj)  in  equation  (10. 5). 


Numerical  examples 

The  usefulness  of  the  method  of  approximation  will  be  investigated 
by  a  few  numerical  examples.  The  accuracy  of  the  results  will  depend  to 
a  great  extent  on  a  suitable  choice  of  the  expression  for  the  velocity 
distribution  according  to  equation  (10.5).  At  any  rate  the  function  f ( -q ) 
must  equal  zero  for  q  =  0  and  have  the  constant  value  1  for  large  T). 
As  first  example  we  select  the  very  rough  assumption  that  the  velocity 
distribution  in  the  boundary  layer  is  represented  by  a  linear  function 
according  to  figure  h3a.  Thus: 


0  —  t)  —  1 :  f(n)  =  T) 
T1  -  1:  f(T))  =  1 


(10.18) 


Hence  the  results  for  the  numbers  a,  3  according  to  equation  (10.8) 
and  (10.12)  are: 


a 


3  =  1 


(10.19) 
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The  formulae  (10.17),  (10.16)  and  (10.17)  can  now  he  evaluated  immediately. 
One  obtains  the  results: 


6 


3.464 


(10.20) 


T 


O 


0.289 


(10.21) 


2W 


1.155 


(10.22) 


A  velocity  distribution  in  the  form  of  a  cubic  parabola  according  to 
figure  43b  is  selected  as  second  numerical  example  in  the  following  manner: 


0  -  T|  “  1:  f(Tl)  =  |  T  ~  |  I3 

t)  -  1:  f(Tj)  =  1 


(10.23) 


This  satisfies  the  conditions: 

T)  =  0:  f  =  0;  tj  =  1:  f  =  1;  f»  =  0; 


that  is,  the  boundary  layer  profile  joins  the  velocity  of  the  potential 
flow  with  a  continuous  tangent.  The  calculation  of  the  numerical  factors 
according  to  equations  (10.8)  and  (10.12)  gives: 


a 


(10.24) 


and  hence  for  the  characteristic  parameters  of  the  boundary  layer: 


6  =  4.64\fe 

K 


(10.25a) 


T 


O 


=  0.323  pU0 


(10.27b) 
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2W  =  1.29 


b 


2 


(10.25c) 


The  exact  value  for  the  drag  is,  according  to  equation  (9.18) 

2W  =  I.328  b  \jjJQ^idpl.  The  simple  assumption  of  a  linear  velocity  distri¬ 
bution  therefore  gives  a  drag  toe  small  by  thirteen  percent,  the  cubic 
velocity  distribution  a  drag  too  small  by  three  percent. 

The  calculation  of  the  displacement  thicknesses  of  the  boundary 
layers  according  to  equation  (10. 10a)  results,  for  the  linear  velocity 

distribution,  in  5*  =  —  6,  and  for  the  cubic  velocity  distribution 

3  ^ 

in  5*  =  g  5.  This  gives,  because  of  equations  (10. 15)  and  (10.25a): 


6*  =  1.732 

5*  =  1.740 


f 

^Vx 


(linear  velocity  distribution) 
(cubic  velocity  distribution) 

Vx~ 


>  (10.26) 


The  agreement  with  the  exact  value  5*  =  1.728  y—  according  to 

equation  (9. 21)  is  surprisingly  good;  this  iB,  however,  more  or  leBS 
accidental. 

The  essential  characteristics  of  the  boundary  layer  according  to 
the  approximate  calculation  described  above  are  once  more  compiled  with 
the  exact  solution  in  the  table  below. 


Characteristics  of  the  Boundary  Layer  on  the  Flat  Plate; 


Comparison  of  Approximate  Calculation  and  Exact  Solution 


Kind  of  calculation 

V  Vx 

_ 

V  vx 

5* 

d 

^(¥)1/2 

Linear  approximation 
(fig.  l+3a) 

0.578 

3.00 

0.289 

1.155 

Cubic  approximation 
(fig.  43b) 

0.645 

2.70 

0.323 

I.29 

Exact  solution 
(Blasius) 

1.729 

0.664 

2.61 

0.332 

1. 328 
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As  one  can  see  from  this  table,  the  agreement,  particularly  of  the  cubic 
approximation  and  the  exact  solution,  is  rather  good.  On  the  whole,  the 
results  of  this  calculation  with  the  aid  of  the  momentum  theorem  .may  be 
regarded  as  very  satisfactory,  especially  in  view  of  the  simplicity,  as 
compared  to  the  exact  calculation. 


Seventh  Lecture  (January  19,  1942) 

b .  The  Momentum  Theorem  for  the  Boundary  Layer 
with  Pressure  Drop  (Plane  Problem) 

Last  time  the  boundary  layer  on  the  flat  plate  in  longitudinal  flow 
was  calculated  by  means  of  the  momentum  theorem.  Today  the  general  case 
of  the  boundary  layer  with  a  pressure  difference  in  the  flow  direction 
will  be  treated.  One  considers  the  flow  along  a  curved  surface,  and 
measures  the  coordinate  x  as  arc  length  along  the  surface;  let  y  be 
the  perpendicular  distance  from  the  surface,  U(x)  the  given  potential 
flow  (fig.  44).  The  fundamental  equation  may  be  obtained  by  a  momentum 
consideration  as  in  chapter  IXc;  now,  however,  the  pressure  difference 
has  to  be  taken  into  consideration.  The  same  result  is  obtained  by 
integration  of  the  equation  of  motion  of  the  boundary  layer  with  respect 
to  y  from  y  =  0  (surface)  to  y  =  h,  the  layer  y  =  h  lying 
everywhere  outside  of  the  friction  layer  (fig.  44). 


The  differential  equations  of  the  boundary  layer  for  the  steady 
case  read,  according  to  equation  (8.3)} 


du  du  1  dp  , 

dx  dy  p  dx 

du  +  dv  _  0 

dx  dy 


d^u 

8y2 


(10.27) 


with  the  boundary  conditions:  y  =  0:  u 
integration  from  y  =  0  to  h  giveB: 


=  v  =  0;  y  =»  :  u  =  U.  The 


nh 


-  —  I  u2  dy  + 
2  dx  1 

lj=o 


du  , 
v  —  dy 

dy 


h  i  &  +  V 
p  dx 


du 

dy 


(10.28) 


In  the  first  term  the  differentiation  with  respect  to  x  and  the  integra¬ 
tion  with  respect  to  y  may  be  interchanged,  since  the  upper  limit  h 
is  independent  of  x.  On  the  left  side  the  second  term  is  transformed  by 
integration  by  parts: 
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vh  representing  the  transverse  velocity  outside  of  the  boundary  layer. 

By  continuity,  —  =  -  and 

Sy  Sc 


(10.29) 


and  hence: 


Insertion  in  equation  (10.28)  gives,  because  of: 


[-*  lh  t 

V  OU  _ _ o 

L^Jo  p 


the  relation: 


(10.30) 


(10.31) 


This  is  the  so-called  Karman  integral— condition,  first  given  by  v.'Karman 
(reference  23).  ' 

For  the  pressure  term  one  now  introduces  the  potential  velocity  U(x) 
furthermore,  equation  (10. 31)  is  to  be  transformed  so  that  the  displacem-n 
thickness  and  the  momentum  thickness  $  appear  in  if  as  defined  by 

equation  (6.31)  and  (6.32),  namely: 
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Hh 


U  5*  = 


(U  -  u)dy 


2 

Ul  = 


nh 


u  (U  -  u)  dy 


(10. j2) 


(10.37) 


According  to  Bernoulli’s  equation: 


1  dp  _  _  dU 
^  cbc  dx 


dU“  dU 

-  +  U  — 

dx  dx 


which  can  also  he  written: 


h  dp _ d_ 

p  dx  dx 


U2  dy  +  U  ~ 
dx 


Jo 


Jo 


U  dy 


(10.34) 


By  substitution  of  equation  (10.34)  into  equation  (10. 31)  there  results: 

th  flh 

T 

—  =  —  (U2  -  u2)  dy  -  IJ  —  /  (U  -  u)  dy 

p  dx  dx 


_4 

dx 


(U  -  u)  u  dy  +  —  \  U 
dx 


J. 


(U  -  u)  dy  >  -  U 


dx 


nh 


(U  -  u)  dy 


and  after  differentiation  of  the  second  term: 

nh  rh 


To  _  d 
p  dx 


(U  -  u)  u  dy  +  — 
dx 


( U  -  u )  dy  (10. 


Jo 
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The  displacement  thickness  and  the  momentum  thickness  can  now  be  introduced 
directly  and  one  obtains: 

—  =  —  (U2*)  +  —  &*  U 
p  dx  dx 


or 


2g+  (2S  +  6.)U§ 


(10.36) 


This  is  the  form  of  the  momentum  equation  for  the  boundary  layer  with 
pressure  drop  that  will  be  used  as  a  basis  for  further  considerations. 
Since  in  it  t„  is  quite  generally  the  surface  shearing  stress, 
equation  (,10.3q)  must  apply  in  the  same  way  to  turbulent  flow,  too.  We 
shall  come  back  to  that  later.  For  the  special  case  of  vanishing  pressure 

drop  —  5  0,  equation  (10. 36)  is  transformed  into  equation  (10.4)  found 
dx 

before  for  the  flat  plate  in  longitudinal  flow. 

The  further  calculation  of  the  boundary  layer  on  the  basis  of 
equation  (10. 36)  is  performed  for  the  laminar  case  according  to  the  method 
of  Pohlhausen  (reference  24)  and  for  the  turbulent  case  according  to  the 
method  of  Gruschwitz  (reference  34)  (chapter  XVIII). 


c.  Calculation  of  the  Boundary  Layer  According  to  the 


Method  of  Karman-Pohlhausen— Holstein 


For  further  calculation  it  is  of  importance  to  find  a  suitable 
expression  for  the  velocity  distribution  in  the  boundary  layer  u(x,  y). 
According  to  our  understanding  of  the  exact  solutions  of  the  differential 
equations  of  the  boundary  layer  this  expression  must  at  least  satisfy 
the  conditions  that  for  y  =  0:  u  =  0,  and  for  y  =  00 :  u  =  U.  Further¬ 
more  the  derivative  du/dy  must  vanish  for  large  y.  Moreover,  velocity 
profiles  with  and  without  inflection  points  must  be  possible,  as  they 
occur  in  the  pressure  decrease  and  pressure  increase  region,  respectively. 


Finally,  a  profile  with 


must  be  possible  in  order  to  have  a 


separation  point  result  from  the  approximate  calculation. 


One  chooses  for  the  velocity  distribution  an  expression  of  the  form 
u  (x,  y)  =  Uf(y/5p),  and  sets,  according  to  Pohlhausen, 
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valid  for  0  —  y  —  5pi  0  —  r\  —  1.  6p(x)  stands  for  the  boundary  layer 

thickness,  the  dependence  of  which  on  x  has  yet  to  he  calculated.  The 
boundary  layer  thickness  of  the  approximate  calculation  6p  is  here 

provided  with  the  index  P  (  =  Pohlhause^  in  order  to  avoid  confusion 
with  the  boundary  layer  thickness  5  used  before.  Whereas  for  the 
exact  solutions  the  velocity  in  the  boundary  layer  asymptotically 
approaches  the  velocity  of  the  potential  flow,  (u  — »U  for  y  — ><»  ); 
the  value  u  =  U  is  to  be  attained  in  the  approximation  at  a  finite 
distance  from  the  surface  y  =  6p,  for  reasons  of  calculation.  This 

modification  of  the  actual  relation  is  physically  insignificant. 

For  the  determination  of  the  free  constants  a,  b,  c,  d  in 
equation  (10.37)  the  following  boundary  conditions  are  prescribed,  all 
of  which  follow  from  the  differential  equation  of  the  boundary  layer 
(equation  (8.3)): 


(10.38) 


Since  the  condition  of  no  slip  u  =  0  for  y  =  0  is  automatically 
satisfied  by  expression  (10. 37),  the  four  free  constants  a,  b,  c,  d 
are  sufficient  to  satisfy  the  remaining  four  conditions.  The  last 
of  the  five  conditions  follows  immediately  from  the  exact  differential 
equation  of  the  boundary  layer  if  one  puts  y  =  0  and  takes  the 
boundary  conditions  into  consideration.  This  condition  is  particularly 
Important  since  it  determines  the  curvature  of  the  velocity  profiles 
near  the  surface  and  assures  that  boundary  layer  profiles  do  not  acquire 
an  inflection  point  in  the  region  of  pressure  decrease  and  do  acquire 
one  in  the  region  of  pressure  increase,  as  required  by  the  exact  solution 
according  to  chapter  VTTI.  From  equation  (10. 38)  follows  for  the 
coefficients  a,  b,  c,  d  the  equation  system: 
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Hence  the  expression  for  the  velocity  distribution,  which  satisfies  all 
boundary  conditions  according  to  equation  (10.37)>  reads: 

“  ~  (2t)  -  2t)3  +  tj1*)  +  |  (tj  -  3t)2  +  3t]3  - 

(10.J+2) 

g  =  f(h)  +  \  g(t)) 

in  which 

F(ri)  =  2i)  -  2t)3  +  q4 
G(tj)  =  5  (tj  -  3t12  +  3l3  -  V*) 
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Due  to  T|  =  y/&p(x)  the  boundary  layer  thickness  8p(x)  is  here  the 
only  unknown.  If  that  is  calculated,  the  parameter  X  follows  immediately 
from  equation  (10.41).  From  equation  (10.42)  one  understands  that  tlje 
velocity  profiles  form  with  the  form  parameter  \(x)  a  one— parameter 
family.  The  functions  F(tj)  and  G(tj)  indicated  in  figure  45  and 
table  4  have  a  universal  character,  that  is,  they  do  not  depend  on  the 
special  body  shape.  The  velocity  profiles  for  various  values  of  X,  are 
plotted  in  figure  46.  The  profile  with  X  =  0  is  obtained  for  dU/dx  =  0, 
that  is,  for  the  boundary  layer  without  pressure  gradient  (flat  plate  in 
longitidinal  flow).  The  separation  profile  with  (du/dj^  =  0,  that  is, 

a  =  0,  has  according  to  equation  (10.40b)  the  parameter  A.  =  -12.  The 
profile  at  the  stagnation  point  has,  as  will  be  shown  below,  X  =  7.052. 

For  X  >  12  there  result  values  of  u/u  >1  in  the  boundary  layer,  which 
physically  does  not  make  sense.  These  values  therefore  have  to  be  excluded. 
Since  behind  the  separation  point  the  boundary  layer  calculation  loses  its 
validity  anyway,  the  form  parameter  X  is  limited  to  the  region 

-  12  -  X  -  +  12  (10.44) 


The  unknown  boundary  layer  thickness  5p(x)  remains  to  be 

calculated.  For  this  the  momentum  equation  (10. 36),  so  far  not  utilized,  is 
at  our  disposal.  Before  performing  this  computation,  a  few  preparatory 
calculations  are  required,  namely  the  determination  of  the  boundary  layer 
characteristics,  displacement  thickness  5*,  momentum  loss  thickness  d, 
and  surface  shearing  stress  tq  on  the  basis  of  the  approximation- 

expression  equation  (10.42):  One  obtains  from  equations  (10. 32)  and 
(10.42): 


o* 

St, 


[l  “  F(tj)  -  XG(tj)]  dq 


*  T|=0 


|1 

[f(t))  +  XG(r))]  [l  -  F(tj)  -  \G(tOJ  dq 

Jrj=o 


The  calculation  of  the  definite  integrals,  with  the  values  of  F(q)  and 
G(t|)  according  to  equation  (10. 43),  gives: 


S^  _  J _ x~ 

Sp  10  120 


(10.45) 
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i_  =  _37  _  _X_  _  \g 

Sp  315  ~  9l5  9072 


(10.16) 


Further,  there  results  for  the  surface  shearing  stress  from  = 
according  to  equations  (10.1-2)  and  (10. I3) 


,  8u  , 


\  s  12  +  X 
H  U  6 


(10.17) 


Now  the  momentum,  equation  (10,36)  is  to  be  used  for  calculation  of 
Sp(x).  After  multiplication  by  d/vU  it  acquires  the  dimensionless  form: 


V  d  /  v  a  u 


(10.I8) 


The  boundary  layer  thickness  Bp  does  not  even  appear  in  this  equation; 

however,  this  is  not  particularly  astonishing,  since  Bp  is  a  rather 

arbitrary  quantity  of  our  approximate  calculation  and  therefore  without 
special  physical  significance.  The  physically  important  quantities, 
displacement  thickness  5*  and  momentum  loss  thickness  d,  appear 
instead  in  equation  (10.18).  Hence  it  suggests  itself  to  first  , 
calculate  d  from  the  momentum  equation  (10.18)  and  then  to  pass  on 
to  Bp  by  means  of  equation  (10.16).  For  this  purpose  one  introduces 

according  to  Holstein  and  Bohlen  (reference  25)  aside  from  the  first  form 
parameter  X  according  to  equation  (10. 11)  a  second  form  parameter  k, 
formed  analogously  with  the  momentum  thickness 


k  =  U* 
v 


(10.19) 


Then  one  sets: 


z  =  £■ 
v 


(10.50) 


Then 


K  =  Z  U* 


(10.51) 
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Between  the  second  form  parameter  k  and  the  first  form  parameter  \ 
exists,  according  to  equations  (10.40)  and  (10.41)  the  universal  relation: 


K  -MI- 4--  A_ 

\315  9^5  9072 


(10.52) 


Furthermore,  for  simplification,  one  substitutes: 

_ iL_ 

fii- - - mf(K) 

__n  _  _x _ L. 

315  9^5  9072 


(10.53) 


l£  1  =  l£  M  ±  =  /12  +  37  _  _x _ x 


n  u  u  u  8P  \  6  A315  945  9072/  2 


=  f0(K) 


(10.51+) 


By  introduction  of  k  and  Z  according  to  equations  (10.49)  and  (10. 50) 
and  by  substituting  from  equation  (10.53)  and  (10.54)  one  now  obtains 

from  the  momentum  equation  (10.48),  because  of  dd*/v  *  ~  ~ 


|u  g  +  j^2  +  f^K) ijie  -  f2(>c)  =  0 


(10.55) 


Finally,  one  sets  as  further  simplification 


2f2(«)  -  4k  -  2fx(K)  k  =  F(k) 


(10.56) 


or  written  in  detail: 


F(k)  =  2(  3§  -  955  -  ^72  )  f 


116  /  2  _i\ 

315  K  +  \9^+5  120/ 


2  P  3 

x  +  §575  1  f(10-57) 


and  thus  obtains  the  following  differential  equation  for  the  momentum 
thickness 


dZ  =  F(k  )  n  =  Z  U* 
'  dx  U 


(10.58) 
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This  is  a  non-linear  differential  equation  of  the  first  order 
for  z  =  $2/v.  The  fact  that  the  function  F(k)  is  rather  complicated 
does  not  constitute  any  appreciable  drawback,  since  F(k)  is  universal, 
that  is,  independent  of  the  shape  of  the  body,  and  hence  may  be  tabulated 
once  and  for  all.  The  functions  F(k),  f^(*),  fg^)  as  aleo  K  =  *(*•) 
according  to  equation  (10.52)  are  given  in  table  5. 

As  to  the  solution  of  equation  (10. 58)  the  following  remains  to  be 
said:  The  calculation  has  to  start  at  the  stagnation  point  x  =  0.  There 

U  =  0;  and  the  initial  slope  —  would  be  infinite  if  F(k)  were  not 

dx 

also  equal  to  zero  at  the  stagnation  point.  The  function  F(K)  actually 
has  a  zero  which  yields  a  physically  significant  initial  value.  This 
zero  of  F(k)  is  given  when  the  second  bracket  in  equation  (10.57) 
disappears.  One  finds: 


F(k)  =0  for  k  =  Kq  =  0.0770;  X.  =  XQ  =  7.052  (10.59) 


The  value  \  =  7*052  therefore  gives  the  value  of  the  first  form 
parameter  at  the  stagnation  point.  Then  the  initial  slope  of  the  integral 
curve  at  the  stagnation  point  now  has  the  indeterminate  value  The 

latter  may,  however,  be  calculated  and  hence  finally  yields  the  initial 
value  and  the  initial  slope  of  the  integral  curve  as: 


z  =  — a  =  Q«Q77Q 
0  »£ 


U" 

=  -0.0652  -2- 


> 


(10.60) 


The  index  o  denotes  the  values  at  the  stagnation  point.  With  these 
initial  values  one  succeeds  easily  in  performing  the  integration  of 
equation  (10. 58),  for  instance,  according  to  the  isocline  method.  A 
calculation  example  is  given  in  the  appendix,  figure  k7,  and  table  6. 

The  calculation  is  to  be  carried  up  to  the  separation  point  \  =  —12: 

K  =  —  O.I567.  Quantities  entering  the  calculation  that  are  given  by  the 
potential  flow  are  the  velocity  U(x)  and  its  first  derivative  with 
respect  to  the  arc  length  dU/dx.*  (Only  at  the  stagnation  point  is 
d%/dx2  also  required  for  the  initial  slope  of  the  Integration  curve.) 

*  In  PohlhauSen*s  treatment  (reference  2k)  a  differential  equation  is 
obtained  instead  of  equation  (10. 58),  for  the  quantity  z  =  5p'/v,  formei 
analogously  to  Z.  Pohlhausen’s  differential  equation  also  contains 
d2U/dx2  which  often  can  be  obtained  from  the  given  potential  flow  only 
by  a  double  graphical  differentiation.  The  representation  of  Holstein 
which  completely  avoids  the  quantity  d2U/dx2  therefore  means  an 
essential  improvement  of  the  method. 
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The  values  of  the  form  parameters  for  the  three  special  cases: 
stagnation  point,  velocity  maximum  (U*  =  0)  (pressure  gradient  equals 
zero),  and  separation  point  are  compiled  in  the  table  below. 


Case 

X 

K 

Stagnation  point 

7.052 

0.0770 

Velocity  maximum 

0 

0 

Separation  point 

-12 

-O.I567 

The  entire  process  of  calculation  takes  the  following  course: 

1.  The  integration  of  equation  (10. 58)  yields  Z(x),  k(x)  and 
according  to  equation  (10.50),  also  d(x);  furthermore  it  yields  the 
position  of  the  separation  point. 

2.  First  form  parameter  X(x)  from  equation  (10. 52) 

3.  Displacement  thickness  5*  from  equation  (10. 53) 

4.  Surface  shearing  stress  tq  from  equation  (10. 54) 

5.  Boundary  layer  thickness  5p(x)  from  equation  (10.45) 

6.  Velocity  distribution  u/U  from  equation  (10.42) 

Flat  Plate  in  Longitudinal  flow 


The  special  case  of  the  flat  plate  in  longitudinal  flow  which  was 
treated  in  chapter  la  with  a  different  form  for  the  approximation  can 
also  be  obtained  very  simply  from  the  present  calculation.  U  =  U  ; 

U*  s  0,  and  hence  k  =  X  =  0,  and  according  to  equation  (10. 58): 

dZ  _  F(0)  _  0.4698 
Uo  Uo 


With  the  initial  value  6=0  for  x  ~  0  there  results 

62  =  O.4698  ^  or  6  =  O.685 \/^- 

o  V  0 


(10.61) 


NACA  TM  No.  1217 


> 


101 


whereas  the  exact  value  according  to  BlasiuB*  calculation,  equation  (9.21b), 
is  d  =  0.66k\jvx/Uo.  Furthermore,  it  follows  that  th9  displacement 
thickness  is  with  b*/$  -  f^O)  =  2.54 


6*  =  1.75 


(10.62) 


The  shearing  stress  becomes,  from  equation  (10.54)  with  f2(0)  =  0.235; 


T 


O 


=  0.343  tiU0 


(IO.63) 


while  the  exact  value  tq  is  tq  =  0.332  fiU0  '/u0/vx  according  to 
equation  (9.17).  The  agreement  with  the  exact  values  is  rather  satisfactory. 
In  figure  50b  the  velocity  distribution  obtained  by  the  approximate  calcula¬ 
tion  also  is  compared  with  the  exact  calculation  in  the  plot  u/U  against 
y/8*.  This  agreement  also  is  rather  good. 


Stagnation  Point  Profile 

A  similar  comparison  can  be  performed  for  the  stagnation  point 
profile  the  exact  solution  of  which  was  given  in  chapter  VI.  For  this 
case  \  =  XQ  =  7.052,  it  =  kg  =  0.0770.  For  the  approximate  calculation 
the  momentum  thickness  is: 


-  Vs  -  °-2f8 


(10.64) 


whereas  according  to 
t3  \|u*  /v  =0.292.  The 
calculation; 


the  exact  calculation  equation  (6.34)  it  is 
displacement  thickness  is,  for  the  approximate 


S*\|F=  Vs  =  °-6u 


(10.65) 


whereas  according  to  the  exact  calculation,  equation  (6.33),  it  is 

5*  \jv  *  /v  =  0.648.  Finally,  according  to  the  approximate  calculation  the 
sheetring  stress  is: 
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T_o  =  i2_lo_  =  0,.332  =  1<19 

nu-vu*  0.278 


(10.66) 


T  I  V 

compared  with  \J—  =  <p”(0)  =  1.234  "by  the  exact  calculation.  Thus 
the  agreement  with  the  exact  values  is  here  also  satisfactory. 

The  velocity  distribution  of  the  approximate  calculation  is  compared 
with  the  exact  calculation  in  figure  ^0—a.  in  the  plot  u/U  against  y/5*. 
Here  also  the  agreement  is  good. 

The  following  table  contains  a  compilation  of  the  comparisons  Just 
given,  between  the  characteristics  5*,  d ,  t0  from  exact  and  from 
approximate  calculation. 


Blasius  Profile 

Stagnation  point  profile 

8*  \F° 

to 

&* 

5#  [f^L 

To 

5X 

d 

Vvx 

M  \jvo3/vx 

d 

*\|  V 

uU\fv/U» 

Pohlhausen 

Approximation 

1.75 

0.685 

O.3I+3 

2.55 

0.641 

0.278 

1.19 

2.31 

Exact  solution 

1.73 

0.664 

0.332 

2.61; 

ti 

0.648 

0.292 

1.234 

2.21 

Of  course,  it  can  not  yet  be  concluded  from  this  good  agreement  of 
the  approximate  with  the  exact  solution  that  similar  good  agreement  would 
exist  for  all  the  boundary  layer  profiles  along  the  body.  Accurate 
comparisons  are  not  easily  performed  since  very  few  exact  solutions 
reaching  from  the  stagnation  point  to  the  separation  point  exist.  However 
one  may  conclude  from  occasionally  made  comparisons  that  in  the  region  of 
pressure  decrease  the  agreement  is  mostly  rather  good;  in  the  region  of 
pressure  increase,  particularly  near  the  separation  point,  some  deviation 
might  occur. 

Since  no  other  serviceable  methods  for  boundary  layer  calculation 
have  so  far  become  known, the  Pohlhausen  method  is  for  the  present  to  be 
regarded  as  the  best.  The  time  required  for  a  boundary  layer  calculation 
for  one  side  of  the  body  immersed  in  a  given  potential  flow  amounts  to 
about  three  hours 4 

The  calculation  described  here  for  the  plane  flow  was  applied  by 
Tomotika  (reference  26)  to  the  rotationally— symmetrical  case. 
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d.  Examples 

A  few  examples  of  'boundary  layer  calculations  will  be  given,  all 
of  which  were  performed  according  to  the  approximation  method  described 
in  the  previous  section. 

The  first  example,  taken  from  a  treatise  of  H.  Schlichting  and 
A.  Ulrich  (reference  35)  concerns  a  series  of  elliptical  cylinders  in  a 
flow  parallel  to  the  major  axis.  The  axis  ratios  are  a^/b-^  =  1 

(circular  cylinder),  2,  4,  8,  »  (flat  plate).  The  velocity  distributions 
of  the  potential  flow  are  given  in  figure  51 •  For  slender  cylinders  they 
are  characterized  by  the  fact  that  the  velocity  has  a  very  flat  maximum 
and  decreases  steeply  toward  the  front  and  rear.  The  result  of  the 
boundary  layer  calculation  is  illustrated  in  figure  52.  Figure  52-a 
shows  the  dimensionless  displacement  thickness  as  a  function  of  the  arc 
length.  It  is  noteworthy  that  for  all  elliptic  cylinders  the  boundary 
layer  on  the  front  half  is  smaller  than  for  the  flat  plate;  this  is 
caused  by  the  pressure  decrease.  The  laminar— separation  point  for  the 
circular  cylinder  lies  at  s/t  =  0.604  (cp  =  109. 5°).  With  increasing 
fineness  ratio  it  shifts  further  toward  the  rear;  it  is  also  plotted  in 
figure  51.  Figure  52-b  shows  the  variation  of  the  form  parameter  X  and 
figure  52-c  the  variation  of  the  surface  shearing  stress.  For  every 
elliptic  cylinder  the  latter  has  a  maximum,  the  position  of  which  shifts 
frontward  with  increasing  fineness  ratio.  The  variation  of  the  shearing 
stress  for  the  cylinder  of  axis  ratio  eight  shows  only  an  insignificant 
difference  from  the  one  for  the  flat  plate  in  longitudinal  flow.  Figure  53 
gives  a  survey  of  the  variation  of  the  boundary  layer  on  the  body  and  the 
velocity  profiles  at  various  points  along  the  surface.  Corresponding 
calculations  for  the  rotaticnally— symmetrical  case  (that  is,  for  ellipsoids 
of  revolution  with  flow  parallel  to  the  axis)  have  been  performed  by 
Pretsch  (reference  27). 

The  second  example  (reference  35)  gives  the  boundary  layer  on  a 
symmetrical  Joukowski  profile  of  fifteen  percent  thickness  for  lift 
coefficients  in  the  region  of  ca  =  0  to  1.  Figure  5^  gives  the  velocity 

distribution  according  to  the  potential  theory  for  the  lift  coefficients 
ca  =  0  ,  0.25  ,  0.50,  and  1.0.  For  the  symmetrical  approach  flow 

(ca  =  0)  the  velocity  maximum  lies  rather  far  toward  the  front  at 

s/t  =  0.141.  With  growing  lift  coefficient  the  velocity  maximum  increases 
on  the  suction  side  and  decreases  on  the  pressure  side.  Simultaneously 
the  maximum  shifts  farther  forward  on  the  suction  side  and  farther  rearward 
on  the  pressure  side.  The  magnitude  of  the  velocity  maximum  and  its 
position  are  of  primary  importance  for  the  development  of  the  boundary 
layer  and  in  particular  for  the  position  of  the  separation  point.  The 
results  of  the  boundary  layer  calculation  are  plotted  separately  for  the 

suction  and  pressure  sides  in  figure  55.  Figure  55— a  gives  the 
variation  of  the  displacement  thickness,  figure  55— b  the  form 

parameter  \,,and  figure  55-c  the  surface  shearing  stress.  The  position 
of  the  laminar  separation  point  is  plotted  in  figure  5^.  Figure  56  gives 
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a  survey  of  the  variation  of  the  boundary  layer  along  the  surface  and  of 
the  velocity  distribution  in  the  boundary  layer  for  c&  =  0,  figure  57 
for  ca  =  1. 

Finally,  the  third  example  gives  a  survey- of  the  influence  of  the 
most  important  profile  parameters  of  a  wing  profile  on  the  laminar 
boundary  layer.  K.  Bussmann  (reference  36)  performed  the  boundary  layer 
calculation  for  a  family  of  Joukowsky  profiles  of  relative  thickness 
d/t  =  0  to  O.25  and  relative  camber  f/t  =  0  to  0.08  for  ca—  values 

from  0  to  1.  Figure  58  shows  the  family  of  Joukowsky  profiles.  Of  the 
very  voluminous  results,  only  the  position  of  the  separation  point  shall 
be  shown  here.  Figure  59  shows  the  position  of  the  separation  point  on 
the  suction  side  as  a  function  of  thickness,  camber,  and  lift  coefficient; 
figure  60  shows  the  same  result  for  the  pressure  side. 

Herewith  the  discussion  of  the  approximation  method  for  calculation 
of  the  laminar  boundary  layer  will  be  concluded. 


Eighth  Lecture  (January  26,  19^2) 
CHAPTER  XI.  PREVENTION  OF  SEPARATION 


For  practical  flow  problems  the  flow  with  pressure  increase 
(retarded  flow)  plays  an  important  role.  It  is  always  desirable  that 
no  separation  of  the  flow  from  the  wall  occur,  because  of  the  resulting 
large  losses  in  energy.  The  wing  presents  a  good  example.  A  pressure 
inert  use  exists  on  the  suction  side  toward  the  trailing  edge  (fig.  6l). 

If  separation  occurs,  the  wing  will  have  an  'undesirably  large  drag  and 
small  lift.  Another  example  is  the  flow  in  an  expanding  passage 
(diffuser)  which  transforms  kinetic  energy  into  pressure  energy  (as  for 
instance  in  the  wind  tunnel  or  in  the  bucket  grid  of  a  turbine). 

Calculations  will  presently  show  that  the  ability  of  a  laminar  flow 
to  overcome  a  pressure  increase  without  separation  is  exceedingly  small. 
Thus  the  pressure  increases  present  in  practical  flows  would,  for  laminar 
flow,  almost  always  lead  to  separation.  The  reason  that,  nevertheless, 
in  many  cases  of  practical  flows  considerable  pre33ure  increases  are 
surmounted  without  separation  is  that  the  flow  is  turbulent.  As  we  shall 
see  more  clearly  later,  the  ability  to  overcome  a  pressure  increase  with¬ 
out  separation  is  very  much  greater  for  turbulent  than  for  laminar  flow. 
Since,  moreover,  the  pressure  increase  always  gives  rise  to  an  early 
transition  from  laminar  to  turbulent*,  one  has  to  deal  almost  exclusively 
with  turbulent  flow  in  practical  flows  with  pressure  increase. 

Nevertheless  it  is  useful  to  clarify  the  fundamental  relations 
regarding  prevention  of  separation  for  laminar  flow,  particularly  because 
Its  phenomena  ler. i  themselves  more  readily  to  numerical  treatment  than 
those  of  the  turbulent  Plow. 


♦Details  are  given  in  chapter  XXI  "Transition  from  Laminar  to  Turbulent." 
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Various  possibilities  exist  for  prevention  of  separation.  The 
simplest  way  is  to  make  the  pressure  increase  so  small  that  separation 
is  avoided.  A  numerical  estimation  in  the  next  section  will  give 
information  about  this  possibility.  Another  possibility  consists  in 
artificially  influencing  the  boundary  layer,  for  instance  by  blowing 
or  suction  of  fluid,  or  else  by  application  of  an  auxiliary  wing  that 
provides  acceleration  at  the  critical  points  of  the  boundary  layer. 
Some  details  will  be  given  in  the  following  sections  of  this  chapter. 


a.  Estimation  of  the  Admissible  Pressure  Gradient 

We  are  going  to  make,  following  Prandtl  (reference  2),  a  generally 
valid  estimation  of  the  pressure  increases  in  a  laminar  boundary  layer 
that  are  possible  without  the  occurrence  of  separation.  We  take  as  the 
basis  the  Karmen— Pohlhausen  approximate  calculation  according  to 
Chapter  X  and  make  the  assumption  that  under  the  effect  of  the  pressure 
gradient  given  by  the  potential  flow  the  boundary  layer  has  developed 
till  near  the  separation  point  (Point  0  in  fig.  62).  From  here  on 
the  pressure  distribution  is  to  be  such  that  the  form  of  the  velocity 
profile  does  not  change  further  downstream,  that  is,  the  form  parameter 
X  is  to  remain  constant.  Since  the  value  corresponding  to  the  separation 
point  is  X  =  —12,  this  constant  X— value  shall  be  chosen  at  X  =  —10*. 

A  definite  value  of  the  second  form  parameter  (according  to  table  5) 
corresponds  to  this  choice: 

X  =  -10;  k  =  -0.1369;  F(k)  =  1.523  (11.1) 


For  the  prevention  of  separation  the  following  relation  between  the 


potential— flow  velocity  U(x)  and 
according  to  equations  (10.50)  and 


IT  =  0.1369 


the  momentum  thickness  d(x)  results 


(10.51): 

0.1369 

-u» 

(11.2) 

1369  ■— 

U,c? 

(11.3) 

=  0.1369  0 

(11.4) 

if  one  put 8  for  simplification 


*At  any  rate,  this  X— value  must  be  negative,  since  otherwise  the  flow 
In  question  is  not  retarded. 
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a 


(11.5) 


On  the  other  hand,  the  momentum  equation  according  to 
equation  (10. 58)  holds  for  the  further  development  of  the  boundary  layer 
for  x  >  0: 


u  —  =  f(k )  =  f(-o.1369)  =  1.523 

dx 


(11.6) 


The  numerical  value  mast  be  substituted  for  F(tc),  if  the  form  parameters 
are  to  remain  constant  at  the  values  given  by  equation  (11. 1).  From 
equations  (11.6)  and  (11.4)  follows  therewith,  for  the  constancy  of  the 
form  parameter  X  =  —10,  the  conditional  equation 


O.I369  =  1.523 

U» 


or 


a  = 


UU" 


U» 


=  11.13  21  11 


(11.7) 


For  cr  >  11  the  boundary  layer  can  still  bear  the  pressure  increase; 
for  a  <  11  separation  occurs;  for  0  =  11  the  boundary  layer  always 
remains  with  X  =  —10,  on  the  verge  of  separation.  Qualitatively,  the 
following  can  be  immediately  said  about  the  distribution  of  the  potential- 
flow  velocity  U(x)  which  gives  no  separation.  Because  of  equation  (11.7) 
a  necessary  condition  for  avoiding  separation  in  retarded  flow  is; 


U"  >  0  for  u*  <  0  (11.7a) 

that  is,  a  negative  velocity  gradient  U’  must  exist,  the  magnitude  of 
which  decreases  in  the  flow  direction.  If,  therefore,  the  curve  U(x) 
in  figure  63  is  curved  downward  behind  the  maximum  (U*  <  0),  separation 
occurs  in  every  case;  if  it  is  curved  upward  (U"  >0),  separation  some¬ 
times  does  not  occur.  The  limiting  case  U"  =  0  for  U*  <  0  always 
leads  to  separation.  The  sufficient  condition  for  avoiding  separation 
is  UU"/(J‘2  >  11. 
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One  now  proceeds  to  calculate  what  potential  flow  and  what  boundary 
layer  thickness  variation  correspond  to  a  =  +11.  From  equation  (11.7) 
follows : 


£  =  XI  UL 
U'  0 


and  after  integration:  log  U*  =  11  log  U  —  log  C|  or 


E-.-c* 
u11  1 


with  as  integration  constant.  Repeated  integration  gives: 

—  U~10  =  C*  x  +  C0  (11.8) 

10  1  2 

For  x  =  0,  TJ(x)  shall  be  U(x)  =  U0,  thus 


(11.9) 


Furthermore,  one  puts 


10 

C[  Uo  =  Cx 


(11.10) 


and  obtains  from  equation  (11.8)  for  the  potential  flow 


U  = 


XL 


(1  +  10  C-ji) 


0.1 


(11.11) 


Thereby  is  found  the  desired  velocity  distribution  that  Just  avoids 
separation.  The  constant  C-j_  can  be  determined  from  the  boundary  layer 
thickness  SQ  at  the  initial  point  x  =  0: 
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According  to  equation  (11.11) 


U*  =  - 


C1  Uo 


(1  +  10  C-ji)' 


1.1 


and  thence 


5  = 


\i 


CTT  (i  +  10  V)0'55 

1  O 


From  5  =  6C  for  x  =  0  follows 


C 


1 


=  10.  v 

2 

V>o 


(11.12) 


and  thus,  as  the  final  solution  for  the  potential  flow  and  the  boundary 
layer  thickness  variation. 


+  ioo 

Vo/ 


-o.i 


(11.13) 


5  = 


B0/l  +  100 


0.55 


U050 


(11.14) 


The  permissible  retardation  (velocity  decrease)  Is  therefore  comparable 
10__ 

to  1 / \y  x  and  Is  thus  very  small.  The  velocity  is  thus  very  close  to 
the  constant  velocity  of  the  flat  plate  in  longitudinal  flow.  For  the 
present  case  the  growth  of  the  boundary  layer  thickness  5  must  therefore 

1/2 

be  somewhat  larger  than  for  the  flat  plate,  where  5  ~  x  .  Here 
0.55 

5  *'*  x  *  ;  thus  the  increase  is  only  slightly  larger. 

The  flow  in  a  divergent  channel  with  plane  walls  (two-dimensional 
problem)  will  be  treated  as  another  example.  In  figure  6k  let  x  be 
the  radial  distance  from  the  origin  0.  The  walls  start  at  x  =  a,  where 
the  entrance  velocity  of  the  potential  flow  equals  U0.  The  potential 
flow  is 


♦Compare  Chapter  IX  e  where  it  was  found,  as  exact  solution  of  the 

differential  equation  of  the  boundary  layer,  that  in  retarded  flow 

.  .  ,  Tr ,  .  — 0.091 

separation  occurs  when  U(x)  =  x 


* 
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U(x)  =  U0  | 


U"  =  2U0 

x-^ 


(11.15) 


Thus  U*  <  0  and  U**  >  0  for  all  x  so  that  the  necessary  condition 
equation  (11.7a)  for  avoiding  separation  is  satisfied.  However,  calcula¬ 
tion  of  the  dimensionless  number  a  according  to  equation  (11. 5)  gives 


a  =  2 


(11.16) 


The  sufficient  condition  for  avoiding  separation,  a  >  11  according  to 
equation  (11.7)  is  therefore  violated.  For  the  divergent  channel  with 
plane  walls  separation  therefore  occurs  for  any  included  angle.  This 
example  shows  especially  clearly  the  low  ability  of  the  laminar  flow  to 
overcome  a  pressure  increase  without  separation.  According  to  a  calcula¬ 
tion  of  Pohlhausen  (reference  24)  the  separation  point  lies  at  (x/a)A  =  1.213 
and  thus  is  independent  of  the  included  angle  a. 


b.  Various  Technical  Arrangements  for  Avoiding  Separation 

It  is  a  favorable  circumstance  for  technical  applications  that  for 
higher  Reynolds  numbers  the  boundary  layer  does  not  remain  laminar  but 
becomes  turbulent.  The  turbulence  consists  of  an  irregular  mixing  motion. 
By  this  mixing  motion  momentum  is  continuously  transported  into  the  layers 
near  the  wall,  and  the  particles  retarded  at  the  wall  are  carried  out  into 
the  free  stream  and  thus  re— accelerated. 

Because  of  this  mechanism  the  turbulent  flow  is  able  to  withstand, 
without  separation,  considerably  higher  pressure  increases  than  the 
laminar  flow;  thus  the  pressure  increases  existing  in  technical  flows 
are  made  possible. 

A  few  technical  possibilities  for  avoiding  separation  will  be 
di scussed. 

1.  Blowing.  For  a  wing  profile  the  separation  of  the  boundary 
layer  for  large  angles  of  attach  (fig.  65)  can  be  prevented  by  blowing 
air  in  the  flow  direction  from  a  slot  directed  toward  the  rear.  The 
velocity  for  the  layer  near  the  surface  is  thus  increased  by  the  energy 


no 
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supplied  and  the  danger  of  separation  is  therefore  eliminated.  It  is 
true  that  in  the  practical  execution  not  much  is  gained,  because  of  the 
large  Jet  energy  required  for  any  considerable  improvement  of  the  flow. 

In  order  to  make  the  energy  output  small,  the  width  of  the  Jet  mast  be 
kept  small.  But  then  the  Jet,  soon  after  its  exit,  breaks  up  into 
vortices. 

2.  Another  possibility  of  avoiding  separation  is  the  arrangement 
of  a  slotted  wing  according  to  figure  66.  The  effect  depends  on  the 
boundary  layer  formed  on  the  slot  AB  being  carried  away  into  the  free 
stream,  before  it  separates,  by  the  flow  through  the  slot.  A  new 
boundary  layer  develops  at  C  which  is,  however,  at  first  very  thin 
and  reaches  D  without  separation. 

The  same  principle  is  used  for  the  Townend  ring  and  NACA  cowling 
(fig.  67). 

3*  Suction.  A  further  possibility  for  the  prevention  of  separation 

is  suction.  For  the  wing,  for  instance,  the  retarded  boundary— layer 

material  is-  sucked  off  into  the  interior  of  the  wing  through  one  or 

several  slots  (fig.  68).  The  point  of  suction  lies  slightly  ahead  of 

or  behind  the  separation  point  so  that  no  reversal  of  the  flow  can  occur. 

A  new  boundary  layer  which  at  first  is  very  thin  develops  behind  the 

suction  point  and  permits  the  pressure  to  increase  further.  In  thiB 

manner  one  can  overcome  considerably  larger  pressure  increases  and 

attain  higher  values  of  maximum  lift  for  the  wing.  Many  different  suction 

arrangements  for  increasing  maximum  lift  have  been  investigated  by 

0.  Schrenk  (reference  28).  Values  for  c_  of  3  to  4  were  obtained. 

'  a  max 


c.  Theory  of  the  Boundary  Layer  with  Suction 

Suction  is  a  very  effective  means  for  influencing  the  friction 
layer  on  a  body  immersed  in  a  flow  and  particularly  for  avoiding 
separation.  This  was  pointed  out  for  th9  first  time  in  1904  by  L.  Prandtl 
in  his  fundamental  work  on  the  boundary  layer. 

Another  possibility  of  application  of  suction,  recognized  only 
recently,  is  to  keep  the  friction  layer  laminar.  Here  the  boundary  layer 
is,  by  suction,  kept  so  thin  that  transition  to  the  turbulent  state  of 
flow  is  avoided.  The  surface  friction  drag  is  thereby  reduced.  Experimental 
investigations  of  this  effect  were  carried  out  by  Ackeret  (reference  39). 

The  laminar  friction  layer  with  suction  can  also  be  subjected  to  a 
numerical  treatment  which  will  be  briefly  discussed.  The  following 
assumptions  are  made  for  th9  calculation: 
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1.  The  suction  is  introduced  into  the  calculation  through  the 

assumption  that  the  normal  velocity  at  the  wall  vQ(x)  is 
different  from  zero.  The  wall  is  therefore  assumed  to  he 
permeable.  A  continuous  distribution  of  the  suction  velocity 
vQ(x)  serves  the  purpose  of  numerical  treatment  best. 

2.  The  suction  quantities  are  so  small  that  only  the  partB  in  the 

immediate  neighborhood  of  the  wall  are  sucked  from  the 
boundary  layer.  This  leads  to  a  very  small  ratio  of  suction 
velocity  vQ(x)  to  free  stream  velocity  UD:  vD/uo  =  °'001 

to  0.01. 

3.  The  no— slip  condition  at  the  wall  u  =  0  is  retained  with 

suction,  likewise  the  expression  for  the  wall  shearing  stress 


The  equations  of  motion  for  the  boundary  layer  with  suction  therefore 
read 


bx  by  dx  ^2 


du  |  bv 
bx  by 


0 


> 


J 


(11.17) 


with  the  boundary  conditions 


y  =  0  u  =  0  v  =  vQ(x) 
y  =  oo  u  =  U 

vQ  <  0  signifies  suction;  vQ  >  0  blowing. 


> 


(11.18) 


As  in  chapter  X  b  the  momentum  theorem  is  again  applied  to  the 
boundary  layer  with  suction.  The  momentum  equation  for  the  boundary  layer 
with  suction  is  obtained  in  exactly  the  same  manner  as  in  chapter  X  b 
(compare  fig.  Uk)  provided  one  takes  into  consideration,  in  addition, 
that  the  normal  velocity  at  the  wall  is  different  from  zero.  In  chapter 
X  b  the  momentum  equation  was  derived  by  integration  of  the  equation  of 
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motion  for  the  x-direction  over  y  between  the  limits  y  =  0  and 
y  =  oo  (compare  equation  (10.28)).  One  imagines  exactly  the  same  calcula¬ 
tion  performed  for  the  boundary  layer  with  suction:  then  the  expression 
for  the  normal  velocity  at  the  distance  from  the  wall  y  =  h  is 
different,  compared  with  the  calculation  in  chapter  X  b.  The  normal 
velocity  now  becomes 


(11.19) 


The  remaining  calculation  is  exactly  the  same  as  in  chapter  X  b  and 
finally  yields  as  the  momentum  equation  for  the  boundary  layer  with 
suction. 


1°  =  U2  M  +  (2.3  +  6*)  U  —  —  v_  U 
p  dx  dx  0 


(11.20) 


The  newly  added  term  — v  U  (compared  with  equation  (10. 36))  gives  the 
loss  of  momentum  due  to  the  suction  at  the  wall. 


We  shall  now  treat  the  special  case  of  the  flat  plate  with  suction 
in  longitudinal  flow  (fig.  69)  (reference  29).  The  free-etream  velocity 
is  U0.  Equation  (11. 20)  then  becomes 


uQ2  —  -  vQU  =  Is 
0  dx  0  P 


(11.21) 


if  one  takes  the  law  for  the  laminar  wall  shearing  stress  into  consider¬ 
ation.  Furthermore,  the  assumption  is  made  that  the  suction  velocity 
(or  blowing  velocity)  — vQ  along  the  plate  is  constant.  In  this  case 
one  can  obtain  from  the  momentum  equation  (11.21),  by  the  following 
simple  calculation,  an  estimate  of  the  variation  of  the  momentum  thick¬ 
ness  along  the  plate.  One  puts 

(!)0  ■ e  r  (11-22) 


0  —  0  signifying  a  dimensionless  form  parameter  of  the  velocity  profile . 
It  may  be  assumed,  to  a  first  approximation,  that  0  varies  only  little 
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with  the  length  of  run  x;  accordingly,  3  will  be  considered  constant. 
Then  equation  (11.21)  may  be  written 


1  dd 

2  dx 


(11.23) 


with  the  initial  condition  d  =  0  for  x  =  0.  For  suction  (v  <  0)  one 
obtains  dd/dx  =  0  for 


d«*  =  3  ——  (suction) 
~vo 


(11.24) 


(that  is,  therefore,  th9  momentum  thickness  reaches,  after  a  certain 
approach  length,  a  constant  asymptotic  value  given  by  equation  (11.24)). 
Simultaneously,  displacement  thickness,  velocity  distribution,  and  all 
other  boundary  layer  coefficients  also  become  asymptotically  independent 
of  x. 


For  blowing  (vQ  >  0)  the  value  dd/dx  is,  according  to  equation 

( 11 . 23 ),  larger  than  zero  along  the  entire  plate;  that  is,  d(x)  increases 
with  the  length  of  run  x  without  limit  so  that  for  large  values  of  x, 
one  can  neglect  in  equation  (11.23)  "the  first  term  on  the  right  side  as 
compared  with  the  second.  One  obtains  therefore,  as  asymptotic  law, 

=  jj2  x  ("blowing)  (11.25) 


On  the  whole,  one  obtains  the  remarkable  result  that  for  the  flat  plate 
in  longitudinal  flow  with  constant  suction  or  blowing  velocity,  the 
boundary  layer  thickness  for  suction  becomes  constant  after  a  certain 
approach  length,  whereas  for  blowing,  it  increases  proportionally  to  the 
length  of  run  x.  In  between  lies  the  case  of  the  impermeable  wall 
where  the  boundary  layer  thickness  increases  with  \Jx. 

For  the  case  of  the'  laminar  boundary  layer  with  the  asymptotically 
constant  boundary  layer  thickness  it  is  also  possible  to  give  immediately 
an  exact  solution  of  the  differential  equations  of  the  boundary  layer  in 
a  surprisingly  simple  form.  In  this  case  du/dx  =  0,  hence  also, 
according  to  equation  (11.17),  dv/§y  s  0  and  therefore 


v(x,y)  =  vq  =  constant 


(11.26) 
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Hence  there  follows  from  equation  (11.17) 


v 

o 


du 

dy 


(11.27) 


and  from  it  the  solution  which  satisfies  the  boundary  conditions 
equation  (11.18) 


u(y)  =  uo^i  - 


(11.28) 


From  this  equation  results  the  displacement  thickness  of  the  asymptotic 
boundary  layer 


(11.29) 


the  momentum  thickness 


00 


(11.30) 


5* 

•0 

and  the  form  parameter  - 

^00 


2.  By  comparing  equation  (11.29)  with 


equation  (11.24)  one  finds  the  factor  0=1.  The  velocity  distribution 
of  the  asymptotic  boundary  layer  profile  according  to  equation  (11.28) 
is  plotted  in  figure  70  together  with  the  Blasius  solution  for  the 
impermeable  wall. 


Herewith  the  considerations  of  boundary  layer  with  suction  will 
be  concluded. 


CHAPTER  XII.  APPENDIX  TO  PART  I 
a.  Examples  of  the  Boundary  Layer  Calculation 
According  to  the  Pohlhausen— Holstein  Method 


For  the  integration  of  the  differential  equation  (10. 58)  it  is 
best  to  use  the  isocline  method.  It  is  expedient  to  calculate  with 
dimensionless  quantities.  The  arc  length  s  is  made  dimensionless  by 
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dividing  by  a  characteristic  length  of  the  body  immersed  in  the  flow, 
for  instance,  for  the  wing,  by  the  wing  chord  t.  The  variable  Z  =  32/V' 

UQ 

is  made  dimensionless  by  multiplying  by  — .  Thus  one  puts: 

t 


Z* 


2 

*  U 

_ o 

v  t 


i'w 

.t  /  V 


X* 


Hence  the  differential  equation  reads: 


> 


(12.1) 


^  =  F(k);  (12.2) 

o 

The  calculated  example  concerns  a  symmetrical  wing  profile  (J  015)  in 
symmetrical  approach  flow  (ca  =  0).  The  prescribed  potential-flow 
velocity  and  its  first  derivative  with  respect  to  the  arc  length  is 
given  in  tfible  6.  The  initial  values  for  the  integration  are  calculated, 
according  to  equation  (10. 60),  to  be,  for  the  present  case; 


ZQ*  =  0.00  11+9 


since  at  the  stagnation  point  d2U/ds2  =  0.  The  auxiliary  function  F(k) 
required  for  the  integration  is  given  in  figure  1+7— a  and  table  5*  The 
calculation  according  to  the  isocline  method  is  shown  in  figure  1+8.  Here 
the  curve  k  =  -0.1567  which  gives  the  separation  point  can  be  calculated 
according  to  the  relation: 


_t_  dU 
U0  ds 


0.1567) 


Jt_  dU 
U0  ds 


The  intersection  of  the  integral  curve  with  this  curve  gives  the  separa¬ 
tion  point.  As  a  result  of  the  integration  one  obtains  at  first  the  varia¬ 
tion  of  the  momentum  thickness.  By  means  of  the  function  5*/9  =  f(K) 


and 


llo 

U  ++ 


=  T2(k) 


given  in  table  5  one  can  also  calculate  the  displacement 
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thickness  and  the  shearing  stress.  The  result  of  the  calculation  Is 
compiled  in  table  6  and  given  in  figure  49.  Moreover,  the  velocity 
distribution  in  the  boundary  layer  can  be  seen  from  figure  56. 


Translated  by  Mary  L.  Mahler 
National  Advisory  Committee 
for  Aeronautics 
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TABLE  I.  -  THE  FUNCTION  <p  OF  THE  PLANE  STAGNATION  POINT 
FLOW  (ACCORDING  TO  HTKMENZ  (REFERENCE  10));  TO  FIGURE  l6 

*  *  ]j^  y;  y  =<p'(i) 


<P 

<P' 

rp” 

0 

0 

0 

1.23264 

0.1 

0.0060 

0.1183 

1.1328 

0.2 

0.0233 

0.2266 

1.0345 

0.3 

0.0510 

0.3252 

0.9386 

0.4 

0.0881 

0.4144 

0.8463 

0.5 

0.1336 

0.4946 

0.7583 

0.6 

0.1867 

0.5662 

0.6751 

0.7 

0.2466 

0.6298 

0.5973 

0.8 

0.3124 

0.6859 

0.5251 

0.9 

0.3835 

0.7350 

0.4586 

1.0 

0.4592 

0.7778 

0.3980 

1.1 

0.5389 

0.8H9 

0.3431 

1.2 

0.6220 

0.8467 

0.2937 

1.3 

0.7081 

0.8739 

0.2498 

1  4 

0.7966 

0.8968 

0.2109 

1.5 

0.8873 

O.9161 

0.1769 

1.6 

0.9798 

0.9324 

0.1473 

1.7 

1.0738 

0.9457 

0.1218 

1.8 

1.1688 

O.9569 

0.0999 

1.9 

1.2650 

O.9659 

0.0814 

2.0 

1.3619 

0.9732 

O.O658 

2.1 

1.4596 

0.9792 

0.0528 

2.2 

1.5577 

0.9841 

0.0420 

2.3 

1.6563 

0.9876 

0.0332 

2.4 

1.7552 

O.9905 

0.0260 

2.5 

1.8543 

0.9928 

0.0202 

2.6 

1*9537 

O.9946 

0.0156 

2.7 

2.P533 

0.9960 

0.0119 

2.8 

2.1529 

0.9971 

0.0091 

2.9 

2.2528 

0.9979 

0.0068 

3.0 

2.3525 

0.9985 

0.0051 

3-1 

2.4523 

0.9988 

O.OO36 

3.2 

2.5522 

0.9992 

0.0027 

3-3 

2.6521 

O.9994 

0.0023 

3-4 

2,7521 

0.9996 

0.0019 

3-5 

2.8520 

0.9997 

0.0014 

3.6 

2.9520 

O.9998 

0.0010 

3-7 

3.0519 

0.9999 

0.0008 

3-8 

3.1518 

0.9999 

0.0004 

3-9 

3.2518 

0.9999 

0.0003 

4.0 

3.3518 

1.0000 

0.0002 

4.1 

3.4518 

1,0000 

0.0001 

4.2 

3.5518 

1.0000 

0.0001 

4.3 

3.6518 

1.0000 

0 
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TABLE  II.  -  THE  FUNCTION  f  OF  THE  BOUNDARY  LAYER  ON  THE  FLAT  PLATE 
IN  LONGITUDINAL  FLOW  (ACCORDING  TO  BLA3IUS  (REFERENCE  3));  TO 

FIGURE  JO  n  --  vj  Ha;  JL  =  f ■  ( « ) 

V  VI  Un 


n 

L  r  _ 

f» 

f" 

0 

0 

0 

0. 33206 

0.2 

0.00664 

0.06641 

0.33199 

0.4 

0.02656 

0.13277 

0.33147 

0.6 

0.05974 

0.19894 

0.33008 

0.8 

0.10611 

0.26471 

0.32739 

1.0 

0.16557 

0.32979 

0.32301 

1.2 

0.23795 

0. 3Q378 

0.31659 

1.4 

0. 322Q8 

0.45627 

0.30787 

1.6 

O.42032 

0.51676 

0.29917 

1.8 

0.52982 

0.57477 

0.28293 

2.0 

1  O.65OO3 

0.62977 

0.26675 

2.2 

0.78120 

0.68132 

0.24835 

2.4 

O.9223O 

0.72899 

0.22809 

2.6 

1.07252 

0.77246 

0.20646 

2.8 

1.23099 

0.81152 

0.18401 

3.0 

1.39682 

0.84605 

O.I6136 

3.2  ; 

1.56911 

0.87609 

0.13913 

3.4 

1.74696 

0.90177 

0.11788 

3.6 

1.9295^ 

0.92333 

0.09809 

3.3 

2.11605 

0.94112 

0.08013 

4.0 

2. 30^76 

0.9r'r'R2 

0.06424 

4.2 

2.49806 

O.96696 

0.0^052 

4.4 

2.69238 

0.97c87 

O.O3897 

4 . 6 

2.88326 

0.98269 

0.02948 

4.8 

3.08534 

0.98779 

0.02187 

c.O 

3.28  329 

0.991c<5 

0.01591 

c.2 

3.48189 

0.9942^ 

O.OH34 

5.4 

3 . 68094 

O.99616 

0.00793 

^.6 

3.88031 

0.99748 

0.00543 

c.8 

4.07990 

0.99838 

O.OO365 

6.0 

4 . 27964 

O.99898 

0 . 00240 

6.2 

4.47948 

0.99937 

0.00155 

6.4 

4.67938 

0.99961 

0.00098 

6.6 

4.87931 

0.99977 

0.00061 

6.8 

5.07928 

0.99987 

0.00037 

7.0 

5.27926 

0.99992 

0.00022 

7.2 

5.47928 

0.99996 

0.00013 

7.4 

5.67924 

0.99998 

0.00007 

7.6 

5.87924 

0.99999 

0 . 00004 

7.8 

6.07923 

1.00000 

0.00002 

8.0 

6.27923 

1 . 00000 

0.00001 

8.2 

6.47923 

1.00000 

0.00001 

8.4 

6.67923 

1 . 00000 

0 

8.6 

6.87923 

1.00000 

0 

8.8 

7.07923 

1 . 00000 
_ - _ I 

0 
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TABLE  III.  -  TEE  FUNCTIONS 


AND 


h  OF  THE  BOUNDARY  LAYER 
5 


ON  THE  CYLINDER  (SYMMETRICAL  CASE)  ACCORDING  TO  HOW ART fl 


(REFERENCE  lr' )  AND  FROSSLING  (REFERENCE  161. 


5 

f3 

f' 

3 

f" 

_ 3 _ 

0 

0 

0 

0.7244 

0.1 

0.0035 

0.0675 

0.6249 

0.2 

0.0132 

0.1251 

0.5286 

0.3 

0.0282 

0.1734 

0.4375 

0.4 

0.0476 

0.2129 

0.3539 

0.5 

0.0705 

0.2444 

0.2780 

0.6 

0.0962 

0.2688 

0.2112 

0.7 

0.1240 

0.2869 

0.1530 

0.8 

0.1534 

0.2997 

0.1037 

0.9 

O.1838 

0.3080 

0.0626 

1.0 

0.2149 

0.3125 

0.0292 

1.1 

0.2462 

0.3140 

0.0028 

1.2 

0.2776 

0.3132 

-0.0173 

1-3 

O.3O88 

0.3107 

-O.O32O 

1.4 

0.3397 

O.3070 

-0.0420 

1.5 

0.3702 

O.3025 

-0.0482 

1.6 

0.4002 

0.2947 

-O.O513 

1.7 

0.4297 

0.2923 

-0.0518 

1.8 

0.4587 

0.2871 

-0.0506 

1.9 

0.4871 

0.2822 

-0.0480 

2.0 

0.5151 

0.2775 

-0.0444 

2.1 

0.5426 

0.2733 

-0.0402 

2.2 

O.5698 

0.2695 

-O.O358 

2.3 

0.5966 

0.2662 

-0.0314 

2.4 

0.6230 

O.2632 

-0.0271 

2.5 

0.6492 

0.2607 

-0.0230 

2.6 

0.6752 

0.2586 

-0.0194 

2.7 

0.7010 

0.2568 

-0.0160 

2.8 

0.7266 

0.2554 

-O.OI3I 

2.9 

0.7520 

0.2542 

-0.0106 

3.0 

0.7774 

0.2533 

-O.OO85 

3.1 

0.8027 

0.2525 

-0.0067 

3-2 

0.8279 

0.2519 

-O.OO52 

3-3 

0.8531 

0.2515 

-0.0041 

3.4 

0.8782 

0.2511 

-0.0032 

3*5 

0.9033 

0.2508 

-0.0024 

3.6 

0.9284 

0.2506 

-0.0019 

3-7 

0.9534 

0.2504 

-0.0014 

3.8 

0.9785 

0.2503 

-0.0011 

3*9 

1.0035 

0.2502 

-0.0008 

4.0 

1.0285 

0.2502 

-O.OOO6 

4.1 

1.0535 

0.2501 

-0.0004 

4.2 

1.0785 

0.2501 

-0.0003 

4.3 

1.1035 

0.2500 

-0.0002 

4.4 

1.1285 

0.2500 

-0.0001 

5 

g5 

g; 

a” 

g5 

0 

0 

0 

0.6348  ^ 

0.2 

0.0114 

0.1072 

0.4402 

0.4 

0.0405 

0.1778 

0.2717 

0.6 

0.0806 

0.2184 

0.1408 

0.8 

0.1264 

O.2367 

0.0483 

1.0 

0.1742 

0.2399 

-0.0106 

Esm 

0.2218 

0.2342 

-0.0431 

Hi 

0.2676 

0.2239 

-O.O567 

IKsHh 

0.3112 

0.2123 

-O.O58O 

1.8 

0.3526 

0.2012 

-0.0522 

2.0 

0.3918 

0.1916 

—0.0432 

2.2 

0.4293 

0.1839 

-0.0335 

2.4 

0.4655 

0.1781 

-0.0245 

2.6 

O.5OO7 

0.1740 

-0.0171 

2.8 

0.5352 

0.1712 

-0.0114 

3.0 

O.5692 

0.1694 

-0.0072 

op 

O.6O3O 

0.1682 

-0.0043 

3*4 

O.6365 

0.1676 

-0.0026 

3.6 

0.6700 

0.1672 

-0.0015 

3.8 

0.7034 

0.1669 

-0.0010 

4.0 

0.7368 

0.1668 

-0.0004 

4.2 

0.7701 

0.1667 

-0.0001 

4.4 

0.8035 

0.1667 

-0.0001 

9 

h5 

b5 

0 

0 

0 

0.1192 

0.2 

0.0017 

0.0141 

0.0249 

0.4 

0.0045 

0.0117 

-0.0436 

0.6 

-0.0010 

-O.O783 

0.8 

-0.0176 

-O.O833 

1.0 

-0.0012 

-0.0330 

-0.0680 

1.2 

-0.0090 

-0.0441 

-0.0423 

1.4 

-0.0185 

-0.0498 

-0.0149 

1.6 

-0.0286 

-0.0503 

+0.0088 

1.8 

-0.0384 

-0.0468 

0.0256 

2.0 

-0.0472 

-0.0406 

0.0351 

2.2 

-0.0K46 

-O.O33I 

0.0380 

2.4 

-0.0604 

-0.0257 

O.O36I 

2.6 

-0.0649 

-0.0189 

0.0312 

2.8 

-O.O68I 

-0.0133 

0.0249 

3.0 

-O.O703 

-0.0089 

0.0187 

3.2 

-0.0717 

-O.OO58 

0.0132 

3-4 

-0.0726 

-0.0036 

O.OO89 

3.6 

-0.0732 

-0.0021 

0.0057 

3.8 

-0.0735 

-0.0012 

O.OO36 

4.0 

-0.0737 

-0.0006 

0.0022 

4.2 

-0.0738 

-0.0003 

0.0012 

4.4 

-0.0738 

-0.0001 

0.0007 
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TABLE  IV.  -  THE  FUNCTIONS  F(y/5p)  and  G(y/5p)  FOR  THE  VELOCITY 
DISTRIBUTION  IN  THE  BOUNDARY  LAYER  ACCORDING  TO  POHLHAUSEN 
(REFERENCE  24)  AND  HOWARTH  (REFERENCE  15) 


y 

Sp 

F 

G 

0 

0 

0 

0.1 

0.1981 

0.01215 

0.2 

O.3856 

0.01725 

0.3 

0.5541 

0.01715 

0.4 

O.6976 

0.0144 

0.5 

0.8125 

0.0104 

0.6 

0.8976 

0.0064 

0.7 

0.95^1 

0.00315 

0.8 

O.9856 

0.00105 

0.9 

O.998I 

0.00015 

1.0 

1 

0 

124 


NACA  TM  No.  1217 


TABLE  V.  -  AUXILIARY  FUNCTIONS  FOR  THE  BOUNDARY  LAYER  CALCULATION 
ACCORDING  TO  HOLSTEIN  (REFERENCE  2e  ) 


0.0885 

0.0920 

0.0941 

0.0948 

0.0941 

0.0920 

0.0382 

0.0831 

0.0820 

0.0807 

0.0794 

0.0780 


0.0767 

0.0760 

0.0752 

0.0744 

0.0737 

0.0729 

0.0721 

0.0713 

0.0706 

0.0697 

0.0689 

0.0599 

0.0497 

O.O385 

0.0264 

0.0135 


-0.0140 

-0.0284 

-0.0429 

-0.0575 

-0.0720 

-0.0862 

-0.0999 

-O.U30 

-0.1255 

-O.1369 

-0.1474 


-0.1567 


-0.1648 

-O.1715 

-0.1767 


F(k) 

=  f 

-0.0657 

2.279 

-0.0814 

2.262 

-0.0913 

2.253 

-0.0946 

2.250 

-0.0911 

2.253 

-0.0806 

2.260 

-O.O608 

2.273 

-0.0332 

2.289 

-0.0271 

2.293 

-0.0203 

2.297 

-O.OI32 

2.301 

-0.0051 

2.305 

0 

2.308 

0.0021 

2.309 

0.0061 

2.312 

0.0102 

2.314 

0.0144 

2.316 

0.0186 

2.318 

0.0230 

2.321 

0.0274 

2.323 

0.0319 

2.326 

O.O365 

2.328 

C.0412 

2.331 

0  0459 

2.333 

0.0978 

2.361 

0.1579 

2.392 

0.2255 

2.427 

0.3000 

2.466 

O.382O 

2.508 

0.4698 

2.554 

2.604 

2.658 

HiiMiDlvK*  bzmsms 

2.716 

0.8698 

2.779 

0.9780 

2.847 

1.0853 

2.921 

1.1981 

2.999 

1.3078 

3.084 

1.4173 

3-177 

1.5231 

3-276 

1.6251 

3.383 

1.7237 

3.5OO 

1.8159 

3.627 

1.9020 

3.765 

1.9821 

3.920 

f2(K)  =  ^ 
*■  U  u 


.331 

.330 

0.330 

0.329 

O.328 

0.327 

O.326 

0.325 

0.324 

O.322 

O.32I 

0.310 

0.297 

O.283 

0.268 

0.252 


-0.019 

-0.037 

-0.054 


TABLE  71.  -  EXAMPLE  FOR  THE  BOUNDARY  LAYER  CALCULATION  ACCORDING 
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Figure  1.-  Simple  shear  flow. 


xyz  dx 


Figure  3.-  The  general  stress  tensor. 


NACA  TM  No.  1217 


Figure  4.-  The  shearing  stress  (to  fig.  3). 


y 


y 


Figure  6.-  Pure  angular  deformation  (e  >  0). 
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Figure  13.-  Velocity  distribution  on  an  oscillating  surface. 


Figure  16.-  The  velocity  profile  of  the  plane  stagnation  point  flow 


Figure  17.-  The  convergent  and  divergent  channel. 


Figure  18.-  Velocity  distribution  in  the  convergent  channel. 
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Figure  19.- 


Figure  20.- 


Uo 


Velocity  distribution  in  the  divergent  channel. 


Separation  in  the  divergent  channel. 


Figure  21.-  Viscous  flow  around  a  sphere. 
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Figure  22.-  Streamline  pattern  of  the  viscous  flow  around  a  sphere 

(according  to  Stokes). 


Figure  23.-  Streamline  pattern  of  the  viscous  flow  around  a  sphere 

(according  to  Oseen). 


Figure  24.-  Concerning  Prandtl’s  boundary-layer  equation.  (Boundary - 

layer  thickness  5  magnified.) 
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Figure  26.-  Velocity  distribution  in  the  boundary  layer  for  pressure 

decrease  f  &  <  0 ). 


Figure  27.- 


Velocity  distribution  in  the  boundary  layer  for  pressure 
increase  >  0^. 


Figure  28.-  Concerning  the  calculation  of  the  friction  drag. 


138 


NACA  TM  No.  1217 


Figure  29.-  The  boundary  layer  on  the  flat  olate  in  longitudinal  flow. 


Figure  30.-  Velocity  distribution  u(x,y)  in  the  boundary  layer  on  the 

flat  plate  (according  to  Blasius). 


Figure  31.-  The  transverse  velocity  v(x,y)  in  the  boundary  layer  on 

the  flat  plate. 


Figure  32.- 


The  boundary  layer  on  a  cylindrical  body  of  arbitrary 
cross  section  (symmetrical  case). 


0 
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Figure  35.-  Velocity  distribution  of  the  potential 
flow  for  a  wing  profile. 


Figure  36.-  Concerning  the  calculation  of  the  friction  layer  on  the 

circular  cylinder. 


Figure  37.-  Wake  flow  behind  the  flat  plate  in  longitudinal  flow. 
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Figure  38.-  Concerning  application  of  the  momentum  theorem  for  the 

flat  plate  in  longitudinal  flow. 


Figure  39.  -  Asymptotic  velocity  distribution  in  the  wake  behind  the 

flat  plate  in  longitudinal  flow. 
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Figure  42.-  Application  of  the  momentum  theorem  for  the  flat  plate  in 

longitudinal  flow. 


Figure  43.-  Velocity  distribution  in  the  boundary  layer  on  the  flat  plate 

in  longitudinal  flow. 

(a)  Linear  approximation. 

(b)  Cubic  approximation  for  the  velocity  profile. 


Figure  44.-  Application  of  the  momentum  theorem  to  the  boundary  layer 

with  pressure  gradient. 
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Figure  45.-  The  universal  functions  F(y/&p)  and  G(y/&p)  for  the 
velocity  distribution  in  the  boundary  layer  according  to  Pohlhausen. 


Figure  46.- 


The  one -parameter  family  of  velocity  profiles  according 
to  Pohlhausen. 


igure  47 
accorc 


F: 


(a).-  Auxiliary  functions  of  the  boundary  layer  calculation 
ding  to  Holstein  (cf.  table  b);  \  and  F(  * )  against  K. 


Stagnation  point 


48 
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Figure  47(b).-  Auxiliary  functions  for  the  boundary  layer  calculation 
according  to  Holstein  (cf.  table  5);  f|(i«)  and  f9(*  )  against  k. 
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'  t 


•■'ij'ure  Intei'ratiun  of  the  * ii f !V* -nt i -j  1  of  tl.f*  tx .ut i- i':  ry  loy  r 

acrorlinp  to  I  ohlhanoon  ;ui.i  iioJo.loiu  (orofilf*  •  01:  ;  o:)  •)). 
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Joukowsky  profile  d/t  =  0.15 

Figure  49.-  Result  of  the  boundary-layer  calculation  for  the  example 
according  to  figure  48  (profile  J  015;  ca  =  0). 
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U  m  Ui  x 


b) 


Plate  profile 

v-ue 


y/&* 


Figure  50.-  Comparison  of  the  approximate  calculation  according  to 

Pohlhausen  with  the  exact  solution. 
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Figure  51.-  Potential -theoretical  velocity  distrioution  on  the  elliptic 
cylinders  with  axis  ratio  a^/b^  =  1,  2,  4,  8  for  flow  parallel  to  the 

major  axis  (A  =  laminar  separation  point),  t'  =  half  the  circumference. 


Semicircumference 


Figure  53.-  Potential  flow  and  laminar  friction  layer  on  the  elliptic 

cylinder  of  axis  ratio  aj/b,  =  4. 
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Figure  54.-  Potential -theoretical  velocity  distribution  for  the  Joukowsky 
profile  J  015  for  ca  =  0;  0.25;  0.50;  1.0. 
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Figure  55.-  Result  of  the  boundary -layer  calculation  for  the  joukowsky 
profile  J  015  (t'  =  half  the  profile  perimeter). 


Suction 


Drofile 


JoukowsKy  profile 
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Figure  58.-  Joukowsky  profiles;  thickness  d/t  =  0  to  0.25;  camber 

f/t  =  0  to  0.08. 


Figure  59.-  Position  of  the  laminar  separation  point  for  the  Joukowskiyt ^ 
profiles  of  figure  58;  suction  side. 


SZ'O-  0, J 


Figure  60.  -  Position  of  the  laminar  separation  point  for  the  joukowsk  y 

profiles  of  figure  58;  pressure  side. 
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Figure  61.-  Pressure  distribution  and  separation  on  a  wing. 


Figure  62.-  Boundary  layer  with  laminar  separation  avoided. 

Ufxi 


Figure  63.-  Potential  flow  with  separation:  U;  <  0;  U"  <  0 
sometimes  without  separation:  U'  <0;  U"  >  0. 
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Figure  64.-  Divergent  channel. 


Figure  6b.-  Prevention  of  separation  on  wing  by  blowing. 


Figure  66.-  Prevention  of  separation  by  a  slotted  wing. 


Figure  67.-  NACA  cowling  for  prevention  of  separation. 
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Figure  68.-  Prevention  of  separation  on  wing  by  suction. 


Figure  69.-  Flat  plate  in  longitudinal  flow  with  suction. 


Figure  70.  -  Asymptotic  velocity  profile  on  flat  plate  in  longitudinal 
flow  with  suction  (I)  6*  =  v/-vQ  (II)  5  *  =  1.73  . 
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(Abstract  on  Beveree  Side)  ;  (Abstract  on  BeverBe  Side) 


NACA  TM  No.  121 V 


